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ABSTRACT 


This paper mathematically defines the course of the Wolf sunspot numbers during 
the sixteen completed cycles since 1755. The usual periodic or ‘‘superposition”’ point 
of view is rejected in favor of the ‘‘outburst’’ hypothesis of Halm and Waldmeier, which 
is here independently confirmed in detail. Each successive cycle or “‘outburst’’ is re- 
garded as a fresh phenomenon, following, however, a standard pattern. 

As a new addition to the ‘Wolf’? sunspot data published by W. Brunner in the 
Zurich statistics, the zero, first, and second moments for each cycle are calculated about 
the times of minimum, and all the observed data are employed to show that the em- 
pirical form R = F6¢e~9 represents, to a good first approximation, the course of spot 
numbers during a single outburst or cycle (of, on the average, eleven years). Here R 
is the Wolf number, @ is the interval in years after the time of beginning, s, of the out- 
burst, and F, a, and } are parameters which differ for different cycles but remain con- 
stant for a single cycle. Values of these four parameters, s, a, b, and F, are tabulated 
for each cycle, and the goodness of fit is discussed. 

There are small random discrepancies, year by year, averaging no more than half-a- 
dozen Wolf numbers, between observed and computed spot numbers; and an important 
systematic discrepancy which is particularly evident when the interval between the 
time of maximum and of centroid is examined. According to the form suggested, the 
centroid should always come after the maximum, but actually this is the case only for 
cycles of higher maximum. For the group of cycles of low maximum the observed cen- 
troid precedes the maximum; the interval between the two shows a strong correlation 
with the height of maximum. This means that in a second approximation the suggested 
formula for R will require modification. 

The computed parameters are examined for possible mutual dependence, and the 
indication is that # and 6 can be taken as rough functions of a. The even-numbered 
cycles have smaller values of a than any odd-numbered cycle. In view of the strong 
family resemblance shown by this study among apparently rather dissimilar cycles, 
the conventional periodogram analysis of spot numbers definitely must be abandoned, 
because such resemblance would be too improbable as the result of the superposition 
of independent periodic factors. Long-range prediction of spots will remain impracti- 
cable, but short-range prediction will be facilitated by the methods outlined in this 
paper. 

* The extensive contributions to this study by Mr. Panofsky represent his Senior 
thesis as a student of astronomy, toward his A.B., Princeton, 1938. 
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I. THE SUPERPOSITION HYPOTHESIS VERSUS 
THE OUTBURST HYPOTHESIS 

Many years ago Halm’ briefly pointed out that Sporer’s well- 
known rule of sunspot latitudes? suggests that each successive spot- 
cycle represents a new outburst independent of previous cycles. The 
appearance, near solar latitudes 30° north and south, of the first spots 
in each successive new outburst, followed by the gradual creeping 
toward the sun’s equator of the two zones of greatest activity, is an 
indication that each spot-cycle represents a fresh recurrence of a 
phenomenon which always follows a standard pattern of develop- 
ment. Each rise in the number of spots from minimum to maximum, 
together with the subsequent fading-back to minimum, constitutes 
a single outburst. Halm’s suggestion was less than conclusive, how- 
ever, and several careful analyses of sunspot numbers have held to 
the older superposition hypothesis, which interprets the rough 
eleven-year period as due to the combination of several continuously 
acting harmonic or pseudo-harmonic variations. These laborious 
analyses without exception have failed to predict the future rise 
and fall of spots. 

In an important paper, Waldmeier* has brought support to the 
outburst hypothesis. He suggested that the course of sunspot num- 
bers during each outburst, from minimum to minimum, could be 
represented by one particular curve of a family of curves. He gave 
the graphical, but not the’mathematical, form of sample curves of 
the family. The present paper suggests a simple mathematical form 
(eqs. [1] and [2]) and applies statistical methods of curve-fitting to 
prove that this form approximately represents all the annual sunspot 
numbers. 

The outburst hypothesis, therefore, becomes so definitely substan- 
tiated that the conventional super position approach must be abandoned. 

The analysis here presented deals with the annual Wolf numbers® 

1 A.N., 156, 44, 1901. 

2G. Abetti, Handbuch der Astrophysik, 4, 96, 1929. 

3H. Kimura, M.N., 73, 543, 1913; S. Oppenheim, A.N., 232, 369, 1928; H. Turner, 
M.N., 73, 549 and 714, 1913. But see Yule, ‘A Method of Investigating Periodicities, 
etc.” Phil. Trans., A, 226, 267, 1927. 

4 Astr. Mitt., Ziirich, 133, 105-30, 1935. (A good bibliography is included.) 


5’ Handbuch der Astrophysik, 4, 95 and 99, 1929. (These are the well-known spot 
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over the available years 1755-1937. These are reprinted in Table 1, 
in asomewhat novel form. 

During these years there were sixteen successive complete spot- 
cycles or outbursts from minimum to maximum and down again to 
minimum. It is shown in Sections V and VI that a fit, rough to 


TABLE 1* 


WOLF SUNSPOT NUMBERS 














: Spot NUMBERS FOR YEARS AFTER MINIMUM 
CycLe} MINI 
NuM-| MUM | —— _ —. | ) 
BER | YEAR | | 
| } oO I 2 3 5 6 7 'o] 9 | fe) II 7s. 1 Es 
2 
I 1755] 9.0/10.2/32.4) 47.0) 54.0) 02.9/55.9/01. 2/45 .1130.4/20.9 
2 1766|11 .4)37.8)69.8)106.1|/100.8) 81 .6/66.5/34.8)30.6 
oz 1775| 7-O}19.8/92.5|154.4/125.9| 84.8/68.1/38.5)22.8 a 
4 1784/10. 2/24. 1/82.9}132.0/130.9/118. 1/89 .9|66.6/60.0/46.9/41 .0}21.3/16.0/6.4 
5 1795} 4.1 ©. S134. 51 344-0) 45-0) 43.3147. 5142. 2/2 2/86. 2) GS. 1} 2.5 
6 1510) 0.0 1.4) 5.0] 12.2] 13.9) 35.4/45.5/41.1/30.4/23.Q]/15.7} 0.0] 4.0 
7 1523} 1.5) &.5|16.6) 36.3) 49.7] 62.5/67 o|71 0147 .5|27.5 
8 1533} 8.5}13.2/56.9/121 . 5/138. 3/103 . 2/85 . 8/63 . 2/36. 8)24. 2 
9 1543/10. 7|/15.0)40 I} 61.5) 98.51124. 3/95 .9/66. 5/64. 5/54. 2/39.0|20.6} 6.7 
IOu«:< 1550} 4.3}22.8154.8) 93.8) 95.7] 77.2159.1|44.0/47.0/30. 5/16. 3 
11 | 1867} 7 - 3137 -3}73-9/139- 1/111. 2/101. 7/66. 3/44. 7/17. 1/11. 3/12.3 
| ‘io - > = - > > ele% 919L 2 2 
12 1878] 3.4} 6.0/32.3] 54.3] 59.7] 63.7/63.5|52.2/25.4|13.1] 6.8 
13...| 1889} 6.3] 7.1/35.6) 73.0] 84.9] 78.0/64.0/41 .8)26.2/26.7/12.1| 9.5 
14...| IQOI] 2.7} 5.0}/24.4] 42.0] 63.5) 53.8/62.0/48.5/43.9/18.6) 5.7) 3.6 
15 | 1913] 1.4] 9.0147 4| 57.1|103.9| 80.6/63.6/37.6/26.1|14.2 
10...| 1923) 5 .8)10.7/44.3}) 63.9) 69.0) 77.8/65.0/35.7)21.2/11.1 
17 | 1933] 5-7| 8.7/36.1| 80.4]114.4 
| 


| | 





* E.g., the spot number for the year 1850 =1843 +7 is 66.5; and in this paper this number is attributed 
to the data 1850.5, since it is the average number for a calendar year. Numbers for 1749.5 through 1754.5 
are, respectively, 80.9, 83.4, 47-7, 47-8, 30.7, 12.2. 


excellent, for each of these sixteen outbursts can be obtained by 
suitable choices of four disposable parameters, F, a, 6, and s, in the 
family of curves, 

R = F(r — s)%e7*'-9) , (1) 
where R is the Wolf number for the varying year r. During a single 
cycle F, a, 6, and s are to be taken as constant. Here s is the time of 
the start of the new outburst; F is a scale factor; a is a pure number; 
the parameter b expresses the amount of damping and has the di- 
mensions of inverse time; 7 runs with the time from s to plus infinity. 
Forr < s, Ris taken as zero. 


numbers of Wolf, Wolfer, and Brunner. Values after 1927 are direct from the Astr. 
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II. RELEVANT PROPERTIES OF THE SUGGESTED FAMILY OF CURVES 
Writing 6 = r — s, the time in years since the start of the outburst, 
R= F@*c. (2) 
This is Pearson’s Type III distribution ;° but for present purposes the 
form must be examined afresh. The following results are readily 
5 - 


verified. 
a) The time at which maximum occurs is 


b) The spot number V at maximum is 


Vy = Fate~* 
—_ h ‘ \4) 
c) The area or zero moment for the outburst is 
ite FY(a+ 1) 
M, mS f F6%e-*°dé@ = pata ’ (5) 
> ) 


where I’ (a + 1) is the ordinary gamma function.’ 
d) The first moment with respect to the arbitrary time ¢ is 


~” . ‘(a+ 2 
M, -{ (7 — )Rdr = EF? + ts — 0). 
oO ) z 


~~ 


In the integration R is taken as zero for values of r < s; we write 
r—th=O0+(s—2). 


The centroid of the outburst comes at time 


— M,_ a+! 
u Re ae ee b? (6) 


since (a + 2) = (a +1) (a +1). 
6 Phil. Trans., A, 186, 373, 1895. 
7 Tabulated, e.g., by J. Brownlee to seven decimal places in Tracts for Computers, 


No. 9, Cambridge University Press, 1923. 
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ec) The second moment with respect to ¢ is 


x ay sf 2) 
M, = f (r — t)?Rdr = ees + 2M,(s — t) + Mi(s — 2). 


Thus, the invariant o? (independent of s — ¢) is 


= We _(# Pos (7) 
"Yar Ce 


This is a measure of the dispersion of the sunspot numbers during 
an outburst. (For a frequency curve o? so defined comes out the 
square of the standard deviation.) 

f) Another useful relation is 


ae lad a 


where 


g(a) = ————.. (9) 


Numerical values of (a), af(a), and (a + 1)¢(a) have been com- 
puted in Table 4. 

The flexibility of the form (1) or (2) is worth emphasizing. When 
y is any positive exponent, 


+ — FryAyao—yb0 
Rv = Fvgr%, yb , 


but, writing PF’ = F’, ya = a’, yb = Db’, this is just the original form 
(2). Accordingly, the rather successful agreement with observation 
described below establishes little as to the physical meaning of the scale 
of Wolf numbers. 


III. THE OBSERVATIONAL CRITERIA 
To fit a particular curve of the family (1) or (2) to any one of the 
sixteen available outbursts since 1749, a suitable value of each of the 
four adjustable constants, F, a, 6, and s, must be determined. In 


* 
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the Zurich statistics* three quantities especially applicable to this 
purpose are published for each cycle, namely: 

The exact time of minimum, which wil! be designated by u 

The time of maximum, designated here by v 

The sunspot number, V, at maximum (the tabulated “maximum 
smoothed monthly number’). 

In addition, the area M, can be readily approximated by summing 
the published annual sunspot numbers (of Table 1) between suc- 


TABLE 2 
OBSERVED DATA 




















ri roy 
Cycle | | 
i s ie: u 1 w Mo Uy M2 |M:1/Mol| 2? 
Number | 
I 1754.00} 55.5 SS.2 | O15 1 °OL 33) 407 2052 |17,710] 5.05 5.603 
2 1766.35) 66.5 | 66.5 | 69.7 | 70.79] 537 2226 |11, 266) 4.14 3.81 
) | | ~ | | : 
3 1775-37) 75-5 | 75-5 | 78.4 | 79.72] O15 | 2502 12,124) 4.07 | 3.17 
4 1784.53) 84.5 | 84.7 | 88.1 | 90.05) 843 | 4555 131,184) 5.40 | 7.19 
> ) | ss | ta ) ee 5 és » 
5 1797-59) 95.5 | 95.3 | 05.2 | 04.00) 277 E462 | 9,227) 5.35 | 4-06 
8 | | >| ed | | 
6 IS10.09| 10.5 [0.0 | 10:4 | 17.30] 230 I5QI |11,550} 6.73 | 4.95 
7 1822.65] 23.5 | 23-3 | 20.9 | 20.35] 302 | 2233 [14,357] 5-70 | 4.19 
8 1833.74 33-5 | 33-9 | 37-2 | 38.37] 653 | 3083 [17,158] 4.72 | 3.97 
9 1843.18] 43.5 | 43.5 | 48.1 | 49 51] 694 | 4068 28,354] 5.86 | 6.50 
Sega |e RE | ae ainie ; al bse 
10 1055-75} 50.5 56.0 | 60.1 OI 59 547 2095 |10,375) 4.93 | 5.01 
11 {1867 03) 67.5 | 67.2 | 70.6 | 71.92} 620 | 2647 |13,920) 4.27 | 4.25 
12 1878.52) 78.5 | 78.9 | 83.9 | 83.78} 382 | 1960 11,757| 5.23) 1 4.44 
13 11889.30} 89.5 | 89.6 | 94.1 | 94 83] 463 | 2395 |14,853] 5.18 | 5.24 
14 TOOL .37] 01.5 | OF. 7 | 00.4 | O7.1T] 373 | 20335 [12,873] 5.40: | 4.07 
15 1913.36] 13.5 | 13.6 | 17.6 | 18.30] 444 | 2105 |11,633] 4.74 | 3-71 
6... [1923.25] 23.5 | 23.6 | 28.4 | 28.35] 410 | 1928 |10,443] 4.70 | 3.38 
~ | 22 Or ab Bh. | 
17 33-5 33.5 | eee os! OE | 








cessive minima, and this we have done. Similarly, we have approxi- 
mated the first and second moments, M, and M,, about an arbitrary 
time /—taken as the time of the minimum near the beginning of an 
outburst (but see Sec. IV). 

Values of these seven quantities, together with values of two 
others, s and w, obtained as explained below, are presented in Table 
2. (Century numbers are omitted before the years, except for values 
of s.) The following specific example shows exactly how the tabu- 
lated values of M,, M,, and M, were obtained. Cycle number 2 
(Table 1) begins with the minimum of about 1766.5 and ends with 
the minimum of about 1775.5. The annual spot numbers for Cycle 


8 Astr. Mitt., Ziirich, 124, 77, 1930. 
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2—years 1766.5-1775.5, inclusive—are 11.4, 37.8, 69.8, 106.1, 100.8, 
81.6, 66.5, 34.8, 30.6, and 7.0. In computing the moments, the two 
minimum values are regarded as shared with the preceding and fol- 
lowing cycles, and therefore (as a rough allowance) are reduced to 
one-half in the computation. So M, is approximately 5.7 + 37.8 
+ 69.8, +....+ 30.6 + 3.5 = 537. The value of ¢ was uniformly 
taken as the year of the initial minimum annual number—here 
1775.5. Thus M, is (0)(5.7) + (1)(37.8) + (2)(69.8) +.... + (8) 
(30.6) + (9)(3.5) = 2226. Similarly M, is pies 7) + (1)(37.8) + 
27)(69.8) + .... + (87)(30.6) + (97)(3.5) = 11266. (All numeri- 
cal computations were carried out on a eieeael Monroe or a 
Marchant computing machine.) In each cycle the computations of 
the moments were carried in this manner from minimum to mini- 
mum annual numbers. (These rounded-off annual minima are not 
to be confused with the Zurich times u estimated for exact minima.) 

Empirical tests with idealized cases showed that the errors in- 
volved in taking the integrals for M,, M,, M, (of Sec. II) as equiva- 
lent to the summations (formed as just described) are negligible 
here. 

It is worth noting that the values of V used in Table 2 doubtless 
are not quite the best values that could have been adopted for the 
particular purposes of this paper. In place of using these Zurich 
maxima from smoothed monthly numbers, maxima could have been 
adopted from smoothed annual numbers. Such maxima would have 
been slightly smaller in many cycles. Their use would have reduced 
a little the systematic discrepancies discussed below—but by no 
means would they have eliminated them. In order to keep our study 
wholly ‘‘objective,’’ however, and since general conclusions would 
not have been affected, we employed the Zurich values. 


IV. OVERLAP AND END CORRECTION 
An outburst of the form (1), although damped out rapidly for 
large values of r — s, still retains a slight vitality even after the be- 
ginning of its successor outburst (Table 7). The spot numbers for 
several years after each minimum represent the sum of the effects 
of the new, growing cycle and the preceding, not yet quite negligible, 
cycle. The computation of the moments, however, ran only from 
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minimum to minimum. Since the theoretical moments represent 
summations to infinity, there are positive end corrections to be ap- 
plied to the values of M, and M, presented in Table 2. These cor- 
rections can be approximated in detail for particular cases (when a 
and 6 are known) by reference to tables for the incomplete gamma 
function.’ (In general, only a negligible correction is indicated for 
M, because the loss of area after the arbitrary cut-off at the second 
minimum is statistically balanced by the gain of area from the tail 


TABLE 3 


OBSERVED DATA: THE OBSERVATIONAL CRITERIA 





























Cycle Number | V Mo | V/Mo v-u v-s W-s w—- 
I | 86.5 467 0.185 6.3 6.90 6.73 Oo 7 
2 | 115.8 537 . 216 3-2 3-35 4.44 | +1.09 
3 | 158.5 615 257 2.9 3.03 4.35 | +1.32 
4 | 141.2 843 167 3-4 3-57 5.52 “a Os 
5 | 49.2 277 173 6.9 7.61 6.41 —1.20 
6 | 48.7 236 . 206 5.8 6.31 7.20 +o.98 
7 77 3902 183 6.6 7 ee 6.70 —0.55 
8 | 146.9 653 225 25% 3.40 4.03 1.17 
9 | 2348.0 694 19gO 4.6 4.92 6.33 +1.4I 
10 97.9 547 .179 4.1 4.35 5.83 +1.48 
i} | 140.5 620 227 3-4 3-57 4.89 +t. $2 
12 | 74.0 382 195 5.0 5.38 5.26 —O.12 
13 | 87.9 463 190 4.5 4.80 5.53 | +0.73 
14 | 64.2 373 -172 4.7 5.03 5-74 0.74 
15 | 105.4] 444 | 237 4.0 4.2 5 03 +0.79 
16 | 78.1 | 410 | 0.190 4.8 6.565 5.10 —0.05 
| ° | 





of the preceding outburst incorrectly included after the initial 
minimum. ) 

It turns out that the additive corrections so indicated for the 
various values of M, are so large as to render the values of M,, and 
consequently of 0”, in Table 2 rather unreliable (see Sec. VIII). 
Fortunately, only slight end corrections are indicated, by detailed 
study, for the values of M,. The uniform correction to M,/M, of 
+o.15 year, while somewhat arbitrary, seemed sound enough for 
use. The values in Table 2 are without this correction, however. 
Further discussion of these points is omitted for the sake of brevity. 

The Zurich data, of course, do not give the times of start, s, of 


9 Tables of the Incomplete V-Function. Karl Pearson, ed., London, 1922. 





























SUNSPOT VARIATIONS 393 


the outbursts but—disregarding random fluctuations in the spots— 
these must always a little precede the minima. Such precedence is 
because the observed spot-curve represents, as has been said, the 
sum of the effects of the fading and the fresh outbursts, and the 
new outburst at first is gaining less rapidly than the old one falls." 
Our examination of the individual cases with reference to appro- 
priate curves of the family (1) showed that an adequate empirical 
first approximation to s is 


$S=u—0.015(v — u)?, (10) 
where 2, it will be remembered, is the Zurich time of maximum, and 
wu of minimum. The values of s in Table 2 were so derived and, 
strictly speaking, are not “observed data’’; but the correction 
—o0.015 (v — w)?is always small. 

Again, detailed discussion is omitted for brevity’s sake. The 
small corrections described in the two immediately preceding para- 
graphs are undeniably in the right direction and are of about the 
right magnitude. Any reasonable quantitative changes in them 
would not make a significant difference. 


V. FITTING THE INDIVIDUAL SIXTEEN CYCLES; THE FIVE CASES 

Table 2 includes values of s computed by means of (10) from the 
Zurich values of « and v. It also presents values of w, defined as the 
times of the centroids of the outbursts, obtained by adding M,/M, 
to ¢ (eq. [6])—-with the inclusion in each case of the adopted end 
correction, +0.15 year, to the tabulated M,/M.,. 

Any four of the five independent quantities, V, M,, s, v, and w, 
suffice in each special cycle or outburst to determine, through the 
equations of Section II, the four adjustable constants, a, }, s, and F, 
of (1). Accordingly, five cases arise. 

1. Adopt from Table 2 M,, V, s, and v. Combination of (3) and 
(8) gives 


(v — ie = adg(a) . (11) 


‘0 Waldmeier, op. cit., p. 107, Fig. 1. 
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But by Table 4 the value of a@(a) so found determines a. Then, by 
(3), 
b= a (12) 





and we compute w by the second part of (6). 


TABLE 4 


AUXILIARY TABLE REPRESENTING EQUATION 9 








a (a) ag(a) (a+1) (a) 
£0: 0.3679 ©. 3679 0.7358 
£2. . 3402 0.4083 0.7485 
ey ae 3180 0.4452 0.7632 
PS; 2915 0.4955 0.7870 
BNO. 5657s vere 2707 0.5413 0.8120 
ee 2537 0. 5835 0.8372 
236... 2396 0.6230 0.8626 
20... 2240 0.6720 ©. 8960 
Cee ae 2082 0.7288 0.9370 
£50: .. 1953 0.7815 0.9768 
ALS... 1846 0.8308 1.0154 
SEIS. stale 1754 0.8774 1.0528 
ee 1675 0.9217 1.0892 
6.0 . 1606 0.9637 1.1243 
6.5 1544 I .0042 1.1586 
TiO; .1490 I .0430 I.1920 
nas 1440 1.0805 2.0115 
8.0 1390 1.1107 rE .25603 
8.5 1355 1.1518 I. 2873 
9.0 1318 1.1858 532570 
05... .1283 1.2189 I.3472 
10.0 1351 1.2511 1.3762 
BIOs 553 1194 53032 1.4326 
BIO os 5s ss .1144 1.3724 1.4868 
13.0.. 1099 1.4292 1.5301 
14.0.. 1060 1.4839 1.5899 
BS.0:.. 0.1024 1.5366 1.6390 














2. Adopt from Table 2 M,, V,s, and w. By (6) and (8), 


V ; 
(w— S)aF = (a+ 1)¢(a), (13) 


determining a by Table 4. Then, by (8), 


oe a+ 1 fen) 
A =F 





and we compute v by (3). 
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Ww 


. Adopt M,, V,v, and w. By (3), (6), and (8) 





and 





= ¢(a) ’ (16) 


(w — v) 


M, 


determining a by Table 4. We compute s by (3). 
4. Adopt M,, u, v, and w. Then (15) applies, but a is found from 





a= =. (17) 


and we compute V by (8). 

5. Adopt V, u,v, and w. Then (15) and (17) apply; so a and 0 are 
the same as in Case 4. We compute M, by (8). 

Values of the quantities w, v, s, V, and M,, so computed, will be 
designated by primes. In every case F is found by (5), which is 


equivalent to 
log F = (a + 1) log b — log '(a+ 1) + log M,. (18) 


As a step in these computations, Table 3 presents values of M,, 
v—u, V—S, w—s, and w—v derived from Table 2, together with 
values of V* and V/M,. The fact that w—v is less than zero for 
Cycles 1, 5, 7, 12, and 16 is definite evidence that the suggested 
family of curves described by equation (1) cannot accurately fit 
these cases, for a negative value of b in (15) is meaningless. But these 
cycles, together with 6, are just the group of cycles of lowest maxima. 
For these the observed centroid precedes the observed maximum; 
but for all curves described by equation (1) the centroid must fol- 
low the maximum. The conclusion is that the cycles of lowest maxima 
systematically deviate from the suggested form—their maxima come too 
late. This discrepancy is discussed further in Sections VI and VII. 
It corresponds in the detailed fits of Section VI to definitely system- 
atic but not very large differences between observed and computed 
spot numbers, especially for years near maxima of the group of low 
cycles. 
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Since the observed time v of maximum and the time w of centroid 
(Table 1) are both assumed in Cases 3, 4, and 5, these cases have 
meaning only when applied to the remaining cycles: viz., 2, 3, 4, 
6, 8, 9, 10, 11, 13, 14, 15. 

Table 5 gives for each of the sixteen cycles the values of a and b 
computed by the methods just described. The determination of the 
second decimal places in a and 0 is not precise. Only Cases 1 and 2 
apply to all the cycles. Cases 4 and 5 give identical values as re- 


TABLE 5 


COMPUTED PARAMETERS 








_ 

> 

wn 
w 





























CASE I CASE 2 CASES 4 AND 5 
CYCLE | eS 
NUMBER | | | 
a b a b log F a b a | b 
I | £0.43 WECSE DS 780) BSc) Oss te css ae 
2 | Fad E09) 1.2.98] £07 lpn Ose 1 2.97 | O.92 | 3.08 | 0.92 
3 | 30001 1.32) | S205. |-25600 | 4 25300 |. -1320 | 6-76. | 2.30: 10.76 
4 2.41 | 0.68 | 3.35 | 0.79 | +1.500 | 1.32 | 0.51 £262 10.81 
5 11.00 PS, GE I a i. 32 =O. 25% Mt Weak 
6 10,00 | 2.91 132.30: | 2.83 | —3.247 | 3074 | E02 | 0.44 | 2B .02 
7 PL. 22 | -0.55 | 7.85. | 1298 | —0-682 ; Pee 
8 i 3507 foots hh AnOO |e ie7 | etd | 20s 1058S) || -S 04-10-85 
9 | 8.63 a? Be 1.28 —0O.103 200° 150-77 | 3.49-| 6:71 
10 | 3.97 | o.g1 | 4.88 | 1.01 | +0.772 | 2.10 | 0.68 | 2.96 | 0.68 
II | 452010220 105.99 || 1738] 4-0 005 E 6x | 0:96. (2.92 |) 0.76 
12 | 7.07 | 1.31 | 4.62 | 1.07 | 0.044 xe | 
13 5.30.) 1.22) | 4705 |] 2-08.) 01622 | 8.00 | 1.37 | 6.58 | 1.37 
14 | 4.88 | 0.97 | 4.14,] 0.89 | +0.839 |10.48 | 1.41 | 7:090.| 3:41 
15 6.54] 1.54 | 7.04 | 1.60 | +0.550 | 4.36 | 1.27 | 5.38 | 1.27 
16 6.20 | £1.20 | 3.92 10.00 | =Fr.108 |... 
me 





gards a, and Cases 3, 4, and 5 give identical values of 6. Values of 
log F are included for Case 2 only; the others may be obtained from 
(138). 
VI. DETAILED CURVES 

The spread in Table 5, for the various cases, in the values of a 
and 6, respectively, for an individual cycle, provides an indication 
of how well the suggested form (1) fits each separate outburst. (But 
see the last paragraph of Sec. VIII). If the fit were exact in every 
detail, there would be no spread. If only random departures ex- 
isted, there would be no systematic spread. A serious systematic 
discrepancy has already been mentioned; it is discussed further 
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in the next section. Because of this discrepancy—-which proves that 
the suggested form (1) requires modification—detailed fits are pre- 
sented in Table 7 for only a few sample cycles and cases. Such fits 
are much better than a superficial inspection of the discrepancies of 
Table 5 would suggest. 

For example, the suggested form (1) represents the annual spot 
numbers of the long Cycle 4 very fairly—-notwithstanding small sys- 
tematic deviations; as Table 7 shows. It will be remembered that 


TABLE 6* 


COMPUTED CRITERIA 





| 
Goat | | 
Cycle Vv’ | Mo’ ieee sf | oe 











$ | =? 
Number 
I a (ey 5.98 
) 109.4 568 4.31 3.50 2.95 
3 118.7 | 822 3.78 3.72 1.58 
| 121.8 | 977 | 5.05 4.25 2.57 
5 8.21 oe! ae CSCI aan a Kt oe 
6 37.4 305 6.90 @.:72 3.07 
7 7.88 SiGe fi iSawsice paaree 
8 125.5 704 4.32 3.83 2.49 
9 102.8 | 889 5.82 6.67 2.90 
10 83.9 | 639 | 5.44 4.33 3:11 
11 i) 785 | 4.41 4.15 2.13 
I: ; 6.10 QoSe |  Roiveek oeeeas 
13 97.5 418 | 5-70 4.58 5-91 
14 77-9 307 6.00 4.95 7-44 
15 95.4 490 4.90 | 4.40 3.44 
10 6.00 4 o8 ; 
* The respective cases (Sec. V), left to right, are 4, 5, 1, 2, and 3. Cf. Table 3. 


this cycle of 1784-98 seemed so abnormal to Newcomb that he re- 
jected it in his statistical examination.’ Such rejection is altogether 
unwarranted from the present point of view, and invalidates New- 
comb’s conclusion that an underlying regularity of period governs the 
suns pot cycles. 

The late observed maximum of such a smallish outburst as Cycle 
12 is associated with very appreciable systematic deviations between 
observed and computed spot numbers throughout the cycle; but 
even so the representation (according to Case 2) computed in Table 
7 for Cycle 12 is in fair agreement with observation. 


ADS ., 13; 1; TOOl. 
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The computed spot numbers in Table 7 are obtained from the 
formula—another form of (2): 


log R’ = log F + a log (r — s) — 0.43429 b(r — 5s), (19) 


all logarithms being to Base to. In Cycle 4, for example, r takes on 
the successive values, 1785.5, 1786.5, ....; and the constant s 
(by Table 2) is 1784.53; log F, a, and 6 also are constant throughout 
the computation, and their values are taken from Table s. 
Examination of residuals O—C for the sunspot numbers, such as 
are presented in Table 7, indicates that a suitable mathematical 
form, replacing (1), which would eliminate the systematic discrepan- 
cies, would represent each outburst with an average random deviation 
regardless of sign of no more than half-a-dozen Wolf numbers. Thus, 
the hypothesis is confirmed that the general course of sunspot num- 
bers follows during each successive outburst an underlying and sur- 


prisingly definite trend. 
VII. SYSTEMATIC DEVIATIONS 

A comparison of the computed values presented in Table 6 with 
the corresponding values of the “observational criteria” given in 
Table 3 brings out the nature of the systematic failure of equation 
(1) to fit exactly the course of sunspot numbers. The overdetermina- 
tion of the four adjustable parameters permitted by the trustworthi- 
ness of the five observed quantities taken as criteria (Sec. III) for- 
tunately makes possible this test of the applicability of the form (1). 

The five cases represented in Table 6 must of necessity lead to 
the same conclusion with respect to the nature of the systematic 
discrepancies. It suffices to consider here in detail Case 2 only. The 
fit in Case 2, it will be recalled (Sec. V), is made by assuming M,, the 
area; V, the maximum smoothed monthly spot number; s, the time 
of beginning of the outburst (from eq. [10]); and w, the time of 
centroid of the outburst (from the first part of eq. [6], with a uniform 
plus end correction of 0.15 year); and computing the time of maxi- 
mum, v’. Table 8 compares these computed times, v’ (Table 6), with 
the observed times of maximum v (of Table 2) in relation to the 
observed heights V of maximum. Definitely, Table 8 shows v—v' is 
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large positively for the lower maxima and takes on negative values for 
the cycles of higher maxima. 


TABLE 7* 



































| 
| 
CASE 2 CASE 2 
r SR rs pee O-C —— ———_| O-C 
R’ R’ r R’ | 
= : 
Cycle 2 Cycle 12 | 
1700.5 0.3 1879.5 3.6 . 
67.5 22.1 TI5-7 80.5 $6.6 + 6.8 
68.5 79-1 =(O.3 81.5 56.8 — 2.5 
69.5 113.5 7-4 82.5 73.6 —tZ.9 
70.5 109.6 == S28 83.5 70.8 on 6 
71.5 4.3 9 84.5 56.4 + 7.1 
72.5 56.2 “BIG. 3 85.5 39.4 “F225 
72,8 34.0 + 8 86.5 25.0 + 0.4 
74.5 19.1 +11.5 87.5 14.8 — 1.7 
ae. & 10.1 Cycle 3 — 2.) 88.5 8.3 — 1.5 
76.5 5.1 8.3 + 6.4 89.5 4.4 | 
77-5 2.5 73-1 +16.9 
78.5 1.2 145-5 ae 
79.5 0.5 153.1 ~ 29.9 
80.5 0.2 Ei2.2 = 27:0 
81.5 O.1 65.6 +24 | 
82.5 32.5 6:8: | 
83.5 14.3 + 8.5 | 
84.5 Cycle 4 5.8 hh | | 
85.5 3423 2.3 + 8.6 Case 3, Cycle 4 
86.5 65.0 0.8 +17.1 1786.5 88.7 — 6.6 
87.5 116.9 0.3 +14.8 87.5 135.2 aad 
88.5 140.3 O.1 — 9.5 88.5 139.0 — 8.2 
89.5 135-5 waa) Ar 89.5 122.2 mee. 
90.5 114.0 — 24.3 90.5 98.4 — 8.5 
QL.5 87.1 26:5 QI.5 75.0 — 8.4 
92.5 62.1 pales. 92.5 55.1 + 4.9 
93-5 41.9 + 5.0 93-5 39.4 + 7-§ 
04.5 Lae +13.8 04.5 27.6 +13.4 
95.5 17.0 4:3 95.5 19.0 + 2.3 
96.5 10.4 + 5.6 96.5 12.9 +g 
OF. 8 6.2 | + 0.2 97.5 8.6 —i2 
98.5 3.6 | | + 0.5 98. 5 5.3 














* Where two computed cycles overlap, as in 1776.5, the two values of R’ were summed before com- 
puting O—C. (For the observed values see Table 1.) 
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The only marked exceptions to this correlation of v — v’ with V 
are Cycles 6 and 14. Cycle 6 in all statistical respects (Tables 3, 
8, 9) behaves like a cycle of high maximum, although it is the lowest 
of all. In view of observational uncertainties for the earlier cycles 
and of possible displacements of maxima due to the undoubted 
“random”’ deviations in spot numbers from their indicated under- 
lying trend, these exceptions are not necessarily significant. The 


TABLE 8 


THE SYSTEMATIC DISCREPANCY BETWEEN EQUATION (1) 
AND OBSERVATION 





| 





| 
V Cycle | Observed Computed O-C 
Number | v | v’ | : 
BS. 7 6 1816.4 1816.81 —0.41 
49.2 } 5 | 1805.2 | 1803.12 | +2.08 
64.2. | 14 | 1900.4 | 1906.02 +0. 38 
if eer (e 7 1829.9 | 1828.55 +-1.35 
74.6 12 1883.9 | 1882.84 +1.06 
78.1 160 1928.4 1927.33 | +1.07 
80.5. I 1761.5 | 1760.58 +o.92 
87.9 13 1894.1 1893.88 | +o0.22 
97.9. 1860.1 | 1860.58 —0.48 
BOS: 4. : 15 | I9gI7.6 | 1917.70 0.10 
cEs.S. . 2 1760.7 1769.85 | -O.15 
131.6.. 9 1848.1 1848.75 | —0.65 
140.5. U1 1870.6 1871.18 | —o.58 
cer. 2. 2. | 4 1788.1 | 1788.78 -0.68 
£40:9..... 8 | 1837.2 | [53757 | =0:.37 
RSG, 52%. | 3 , 1778.4 1779.09 | —0o.69 
} | 





observed spot numbers for Cycle 14 obviously are a bit ragged 
it seems that “random”’ deviations in the actual number of spots 
from their underlying trend are pronounced for that cycle. The 
discrepancies for Cycle 6 would be avoided if the observed maximum 
came later. On a flat-topped curve with random deviations such an 
error would not be surprising. 

A rough empirical representation of Table 8 is 


s 108 — J 
Ee gr = “<n Seaae ay (20) 


40 
Search for an improved mathematical form to be substituted for 
(1) should be guided, therefore, by (20). That is to say, a new family 
of curves should be looked for such that, when the height of maxi- 
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mum, V, area, M,, time of start of outburst, s, and time of centroid 
of outburst, w, are fitted exactly, the new computed time of maxi- 
mum v” will follow the corresponding Case 2 computed time v’ 
by roughly the quantity (V — 108)/40. This linear correlation with 
V need not be anything like exact. 

The existence of the systematic deviations shown in Table 8 is 
definitely confirmatory of the hypothesis that a suitably modified 
mathematical form rather similar to (1) will represent with fair ac- 
curacy the general course of sunspot numbers. The various numerical 
tables already presented lead to the very tentative suggestion that 
the true form replacing (1) likewise needs to have four adjustable 
parameters. The wide variations in a, b, and log F (Table 5) per- 
mitted in almost equally good representations of individual cycles 
suggest that (1) uses four parameters to make adjustments for which 
only three would be sufficient—resulting in some indetermination. 
On the other hand, the inability of the form (1) to adjust itself 
generally to the observed relation between maximum and centroid 
(Sec. V, Cases 3, 4, 5; Cycles 1, 5, 7, 12, 16) shows that a fourth 
parameter is required in the true form to meet this condition. 


VII. THE SECOND MOMENTS 

Were it not for the large, uncertainly determined end corrections 
which are required (Sec. IV, because of overlap, for the values of M, 
presented in Table 2, equation (7) would give a sixth control for the 
determination of a and 6, and so would give an excellent further test 
of the applicability of form (1). In the actual circumstances this 
test is a weak one. Table 9 compares the computed values of 
(a + 1)/b? (from Table 7 for Case 2 only) with the observed values 
of o? from Table 2. If form (1) applied without systematic discrep- 
ancy, the observed o’ of the cycle could not be systematically greater 
than the computed (a + 1)/b?, for the expected end corrections to 
o? are positive. But for Cycles 1, 3, 4, 6, 8, 9, 11, 13, 15, o? is greater 
than (a + 1)/b*. Of these only Cycle 1 is among the group of cycles 
of low maximum for which the centroid precedes the maximum. 
That group comprises (as has been said in Sec. V) Cycles 1, 5, 7, 12, 
16. The condition that the observed o? be less than the computed 
(a + 1)/b? is indeed met (when the fit is by Case 2) only by Cycles 
2, 5, 7, 10, 12, 14, 16. (If Case 1 instead of Case 2 had been exam- 
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ined, the only change is the reversal of size of the O—C for Cycle 4.) 

Accordingly, Table 9 gives further, although slighter, evidence, 
supplementing Table 8, that the suggested form (1) is systematically 
somewhat in error. 

It is worth mentioning that the O—C in Table g for Cycle 2 is 
about what would be expected from end correction. Also, in Table 8 
Cycle 2 shows only a small discrepancy. Therefore, equation (1) 
applied to Cycle 2 fits well all six criteria. The height of maximum 

TABLE 9 


THE SECOND MOMENTS 

















Cycle (a+1)/6? . 5 

: ‘ ao Oo—( 

Number (Case 2) 
I.. 5.53 5.68 +0.57 
2 4.15 3.81 =0..34 
3 2.71 3-17 +0. 46 
4 6.97 7.19 +0.22 
5 5-75 4.68 1.07 
6 3.97 4.95 +o0.98 
7 eee 4.19 —1.02 
8 3-64 3-97 +0. 33 
Q. 4.96 6.50 +1.54 
IO. 5.76 5.61 —=—O.55 
II. 3-53 4.25 +0.72 
12 4.91 4.44 —0.47 
13 5.10 S29 +o.14 
14 6.49 4.67 O52 
15 | 3-14 3-71 +0.57 
16. | 5-34 2.38 —1.96 





of this cycle is intermediate between the two extreme groups, which 
indicates that equation (1) 1s exact for cycles of intermediate height of 
maximum, and the run of values of v-v' in Table 8 confirms this. 
IX. REPRESENTATION OF THE COURSE OF SUNSPOTS 
BY A TWO-PARAMETER FAMILY 

Referring to Table 5, Case 2, 6 and log F can be expressed with 
fair statistical approximation as functions of the corresponding a. 
The empirical approximations 


b 


1.6ologa+ 0.03, (21) 


log F = —0.537a + 3.63, (22) 


represent most of the sixteen cycles not badly (Table ro). 
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A series of curves R as a function of 6 (eq. [2]) for different rele- 
vant values of a, and having 6 and F related to a by (21) and (22), 
could be computed and plotted on tracing paper. If all the sunspot 
numbers of Table 1 are plotted, then by sliding in time sample curves 
of this series the best fit with respect to values of a and s can be 
readily found for each cycle. These two-parameter representations 
so found of the whole course of the sunspot numbers would be, in general, 
fair, although in certain cycles not so good as the four-parameter 


TABLE 10 


CORRELATION OF a, b, AND F (CASE 2) 


| | 








Cycle 





| 
Number | $ | log F | log ¢ | . 
: | eee 
4 3:35 +1.500 0.53 °.79 
2 3-75 +1.650 0.57 1.07 
16 3.92 +1.198 °.59 °.96 
14 4.14 +0 .839 0.62 o.89 
12 4.62 +0.944 0.66 1.07 
8 4.86 +1.447 0.69 : a7 
£0... . 4.88 +0.772 | 0.60 1.01 
13 4.95 +0.823 0.69 1.08 
Pf. 5.92 +1.095 0.76 1.38 
3 5.05 +1.390 0.77 1.60 
5. 0.21 =), 234 ©.79 E.32 
15 7.04 +0.550 0.85 1.60 
0... 7, 04 —0.103 0.85 1.28 
a 7. $5 —o.682 0.88 1.28 
- 7.80 —0.763 0.90 1.32 
6 12.30 — st Tay 1.09 1.83 











fits of Table 5. A careful study of the overlaps so found doubtless 
would result in a slight modification of equation (10) and of the 
correction to M,/M,. Of course, all these representations, which 
follow form (1), are subject to systematic deviations described in 
sections 7 and 8. 

Waldmeier’s suggestions” related implicitly to a two-parameter 
representation of each outburst. A definitive decision as to the 
validity of such a representation and of the physical usefulness of 
formulae such as (21) and (22) must be postponed until an improved 
substitute for (1) has been found, and such a substitute is not pre- 
sented here. In any event, with so few cycles available, the statisti- 


2 Waldmeier, op. cit., Figs. 2 and 5. 
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cal situation may remain too indefinite to permit the decision for 


decades to come. 


X. THE PROBLEM OF PREDICTING THE TIME OF START OF AN OUTBURST 

A similar study relates to the possible correlation between the 
time of the start of an outburst and characteristics of that outburst 
or of the preceding one. Jf such a correlation exists, then s ceases to 
be an independent parameter; and a one-parameter representation is 
possible. Evidently, the new outburst never begins until the old one 
has faded rather low. For example, one finds evidence to support 
the very tentative rule that the new outburst begins about a year after 
the old one has dropped to one-tenth of its maximum. This problem 
could be readily given further study, employing data in the fore- 
going tables as a basis for computation. 


XI. THE AREA AHEAD OF MAXIMUM 


Waldmeier has stated the proposition that in every outburst the 
sum of all the monthly spot numbers between the preceding mini- 
mum monthly number and the maximum monthly number of the 
outburst is constant. Observed data which he presents roughly 
bear out this suggestion, but in this part of his study he has omitted 
Cycles 1-6, inclusive—and Cycles 5 and 6, with the smallest maxima 
of all, if included, would definitely disprove his suggestion. There 
seems to be no adequate statistical reason for omitting the earlier 
cycles—as examination of the tables of the present paper shows. Of 
course there is no denying that the older data are uncertain, but 
their inclusion seems necessary. 

From the standpoint of the present paper, a brief treatment of the 
area ahead of maximum is of some interest. The fraction a of the 
whole area M, of the outburst which is included ahead of the maxi- 
mum may be found from tables for the incomplete gamma function. 
Equation (3) shows that the maximum is always at 6 = a/b. The 
fractional area included between @ = O and @ = a/b is readily 
proved to be a function of a only, independent of 6; and, referring 
to Pearson’s tables? of I (u, p), ais I(a/Va + 1), and is independent 
of 6. Values are presented in Table 11; a is nowhere as large as 0.500. 


13 [bid., p. 120. 
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For any given outburst, the computed area ahead of maximum 
would beaM,. The observed values given by Waldmeier are subject 
to much uncertainty, in that a slight shift of the times of least and 
greatest monthly numbers obviously changes them relatively more 
than the intervals vw are changed. 

XII. CORRELATIONS AMONG THE OBSERVATIONAL DATA 

a) The very significant correlation between height of maximum, 
V, and the time of rise, v-u (Table 3), of a cycle, has been redis- 
covered independently by several investigators, but even yet is 

TABLE 11 


FRACTIONAL AREA, a, AHEAD OF MAXIMUM 











a a a a 

1.0 0.204 2.9 ; 0.351 
i.e: 281 3:2 357 
B:3 289 3.5 303 
1.4 294 3.9 309 
E.5 300 4-3 370 
¢. 308 4:7 381 
EO. . 310 5.0 | 354 
a0 . 320 0.0 | 3905 
2.3 335 7-0 399 
2:5 341 49.0 0.462 
2.7 0.345 

















frequently ignored." This relation is well satisfied by the methods 
of fit based on equation (1). 

b) The relation between the time of rise and the time of fall often 
has been inadequately stated. The ratio of these depends upon the 
height of maximum. Similarly, an excellent correlation—to which 
explicit attention appears not to have been directed prior to the 
present study—holds between V and w-v: that is, between height of 
maximum and interval between maximum and centroid. This same 
phenomenon, however, previously has been implicitly recognized.’ 

c) A third, rougher, correlation probably exists, when annual 
spot numbers are studied, between the height of a minimum and 

14 A synthetic treatment of the spot cycle by the authors (Pub. A.A.S., 9, 21, 1936) 
adopted the outburst hypothesis but did not take this correlation into account. 


1s Waldmeier, op. cit. 
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the height of the following maximum. Minima which are higher 
than average tend to be followed by especially high maxima. To the 
extent to which equation (1) applies, this tendency is to be expected 
because a rapid initial increase in spot numbers at the beginning of 
an outburst is likely to result in a relatively high spot number for 
the year which includes that minimum and because such a rapid 
initial increase is associated with a high maximum. 

d) Section X tentatively suggests an additional correlation re- 
lating to the time of start of an outburst. 

e) Table 10 exposes an additional and remarkable correlation, 
which scarcely can be the result of chance: the even-numbered cycles 
all have lower values of a (Case 2 solutions) than any odd-numbered 
cycle (the only exception is Cycle 6). This result confirms the system- 
atic difference in alternate cycles shown by Ludendorff"® and gives 
a basis for the twenty-two-year “‘period’’ which sometimes has been 
suggested. The well-known reversal of polarization of the spots in 


alternate cycles is of great interest in this connection. 


XIII. FAILURE OF THE SUPERPOSITION HYPOTHESIS 

A Fourier analysis—and a fortiori periodogram analyses which 
do not restrict themselves to harmonic components—can account 
for any course of sunspot numbers whatsoever, once the sunspots 
have passed into history. The improbability, obviously, is tremendous 
that such components would ever combine to give a succession of super- 
ficially dissimilar outbursts each of which nevertheless exhibits common 
family characteristics to the degree that the present study reveals. 

In addition to this argument, the failure of superposition treat- 
ments to predict successfully seems conclusive. The outburst 
hypothesis evades the test of long-time prediction because, until 
many more sunspot cycles have elapsed, data will be lacking for 
discussion of what weak statistical relations there may be between 
the parameters of successive outbursts. Short-range prediction 
ought to be fairly successful for the remaining course of an individual 
cycle which has reached or passed its maximum. The methods indi- 
cated in Section [IX would apply to this problem—for example, to 
prediction of the remaining years of the present cycle (No. 17). 


6 Zs. f. Ap., 2, 370, 1931. 
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Further study along the lines which we have indicated may dis- 
cover another direct analytic formula which retains the many ad- 
vantages and avoids the defect of equation (1). If not, it should not 
be very difficult, retaining equation (1) as a first approximation, to 
modify the treatment empirically in such a way as to eliminate the 
discrepancies described by Tables 8 and g. The definitive study 
doubtless ought to be made directly with the monthly, instead of 
with the annual, sunspot numbers. When either of these alterna- 
tives has been accomplished, it will be desirable to publish full tables 
for facilitating short-range prediction of the smoothed run of sun- 
spot numbers. 

PRINCETON UNIVERSITY OBSERVATORY 
June 25, 1938 








RADIAL VELOCITY-CURVE OF THE RR LYRAE 
VARIABLE W CANUM VENATICORUM* 


ALFRED H. JOY 


ABSTRACT 

Twenty-four spectrograms of the short-period variable W Canum Venaticorum 
were obtained, mostly by the use of a 6-inch camera and a one-prism spectrograph 
giving a dispersion of 120 A/mm. The velocity-curve closely resembles the light-curve 
by Detre, but a lag of 6 per cent of the period is noted. The normal velocity is 
+26.0km/sec, with a total range of 70 km/sec. 

The spectroscopic absolute magnitude of —o.4 and a distance of 1400 parsecs were 
determined. 

Of the short-period variables of the RR Lyrae type the velocity 
variations have thus far been published only for the type-star itself. 
The lack of observations may be attributed largely to the faintness 
of the stars, all of which fall well below the ninth magnitude in 
brightness at minimum light. Another difficulty is that the change 
in velocity is so rapid at certain times that exposures of an hour or 
more would involve doubtful corrections depending on both change 
in velocity and change in light. 

Observations of W Canum Venaticorum were undertaken in order 
to test the applicability of low dispersion to this problem. A 6-inch 
camera giving a dispersion of 120 A/mm at Hy was used with the 
one-prism spectrographs at’ the 60- and too-inch telescopes. With 
this camera the exposure times varied from 8 to 50 minutes. Four 
additional plates, C2307, C2624, C7162, and 721356, with a disper- 
sion of 75 A/mm at Hy, are also available. The exposure times for 
the latter were 120, 165, 165, and 65 minutes, respectively. 

In Table 1 are listed the data concerning 24 plates of W Canum 
Venaticorum. The elements, 


Max. = J.D. 2421402.4206 + 0.55175981E , 


by Detre’ were used. The phases are expressed in decimal fractions 
of the period, and the estimated relative weights are given in the 

* Contributions from the Mount Wilson Observatory, Carnegie Institution of Washing- 
ton, No. 595. 


© A.N.; 294; 21, 1034. 
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last column. The velocities are plotted in Figure 1, and a rough 
curve is drawn to picture the course of the velocity variations. 
Detre’s visual light-curve (range 9.9-10.7 mag.) is sketched above for 
comparison. A remarkable resemblance between the light- and ve- 
locity-curves is evident, but, as was found for variables of the 
6 Cephei type? and for RR Lyrae,’ there is a distinct lag of the veloc- 


TABLE 1 


OBSERVATIONS OF W CANUM VENATICORUM 











‘ _ } Vel. 

Plate Date G.M.T. J.D. 242+ Phase (Km/Sec) Wt 

GC. 2367... 1923 June 23 3504.757 0.353 +24.9 Es 
2624 1924 Jan. 16 3801 .o81 .Q21 +24.9 1.0 
721286 1938 Feb. 5 8935 .899 547 +42.2 1.0 
21287 5 8935.928 600 +52.8 0.3 
21288 5 8935 .958 -654 +43.3 0.6 
21289 5 8935.988 .707 7.6 1.0 
21290 6 8936.021 . 768 +41.4 1.0 
21308 Apr. 0 8998 .922 768 +34.2 1.0 
21309 9 8998 .958 833 +63 .6 1.0 
21310. 9 8998 .9QI1 804 +40.6 1.0 
21311 10 8999 .017 .940 +11.0 0.6 
21312 IO 8909 .033 .970 + 9.1 O.3 
2TSE2.. 10 8999 .640 .070 —16.3 0.6 
21314 10 8999 .655 .097 —1I5.9 o.3 
21315 10 8999 .672 527 — 32.0 0.3 
21316 10 8999. 6902 164 — 626 0.6 
21318 10 8999.751_ | .271 +26.1 1.0 
21319 10 89909 .774 JSES +18.9 1.0 
© 7162 May 6 go25.908 677 +48.0 ¥.5 
7163 6 9025.990 826 +38.0 0.6 

¥ 21344 fe) 9028.79! .9O3 +63.1 0.6 
21345 9 9028 .822 959 + 9.0 | 
21346 9 9028 .831 .975 7.9 1.0 
21356 July 11 QOQI.755 0.017 + 8.6 rs 

















ity-curve as compared with the light-curve. The lag amounts to 
about 6 per cent of the period at both maximum and minimum 
phases and is closely the same as that found for RR Lyrae by Kiess 
and by Sanford. The normal velocity line is drawn so as to give equal 
areas of the velocity-curve above and below the line. 

The normal velocity is found to be +26.0 km/sec, and the range 

2 Mt. W. Contr., No. 578; Ap. J., 86, 435; 1937. 

3 Kiess, Lick Obs. Bull., 7, 148, 1913; Sanford, Mt. W. Contr., No. 351; Ap. J., 67, 


325, 1928. 
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is 70 km/sec (from —22 to +48 km/sec). In Mt. W. Contr. No. 578 
the mean range for Cepheids was shown to be about 42 km/sec. If 
we may judge from RR Lyrae and W Canum Venaticorum, whose 
ranges are 61 and 70 km/sec, respectively, the velocity range of the 
short-period variables is at least 50 per cent greater than the mean of 
the 6 Cepheids and double that found for stars with periods between 
one and two days. 
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Fic. 1.—Light- and velocity-curves of W Canum Venaticorum. The visual light- 
curve (above) was observed by Det're with a wedge photometer. In the velocity-curve 
(below) the individual observations are plotted. 


The spectral type of W Canum Venaticorum is found to be A6 at 
maximum and F6 at minimum. The spectroscopic absolute magni- 
tude determined from the usual A- and F-curves is —o.4, which cor- 
responds to a distance of 1400 parsecs. The galactic co-ordinates are 
1 = 36°,b = +70. 

The results of this investigation indicate that much interesting 
and valuable information in this unexplored field may be obtained 
with a moderate expenditure of time and effort by the use of low- 
dispersion spectrographs. 

CARNEGIE INSTITUTION OF WASHINGTON 
MowntT WILSON OBSERVATORY 
July 1938 














PHOTOGRAPHIC LIGHT-CURVES OF THE TWO 
SUPERNOVAE IN IC 4182 AND NGC 1003* 


W. BAADE AND F. ZWICKY 


ABSTRACT 

Photographic light-curves of the two supernovae in IC 4182 and NGC 1003 have 
been derived from plates obtained at the Palomar and the Mount Wilson observatories. 
Since the distance moduli of the two nebulae are known, the absolute photographic 
magnitudes of the supernovae can be computed. At maximum their values were 
M=—16.6 and —14.0, respectively. An integration of the light-curve of the super- 
nova in IC 4182 shows that during the outburst this star emitted radiation in the 
photographic region equivalent to 4-107 years of solar radiation within the same 
spectral range. 

When we discussed some years ago the information then avail- 
able about supernovae," we realized that further advances in this 
field could be achieved only by a systematic search for these ob- 
jects. Although fourteen supernovae had been discovered up to 
that time, the resulting data were so fragmentary that little more 
could be derived from them than the fact that supernovae as a class 
do exist. Since supernovae are very rare objects, even if large num- 
bers of nebulae are kept under control, it was clear from the begin- 
ning that the time and effort necessary to carry through a system- 
atic search program could be justified only if a proper search in- 
strument were provided. As to the choice of instrument, there was 
no doubt that a Schmidt telescope of approximately 20 inches aper- 
ture would best meet the requirements of the proposed program. 
Sooner than expected, such an instrument became available; for 
early in 1935 the Observatory Council of the California Institute 
of Technology, whose members took an active interest in our plans, 
authorized the construction of an 18-inch Schmidt telescope? for the 
Palomar Observatory. The instrument, which has been constructed 


* Contributions from Mount Wilson Observatory, Carnegie Institution of Washington, 
No. 601. 

1 Mt. Wilson Comm., No. 114; Proc. Nat. Acad. Sci., 20, 254, 1934. 

2 The telescope, of the focal ratio F/2, has a correcting plate of 18 inches and a 
mirror of 26 inches aperture. The field covered by the photographic film has a di- 
ameter of 9°5. 
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under the direction of Dr. J. A. Anderson, was ready in the fall of 
1936. Since that time it has been used in an extensive search for 
supernovae by Zwicky. So far, three have been discovered by ob- 
servations made on Palomar Mountain. Of these, the supernovae in 
IC 4182 and NGC 1003, discovered in the fall of 1937, have yielded 
valuable information because they were relatively bright objects 
which could be followed photometrically and spectroscopically over 
a long interval of time. Both are still under observation and can be 
followed probably through another season before they pass out of 
the reach of the roo-inch telescope. To avoid undue delay in the 
publication of the photometric results thus far obtained, we present 
in this paper the photographic light-curves* of the two supernovae 
from the time of their discovery up to the time of our latest observa- 
tions. 
THE SUPERNOVA IN IC 4182 

The spiral nebula IC 4182.—This nebula (R.A. 13"2™g9, Dec. +37° 
56’, 1938; gal. long. 67°8, lat. +78°9) is a late Sc spiral of 4/1 diam- 
eter, oriented nearly perpendicular to the line of sight. It has a very 
low surface brightness, which necessitates long exposures to reveal 
its structure.4 The integrated photographic brightness of the nebula 
is correspondingly low, being only of magnitude 13.5.4 Nevertheless, 
it is easily resolved on plates obtained at the Mount Wilson re- 
flectors from which the value m, = 19.2 has been derived for the 
mean apparent photographic magnitude of the five brightest stars 
in the system. Since for the average late-type system the correspond- 
ing value of the absolute magnitude is M, = —6.3,5 we obtain the 
modulus m — M = 25.5, whence, for the distance of IC 4182, D = 
1.02:10° parsecs. The resulting values for the linear diameter (1200 
parsecs) and the luminosity (M = —12.0) clearly indicate that IC 
4182 is a dwarf system. For such systems, however, the absolute 
photographic magnitude of the brightest stars is decidedly fainter 
than —6.3, and rather of the order M, = —5.6. Substituting, 

3 Light-curves based on observations in the red region of the spectrum will be 
published after a sufficiently extended scale of red magnitudes has been established. 

4 For an illustration of IC 4182 see Mt. W. Contr., No. 600; Ap. J., 88, 285, 1938. 


5’ E. Hubble, Mt. W. Contr., No. 548; Ap. J., 84, 158, 1936. 
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therefore, this smaller value as a second approximation, we obtain 
the following final data for IC 4182: 


Modulus...... 24.8 

Distance. 912,000 parsecs 

Diameter... 1100 parsecs 
Integrated absolute magnitude ... .—11.3 (photographic) 


It was a very fortunate circumstance that these data were known 
to us when we compiled the observing list for our search program.° 
Since any supernovae appearing in IC 4182 should have, in the 
mean, the maximum apparent brightness m,, = 10.5, we included 
the nebula in our observing list, although the chance for the dis- 
covery of a supernova did not seem to be very high, on account of 
the small stellar content of the system. Within less than a year, 
however, one of the brightest supernovae on record appeared in 
IC 4182. It offered an unusual opportunity for extensive photo- 
metric and spectroscopic observations. Because of its brightness and 
high declination, the nova could be followed through the conjunc- 
tion with the sun, so that an uninterrupted record of its history has 
been obtained. 

Photographic magnitudes of the comparison stars.—On account of 
the large observable range of the supernova, an extended sequence of 
comparisons had to be established, involving instruments of widely 
different apertures. To cover the interval from magnitude 8 to 11, 
polar comparisons were made with the 5-inch Ross lens (intrafocal) 
and the ro-inch Cooke refractor of the Mount Wilson Observatory. 
Magnitudes between ro and 13 were derived from polar comparisons 
made with a schraffierkassette at the 60-inch reflector, which had 
been diaphragmed down to 4o inches. For stars fainter than 13, 
photographic magnitudes have been determined from focal exposures 
at the 60-inch reflector, using either the Polar Sequence or S.A. 57 
as standard regions. The magnitudes of the comparison stars are 
therefore on the international system. Because of poor seeing con- 
ditions during the last spring season, a rather large number of com- 
parisons had to be made to reduce sufficiently the influence of the 
systematic errors with which the individual plates may have been 

6 The nature of IC 4182 was discovered by Baade in 1936 during a search for near-by 


extragalactic systems of low luminosity. 
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affected. The final photographic magnitudes of the comparison stars, 
together with the residuals of the individual determinations, are 
given in Table 1. For identification of the stars see Plates XVII and 
XVIII. The sequence extends at present to photographic magnitude 
17.0, which should be sufficient for observations of the supernova 
during this season. An extension of the sequence down to magnitude 
21 is under way and will be published later. 


TABLE 1 
SEQUENCE FOR SUPERNOVA IN IC 4182 
(Unit for Residuals, 0.01 Mag.) 












































| 
| RESIDUALS 
. MEAN r 
STAR Mg Schraffier- BD 
5-Inch 10-Inch kassette, 60-Inch 
60-Inch 
eae 8.37 —4i-+ 2|+ 2 | ; | ; | | fe + 38°2306 
a......| 8.66] +3/+ rl— 5 | | 38 2403 
Bic ais 8.93} +4] of- 3 37 2371 
Oiss<es Q.1r] +1/— 31+ 3 37 2360 
d 9.74, —1;—11/+ 8/+ 5 ° 37 2307 
e. 9-79] +2]/— 2]+ 8+ 4]/+ 2/— 5 o]— 9 38 2308 
f 10.03] +4]/— 5|/+10/+ 3 ol— 5 oj— 7] 38 2402 
g. 10.79]... sg) (OS cece". Numa hy ete aS fae ge 38 2307 
h. 10.85] —3]—20 -[— 8+214/+ 11+ 4/41 38 2405 
Beis 10.78] +7/— 5 mae ice | Ree) oe yeh 37 23608 
5 11.46 —18/+ 5|/-— 1 3 5 
i II.99 ‘ — 21-1 4 
m 12.07 — 4/— I+ 5 
Be. . 12.57 ° oj+ 1 | 
o 13.25 [—rr]+ r/— 6/— 1/+ 5 + 4)+12 3 
| 

p 13.72 meee 3}—19}+ 3}/+ 8/+13)/+19|/ +20 2} 
q 14.02 | + 21+ 2]/— 2]/+ 2/— 2 o} 
r 14.21 | TIai— ti— 41+ Ji~> 7 5 
s 14.73 ATIOl— Fit I+ Sim Fim 5S 
t 15.22 ..-|—12]+ 4/+ 8]/-— 1 I 
u 15.78 | ti O1—I51-+ 4 
Vv 16.07 | on rr2|—11 I 
Ww 16.34 ae II 
oe 16.70 | —10/+ 9 
Vissencp Sa0Od | | | | — sit 5 









































Photographic light-curve of the supernova.—The observed photo- 
graphic magnitudes of the supernova from August, 1937, until June, 
1938, which have been obtained with the Palomar and the Mount 
Wilson instruments, are collected in Table 2. 

The resulting light-curve is reproduced in Figure 1, in which the 
individual observations are marked by dots. It shows that the super- 
nova was discovered when it was already past maximum. A pre- 


























PLATE XVII 


COMPARISON STARS FOR SUPERNOVA IN IC 4182 
Scale, 1 mm g2"5. Enlarged 2.5 times from a film taken with 18-inch Schmidt tele 
scope on Mount Palomar. 




















PLATE XVIII 


ee 
t— 
o— 
s— 
— 
u-— 
—vV 
—x 


—Z 


COMPARISON STARS FOR SUPERNOVA IN IC 4182 


Scale, 1 mm 
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TABLE 2 


PHOTOGRAPHIC MAGNITUDES OF SUPERNOVA IN IC 4182 



































see SSS 
| | 

J.D. G.M.T | _™pg | Instr II _# ssc wee: - Mog ‘ — 

2428774 .7 8.58 S* || 2428808 .6 11.36 | § 
275.7 8.63 S I 808 . 6 } 11.35 | S 
776.7 8.62 S | 810.6 [ Rese | S 
98 a 8.64 S | 813.6 11.46 | S 
778.7 8.74 S | 815.6 11.62 S 
779.6 8.83 S | 818.6 Fa: 79 S 
770.7 | 8.98 S | 820.6 11.86 S 
780.6 9.25 S {| 831.0 | IIl.94 S 
780.7 |} 9.14 5 | 833.0 | 12.03 S 
781.7 | 9.13 5 i 837.0 le 3202 = 
781.7 | 9.09 S || 842.0 | 12.05 S 
781.7 | 9.24 S | 862.9 | 32 38 S 
782.7 9.50 2 H} 865.0 | 12 40 10 
783.7 | 9.54 S | 865.0 12.46 10 
783.7 | 9.55 S | 865.0 12.42 IO 
784.6 | 9.58 > 865.0 12.43 fe) 
784.6 | 9.55 S 866.0 | 12.49 60 
785.6 f  O5.78 S | 866.0 | 12.43 60 
785.6 | 9.977 S | 871.9 | 12.47 S 
786.6 9.81 3 1 873.9 | 12.47 S 
786.6 | 9.89 S | 880.9 | 12.64 S 
787 6 9.93 S | 805.9 | 34,84 60 
787.6 | 9.93 S | 596.0 | 13.15 60 
788.6 | 10.08 S 896.0 53.35 60 
788.6 | 10 30 S | 896.0 | 13.41 60 
788.7 | 9.092 5 | 897.0 ae ee S 
789.6 | 10.11 S | 899.0 [ F233 2 
789.6 | 10.45 > | 904.9 | -¥3.373 S 
789.7 | 10.2 5 | 906.9 | 13.05 S 
790.6 10.22 Ss i 907.9 | 3.37 S 
790.6 | 10.37 S | 908 .9 | mz. 30 S 
790.7 | 10 32 5 | Q10.0 | 13.2 S 
791.6 | 10.37 S | Q10.9 | -¥3:30 S 
7901.6 | 10.54 S i 921.9 13.48 S 
792.6 | 10.56 S | 947.8 | 3¥.65 S 
792.6 | 10.41 S | 949.8 | 13-74 s 
793.6 10.52 S i 950.8 13.72 S 
793.6 | 10.56 S | 951.8 | 13.71 S 
794.6 | 10.79 S I] 952.8 14.13 S 
794.6 | 10.68 S | 052.8 14.13 S 
799.6 | 3.0 S 954.8 | 13.92 100 
799.6 | 11.3% SS | 955.8 | 13.94 100 
800.6 | trac S | 963.9 | 14.10 S 
800.6 ! 76 S ] 981.8 | 14.32 S 
801.6 | oa 30 S } 985.7 14.35 S 
801.6 11.16 S | 988.7 14.34 60 
802.6 | BP 27 S } 988.7 14.44 60 
802.6 | 411.25 S | 989.7 14.21 S 
803.6 | hay S 990.8 14.57 60 
803.6 II.26 S | 992.8 14.66 S 
804.6 | Ir.33 S | 2429010.8 | 14.68 S 
804.6 11.26 S O12.7 14.90 S 
806.6 | 11.29 S 024.6 | 15.07 100 
806.6 TT. 31 S 040.8 15.28 S 
807.6 II.42 S | 048.8 15.57 60 
807.6 II.39 S 049.8 <3) Sg SF 60 


| 





* 18-inch Schmidt. Figures in this column indicate the apertures of the various telescopes used. 
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maximum observation by E. Leutenegger’ makes it possible, how- 
ever, to establish with sufficient accuracy the date and the bright- 
ness of the maximum. Leutenegger, who had photographed comet 
Finsler 1937 on August 16.9, U.T., found, upon a later examination 
of this plate, that it included IC 4182 and that the supernova was 
already present as a star of the ninth magnitude. At our request he 
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Fic. 1.—Photographic light-curve of supernova in IC 4182 
kindly redetermined the magnitude of the supernova, using our 
system of magnitudes, with the following result :° 


JD 2428762.4 G.M.T. mpg = 8.85. 


The corresponding point is marked by a cross in Figure 1. There is 
hardly any doubt that the observation belongs to the ascending 
branch of the light-curve. We thus obtain for the date of the maxi- 
mum August 22, 1937, and for the apparent maximum brightness, 
Mpg = 8.2. The light-curve itself resembles in many ways that of 
an ordinary nova. The decline is rather steep during the first 30 

7A.N., 263, 446, 1937. 

’ We are much obliged to Dr. Leutenegger for his co-operation. An enlarged print 
of his plate which he kindly sent us leaves no doubt as to the correct identification of 


the supernova. 
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days following the maximum, with a gradient of 0.10 mag. per day. 
After that a gradient of 0.018 mag. per day represents the further 
decrease in brightness to the end of June, 1938. 

Absolute photographic brightness of the supernova at maximum. 
From mp, = 8.2 for the apparent brightness of the supernova at 
maximum and the known modulus of IC 4182, m — M = 24.8, we 
obtain for the absolute magnitude of the supernova at maximum 


Murex == 1050). 


This is the highest value so far observed among supernovae, and 
2.3 mag. brighter than the average supernova at maximum. Since 
the integrated absolute magnitude of IC 4182 is only —11.3, the 
supernova at its maximum surpassed in luminosity the nebula in 
which it originated by a factor of 100. The data already given fully 
explain the unusually large value Mmax —M!neb = +5.3. It is due 
simply to the fact that a supernova of very high luminosity appeared 
in an extragalactic system of very low luminosity. Moreover, the 
difference is numerically well within the range allowed by the known 
dispersions in the luminosities of spiral nebulae and supernovae. 
Finally, we compare the radiation of the supernova with that of 
the sun. Since the photographic absolute magnitude of the sun is 
M = +5.3, the supernova at its maximum surpassed the luminosity 


of the sun by a factor of 6-10’. 


To obtain the total energy in the 
photographic region of the spectrum which the supernova radiated 
away during the outburst, we have integrated the light-curve over 
the first 225 days after the maximum. The total corresponds to 
4:107 years of solar radiation within the same range of the spectrum. 
These two figures represent what may be said at the present time 
about the energy involved in a supernova outburst without invok- 
ing assumptions. Obviously, they are a mere fraction of the total 
energy emitted between \ = o and \ = ©, which cannot even be 
guessed at because we do not know the intensity distribution in the 
spectrum. In the observed region from \ 7000 to X 3300 the distribu- 
tion is so peculiar that the assumption of black-body radiation for 
the supernova can be ruled out at once. In the photographic region 
a strong band at \ 4600, approximately 200 A wide, is the dominant 
feature. The fact that about half of the radiation emitted in the 
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photographic region between X 5000 and A 3300 originates in this 
relatively narrow band illustrates in a striking way the unusual fea- 
tures presented in supernovae. 
THE SUPERNOVA IN NGC 1003 

The spiral nebula NGC 1003.—The nebula (R.A. 2"35™3, Dec. 
+ 40°36’, 1938; gal. long. 112°1, lat. —16°7) is an Sc spiral of 4/3 di- 
ameter, seen edgewise. The brightest stars of the system, which have 
the mean photographic magnitude m, = 20.5, are just above the 
limiting magnitude of 1-hour exposures with the roo-inch telescope. 
Adopting M, = — 6.3 for the corresponding absolute magnitude, we 
find the modulus of NGC 1003 to be m — M = 26.8. The integrated 
apparent photographic magnitude of the nebula is mye = 13.1,° 
which gives M,.. = — 13.7 for the absolute magnitude of the system, 
slightly below the mean value My nev = —14.2 for the average 
extragalactic system. The distance and the linear diameter of NGC 
1003, derived from the modulus m — M = 26.8, are 2.29: 10° and 
2980 parsecs, respectively. In view of the low latitude of the nebula, 
these figures should be corrected, however, for the general galactic 
absorption. The factor by which the uncorrected values of the dis- 
tance are to be multiplied is, according to Hubble,’ 


log f = —o.05|cosec b|, 


or, in the present case, f = 0.67. Thus, we have the following data 


for NGC 1003: 
Modulus, 26.8 , 


Integrated absolute magnitude, — 13.7 (photographic), 


Distance 1.51 + 10° parsecs 
Diameter 2000 


corrected for latitude effect . 


We use intentionally the uncorrected value of the distance modulus. 
Since in the following discussions only magnitude differences are 
used, the latitude effect is unessential. 

Photographic magnitudes of the comparison stars.—The sequence of 
our comparison stars is marked on Plate XIX. Their photographic 

9 Mt. W. Contr., No. 600; Ap. J., 88, 285, 1938. 

to Mt. W. Contr., No. 485; Ap. J., 79, 8, 1934. 
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COMPARISON STARS FOR SUPERNOVA IN NGC 10032 


Scale, 1 mm 8"o 
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magnitudes have been derived from comparisons with S.A. 22 and 
S.A. 68 made at the Mount Wilson reflectors. For stars fainter than 
17.7 the scale is that of S.A. 68, for which photographic magnitudes 
down to Mpg = 21.0 are now available. Since careful checks have 
shown that the magnitude system of S.A. 68 agrees in zero point and 


TABLE 3 
SEQUENCE FOR SUPERNOVA IN NGC 1003 


(Unit for Residuals, 0.01 Mag.) 























RESIDUALS FOR INDIVIDUAL PLATES 
Mrzan j|_._._- Se ee ee 
STAR _ an ae ] 
Mog | | 
608 | 609 | 616 | 874 875 | 629 | 633 
\ £2.50 | 
1 t3.12 +12 | io} —1| 
b 13.22] +18 | —13] 4 | 
¢ 13.31 +23 i eons = 1 
d 13.82 + 9 |} —8] —1 + 2} 
€ 14.46 + 9 — 6] — 4] +3 ear 
| | 
f 14.74 | +10 sel alge aes — 4] 
g 15.46) +1 — 5| “~21 SIG) Se 
h 15.55 a; 2] Sg 
i 16.00 + 5 — | — ] — J — 6 + I 
k 16.22 + 9 —I1I = aa =26 
] 16.56 + 7 ° Lay —= $4 9 
m C7542 + 5 aa ra ° 
n 17.50 — 6] 3 EA “iy 
0 17.80 — 23 +12 — 9 +18 
p 18.22 f +4 —4 
| 
q 18.54 31 Fy 
r 18.81 a = § 
S 19.25 4 <9 
t 19.82 —12 +12 
u 20.03 | 











scale with that of the Polar Sequence and since there is no percepti- 
ble systematic difference between the magnitudes of S.A. 22 and 
S.A. 68, the magnitudes of our comparison stars should be on the in- 
ternational system. Their values are given in Table 3, together with 
the residuals for the individual determinations. 

Photographic light-curve of the supernova.—Our measurements of 
the supernova, from the date of its discovery to the time when ap- 
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proaching conjunction with the sun prevented further observations, 
are given in Table 4. 
TABLE 4 


PHOTOGRAPHIC MAGNITUDES OF SUPERNOVA IN NGC 1003 






























































J.D. G.M.T Mog Instr. 5D. GMT. Mog Instr 
2428787 .8 12.98 i 2428808 .9 14.27 S 
788.8 13.01 S 812.7 14.68 S 
788.9 13.52 100 812.9 14.61 60 
789.8 13.16 S 813.9 14.61 60 
790.8 ; 12.85 5 814.0 14 64 60 
790.9 12.84 5 814.7 TAF S 
791.9 12.67 S 820.8 15.03 S 
791.9 12.89 5 838.9 15.80 100 
792.9 12.90 S 838.9 15.68 100 
793-9 } 3.5% 5 840.8 15.69 60 
7904.9 13.00 = 840.9 ee ke 60 
795-9 13.00 S 863.7 16.42 60 
799.7 | 13.19: S 865.7 16.47 60 
799.7 | 12.93: S 872.6 10.55 60 
800.7 13.10 S 895.6 16.93 60 
801.7 13.39 S 895.8 16.95 60 
801.8 13.20 5 897.8 16.91 100 
802.7 13.69 S 898 .7 16.99 100 
803.7 13.66 S 929.8 17.05 100 
804.8 13.79 S 954.6 7.41 100 
805.8 13.90 S 955.6 17.42 100 
807.9 14.26 S 990.6 17.43 100 
808.8 rA..27 S 
* See note to Table 2. 
m M 
12 —14.8 
14 —12.5 
16 —10.8 
4 a) a a —i— 8.8 
JD2428780 820 860 goo 940 980 


Fic. 2.—Photographic light-curve of supernova in NGC 1003 
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The resulting light-curve is reproduced in Figure 2. Inciuded in 
this curve are the following two observations by A. A. Wachmann"™ 
which precede our earliest plates: 


JD 2428782.6G.M.T. 13.73 
784.5 G.M.T. 13.28 


Since Wachmann’s magnitude system differs systematically from 
ours by +0.16 mag., his published magnitudes have been diminished 
by this amount. The light-curve shows the importance of Wach- 
mann’s early observations, which clearly belong, as he already point- 
ed out himself, to the ascending branch. The supernova had a rather 
extended maximum lasting 10-12 days, during which minor fluctua- 
tions are indicated. The decline in brightness following the maxi- 
mum is closely similar to that of the supernova in IC 4182. Since the 
apparent brightness of the maximum was Mmax = 12.8, we obtain 
for the absolute photographic brightness of the supernova at maxi- 
mum 
Mimax = —14.0, 


a value practically identical with the mean, Mma = —14.3, for 
supernovae. 


Further discussion of the two light-curves is postponed until the 
observations of another season and the light-curves in the red region 
of the spectrum are available. 

CARNEGIE INSTITUTION OF WASHINGTON 
Mount WILSON OBSERVATORY 


June 1938 


1 Beob. Zirk., 19, No. 42, 1937. 








PHYSICAL PROCESSES IN GASEOUS NEBULAE 
V. ELECTRON TEMPERATURES 


JAMES G. BAKER,’ DONALD H. MENZEL, AND 
LAWRENCE H. ALLER 


ABSTRACT 

The relation between the electron temperature of a gaseous nebula and the temp- 
erature of the exciting star is derived and applied to an optically thin nebula. The 
b,’s have been recomputed both with the radiation field included and as a function of 
the temperature of the exciting star. The Balmer decrement proves to be steeper when 
the influence of Lyman line absorption is included. 

The equations derived in papers I, II, and IV of this series may be 
applied to the problem of calculating the electron temperature of the 
nebular gas. We make the following assumptions: the electrons are 
produced by photo-ionization; radiationless encounters between the 
various components of the assembly are sufficiently frequent that 
the electrons attain a Maxwellian velocity distribution; the nebula 
is in a steady state with respect to the distribution of atoms in vari- 
ous levels; the energy emitted per unit time must equal the energy 


absorbed. 
The calculations of paper III gave the values of the quantities 
b,,b;,.... etc. The results are based on the equations for statistical 


equilibrium of quantum levels 2 and higher. These equations are 
sufficient to determine the absolute nebular emission and the Balmer 
decrement when JN, and 7, are given; we have still to predict the 
values of the two latter quantities from the physical parameters of 
the exciting radiation. In addition, there is the third parameter, 0, 

which defines the number of atoms in the ground level—to be de- 
termined. We may suppose, for the sake of simplicity, that the radia- 
tion from the nuclear star is “‘black,” at temperature 7,, with the 
dilution factor W specified. There are still two equations from paper 
I or paper IV that have not been employed, viz., the one defining 
the statistical equilibrium of level 1 and the one expressing conserva- 
tion of emitted and absorbed radiation. A third relationship is ob- 

* Society of Fellows, Harvard University. 
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tained when the density of the nebula is given, hence the problem 
is determinate. One may substitute the previously evaluated b,’s in 
the two remaining equations and calculate 6,, V., and T,. 

Fortunately, for the gaseous nebulae, where the radiation is very 
dilute one may determine 7, directly without first evaluating the 
intermediate b,’s. T, proves to be independent either of the density 
or of the dilution factor, as long as the latter is small. The simplifica- 
tion is obtained by means of the mechanistic approach developed in 
paper IV. The numerical results are given for an optically thin 
nebula, by which term we mean a nebula of optical thickness small 
enough that the spectral energy curve of the star, beyond the Lyman 
limit, is not seriously depleted by absorption effects. 

It now seems advisable to write out the necessary equations in the 
form in which they are to be used. We shall use the equations and 
the notation of papers I and II,’ writing out first the explicit forms 
of the equations IV (3), IV (4), and IV (9). Adopting as before 
b, 1, we find that IV (3) becomes, in view of the definitions of 


II and equations IV (14), I (25), 


oN shy, /kT 1 ‘ S ie, s—h(v-vy)/kT 1 
: em € dv Vion s ¢ dv 
bf WW a ——— == +f Wel Soe J) —- 2) 
SJ" Gree == ee aw Six Cn ee 
I 


It must be borne in mind that W,, is a very small quantity of the 
order of to '°. The second term on the right of the foregoing equa- 
tion represents the negative absorptions. It follows that 


bie /kT Wye /*T dy ehr/kTe pee hy/kT dv 
erin /kle a ee et A sees W ye he/ BF e — 
Jv, CT — 1 y eal/tt: — 1 J, y 


kT. kT, 
hy, aamee hy, ees ‘ 
< oo iT. — ; Wa < ohn [kT a Wa S&S 3Wu- 


Therefore, the effect of induced emissions, i.e., induced captures, is 
negligibly small. From (1) we see that 6, is a very large quantity 
determined principally by the magnitude of W,,. Terms that con- 

2D. H. Menzel, Ap. J., 85, 330, 1937; D. H. Menzel and J. G. Baker, Ap. J., 86, 


70, 1937; J. G. Baker and D. H. Menzel, Ap. J., 88, 52, 1938; D. H. Menzel, L. H. 
Aller, and J. G. Baker, Ap. J., 88, 313, 1938. 
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tain both W,, and 0, as a product cannot be neglected, but the terms 
that contain W,, without the 0, factor can readily be neglected. This 
amounts to a formal proof of our previous statement that the nega- 
tive absorptions from all levels and the positive absorptions from all 
levels but the first are negligible. Our equations may now be writ- 
ten 


bo 
— 


; rae dy —S, : 
ek W, ~~ « SS a Ge, ( 
Jy: Vey — tT) : 


ee ae 
I 


which defines Gr,. The equation of radiative equilibrium for the 
assembly is 


T’2 k? va Pr ~ I I | 
t he (=z) + akT, ot 2hrZ — eae” ins 


Tl 2hRZ? I ee yee dy | 
= b,e*: > WW in” > os Te beXikT, i a ‘ = Sa: I 
hy/kT 1 
manent e + a J "alt | 
2 


where 


The equation of radiative equilibrium for the continuum becomes 


k ee ( : 
»X y £ = 2 
ee (az) rf : “ev — 1 el i 
k aes &*S. 
r ee ae ioe, Te ns’ 


Equations IV (1) and IV (2) may be discussed in a similar man- 
ner, and we shall treat the problem later in connection with the der- 
ivation of the formula for b,. It is to be pointed out that all these 
equations contain b,, T,, and 7,. One might expect, therefore, that it 
would be possible to choose any two of them, eliminate ,, and derive 
a relation between 7, and 7,. In particular, one might suppose that 
we could choose any two of equations (2), (3), or (4), solve each for 
b,, and equate. Our choice is really somewhat more limited than 
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that. If we suppose, for example, the star to radiate as a black body 
with no absorption lines in the Lyman series, we do not know the 
b,’s and therefore cannot use equation (3). Our procedure in this 
case would be to solve (2) and (4) for b, and equate. In the problem 
where we supposed no radiation to leave the star in the Lyman lines 
but where the star is supposed to shine as a black body beyond the 
Lyman limit, we know the 0,’s and can therefore apply any two of 
these three equations. Numerical calculations based on (2) and (3) 
gave exactly the same relation between 7, and T, as those based on 
(2) and (4). The agreement provides an excellent check on the for- 
mulas and computations. Note that the relation between 7, and 
T, does not depend on whether the central star has any line radia- 
tion or not. It is to be mentioned that the expressions for b, from 
(2) and (4) are equally good for a static or an expanding nebula 
because absorption in the continuum is not seriously affected by a 
Doppler shift. From (2) and (4) we have for the optically thin 
nebula, where W,, does not depend on », 


( k )r { dy 
2hRZ?) "Jy, ev = | 
{ dy | 
Py MO = 1) | 


i Vey cf lee ae, ee 
(ae) fro (sa) LT. + 2G7,— 2 : 





or, symbolically, 


since the left side of the equation depends only on T, and the right 
side only on 7,. For any assigned T, the corresponding 7, may be 
readily found. Each side of the equation is plotted as a function of 
its appropriate temperature. A specified value of 7, determines 
f.(T,), which must equal f,(7.). From the graph of f,(7.) the appro- 
priate 7, may then be ascertained. It is to be mentioned that the 
quadratic f, — f, gives but one real positive root. 
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The relation between T, and 7, is given in the following table for 
the optically thin nebula. The Gaunt g factor was taken equal to 
unity in the calculations. 


T; Te 
MOM Sh crenatinu trek Bit Aoree 5,000 
FOROOO occ Maribel thie’ derine “OG5CO 
ROFOOO < ...54s hs Eo . 18,000 
BOUOSO os. dlehaicca dake eeatae BaG6o 
BOVOGO: 6s sues Neigh dai wtewe. sO0o 
MOOMOOO so 5. de acndinisie as ue, O2EOOO 
BIO OO! «5 clit. dain dled ce 1EG2ZO0O 


Note that for the higher temperatures the electron temperature de- 
parts farther and farther from the temperature of the source. For 
very high temperatures the effect of the neglect of the g factor pro- 
duces an error of at most 5 per cent in the foregoing results. At low 
temperatures the error is negligible. 

From equation (2) or (4) we can determine b,. The fact that its 
order of magnitude must be determined by the dilution factor fol- 
lows at once from elementary considerations. Comparing the gaseous 
nebula with an assembly in thermodynamic equilibrium at tempera- 
ture 7., we see that the captures must go at essentially the same rate 
but that the photo-ionizations are fewer by a factor of about IW ,,. 
As a result, there appears an enormous excess of atoms in the ground 
state as compared with the assembly in thermodynamic equilibrium. 

The ionization equation for hydrogen follows at once. We shall 
defer discussion, however, because ionization is unimportant so far 
as the present investigation is concerned. 

We now write equation IV (2) in terms of the physical quantities 
involved. In view of I (7), and using the notation of II (3), we find, 
neglecting the terms multiplied by W and not by 0,W, 


a) . | 
byveXn Bn’ nUn''n + gS n 
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TABLE 1 


-LAR TEMPERATURE, ELECTRON TEMPERATURE, AND ASSOCIATED 
VALUES OF bn 
































T; 
5,000 10,000 20,000 | 40,000 80,000° 160,000° 
Te 
- | 
5,000 9,500 18,000 | 4,000° 57,000° 92,000° 
| 
b» (Neglecting Radiation Field) 
| | as 
O.00003 0.0029 0.027 0.090 | O.220 | 0.388 
0033 025 forei6) 105 324 471 
o1s of 8 I5I 270 392 526 
043 109 210 324 444 . 500 
072 150 251 366 482 585 
000 182 290 400 511 605 
126 213 319 425 531 621 
150 237 341 449 549 639 
171 201 303 4609 563 6050 
1QI 281 381 481 576 .659 
209 299 396 .492 .588 666 
22: 313 409 503 598 .073 
238 328 420 511 604 679 
0.251 0. 339 0.430 0.520 0.610 0.685 
b, (Including Radiation Field) 
P 268 2.567 2.369 | 2.281 2.181 2.3093 
| S592 1.608 1.603 1.561 I.429 1.551 
| 1.48 1.414 1.411 1.367 1.258 1.356 
| 1.40 E24 "433 1.29 1.19 ey 
| E37 1.29 5.27 t. 25 1.16 r.22 
1.36 E25 E.25 Tu23 I.14 1.19 
| 1.34 £29 1.22 I.20 1.52 B49 
bh. 32 20 I.20 1.19 r.21 1.16 
E. 3y 1.18 1.18 1.18 E16 1.15 
1.28 E257 1 ee 1.17 1.09 1.14 
1.260 r.55 1.10 1.10 I.09 I.13 
i323 I.14 $25 es | 1.08 1.13 ° 
1.20 r.33 1.14 I.14 I.08 I.12 
1.19 E32 E33 I.14 I.07 I.11 
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This equation may be solved in a manner similar to the method ex- 
plained in paper I, once the W,, factors are known. We have 





Sar V 
b, = es ii IE I : (7 
7 ex nt, 7) 
where 
; : bye/FTe \ anteXn 
has — Me + WW In (=! 2 aS | 
etv/KIT; — ~} yn? — 1 | 
2n? | | 
| G =i » (8) 
= 1S, + Win "et, jie | 
ais W xx { dy (eM/kT:1 — 7) 
. Ser >. x) 





For an optically thin nebula surrounding a star radiating as a 
black body, W,, = W,,. When we consider the transfer problem, 
this relation will no longer hold. V,, is defined in terms of S,, in the 
same manner that V,, was defined in terms of Sy). 

The mathematical grouping of these terms must not be interpreted 
as due to their physical similarity. S, arises from captures on the 
level in question; the W,, term originates from transitions from the 
ground level to the mth level. 

If the star is assumed to radiate as a black body, the b,,’s are con- 
siderably different from the 0,’s discussed in paper III. For small 
values of they are larger than unity and gradually approach 1 as 
n—>o,as indeed they must. For purposes of comparison, 0, is tabu- 
lated as a function of the temperature of the central star for the two 
cases considered. 

Computation (Table 1) shows the Balmer decrement to be very 
insensitive to the excitation temperature. It should be pointed out 
that the optically thin nebula here considered will have a brilliance 
much less than that of the nuclear star, so that the present results 
are not directly applicable to planetary nebulae. But the enormous 
effect of the nuclear radiation suggests that even in an optically 
thick nebula the nuclear radiation might show some influence on the 
Balmer decrement, which would tend to be steeper in consequence. 


In the next paper the theory of radiative transfer for the Lyman 
levels and continuum will be discussed. 
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THE MAGNITUDE OF THE SUN, THE STELLAR 
TEMPERATURE SCALE, AND BOLOMETRIC 
CORRECTIONS 


G. P. KUIPER 


ABSTRACT 


Three connected subjects are treated in this paper. The first section deals with the 
effective temperature, the absolute photovisual and bolometric magnitudes, and the 
spectral type of the sun. The quantities adopted are: 7, = 5713° (Unséld); Mpv = 4.73; 
Mvpol = 4.62; spectrum dGo2. 

The second section contains a discussion of 7, of the O, B, and A stars. The result 
is given in Table 3 and shown in Figure r. 

In the third section, tables of bolometric correction are constructed, as far as pos- 
sible independently of assumptions about the shape of the spectral energy-curves. For 
spectral types earlier than Fo, theoretical spectral energy-curves are used; the results 
are found in Tables 6 and 3. For later types, empirical bolometric corrections are de- 
rived from Pettit and Nicholson’s radiometric observations. The discussion is based on 
Table 7 and Figures 2 and 3; the result is given in Tables 9 and to. 

In the fourth section, effective temperatures are derived for late-type stars with 
suitable data: Castor C and the giants for which diameters were measured interfero- 
metrically. Limb darkening has been taken into account. The result for Castor C is 
3550 + 110°; for the giants, the values in Table rr. 

In the last section, color temperatures are collected and reduced to the system of 
effective temperatures. Particular attention is paid to the observations of the dwarfs 
(Tables 12 and 14), on which the reduction to effective temperatures depends to a large 
extent (sun). Table 13, last column, gives the derived effective temperatures for classes 
A-M;; they are shown graphically in Figure 4. Directions are given for the calibration 
of temperatures derived from W. Becker’s Catalogue of color indices and Pettit and 
Nicholson’s water-cell absorptions, to effective temperatures (p. 469). Two applications 
are given. The result on e Aurigae reduces a difficulty previously encountered by B. 
Stromgren. 

INTRODUCTION 

The purpose of this article is the selection or derivation of the 
most probable set of effective temperatures and of bolometric cor- 
rections that can now be obtained. These data are required in the 
discussion of the mass-luminosity relation given in the following 
paper. Only a few new results are derived here, but the main object 
is to select and relate the most reliable observations and results thus 
far obtained. 

Among the recent investigations of the temperature scale, those 
by Pannekoek! and Pilowski’ are of special interest. The results of 
Pannekoek’s paper are extensively used in our article. Pilowski’s pa- 


1 Ap. J., 84, 481, 1936. 2Zs.f. Ap., 11, 265, 1936. 
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per is an attempt to use, for main-sequence stars only, the statistical 
relation between radius and spectral type, and between spectral 
type and absolute magnitude, from which the relation between spec- 
tral type and surface brightness is derived. The first relation is 
found from spectroscopic binaries; the last, from stars in general. 
But as Pannekoek has pointed out,’ the relation between spectral 
type and temperature is much closer than the statistical relations 
from which this relation is derived by Pilowski, so that it seems pref- 
erable to restrict the discussion to those stars for which all the data 
are well determined. Furthermore, Pilowski derived ‘‘radiation tem- 
peratures,”’ as defined by Brill, and not effective temperatures (“‘ra- 
diation temperatures” are formal temperatures related to the visual 
surface brightnesses by Planck’s formula; since the stars do not ex- 
actly radiate as black bodies, the “‘radiation temperatures’? depend 
on the wave length used in measuring the surface brightness). It 
will be our aim to obtain the more fundamental effective tempera- 


tures, defined by 
L = 4rR’- oT, (1) 


in which L is the luminosity and R the radius of a star, and o the 
radiation constant. If L and R are expressed in terms of the sun, we 


. ; 7 4 
L =o (7) ‘ (2) 


If L and R are known, 7, may be found from equation (2). But 
since L is never known directly from observation, but usually only 
the absolute visual (or photovisual) magnitude, equation (2) can 
be used in practice only if a table of bolometric corrections is avail- 


have 





able. The bolometric correction is defined by 


B.C. = Mtoi — Mpy = mbo1 — Mpy , (3) 


Mbo1 = Mooi (Sun) — 2.5 log L . (4) 


3 Op. cit., p. 505. 4 E.g. Pannekoek, of. cit., p. 483. 
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The definition of the quantity M,,, (sun) is given in the course of 
this paper. It involves the zero point of the table of bolometric 
corrections. The photovisual magnitudes (m,,) used in this paper 
are international photovisual magnitudes (IPv) on the system of the 
North Polar Sequence of 1922. Other visual magnitudes are reduced 
to this system before they are used as is explained on p. 477. 

In the following sections it will appear that a satisfactory table 
of bolometric corrections can be constructed without recourse to the 
assum ption of black-body radiation, and that also the scale of stellar 
effective temperatures can be brought into a fairly satisfactory 
state. 

THE TEMPERATURE AND THE LUMINOSITY OF THE SUN 

1. The effective temperature of the sun is one of the very few 
absolute determinations of stellar temperatures which may be ob- 
tained. The luminosity, expressed in stellar magnitudes, is of im- 
portance in fixing the unit of stellar luminosities. 

The solar effective temperature is derived from the solar constant; 
since the atmospheric extension depends on the wave length, the de- 
termination of this constant involves a spectrophotometric study 
of the whole accessible solar spectrum. Unséld has just published’ 
a careful rediscussion of the available data, and we adopt his result: 


T. = 5713 + 30° K (estimated mean error 
O/435 J ’ 


or 
log f Pa — +. 969 + 0.002 (m.e.). (5) 


2. The stellar magnitude of the sun is usually taken from Russell’s 
discussion,® which gave my = —26.72 and M, = +4.85. In the 
present paper, as in the next one, an attempt is made to reduce 
all magnitudes to the IPv system. Accordingly, in the comparisons 
tabulated by Russell IPv magnitudes will be used: 

a) Zoéllner found the sun 26.87 + 0.05 mag. brighter than Capella, 
or — 26.68 IPv. 

b) Fabry found it to be 26.94 mag. brighter than Vega, or — 26.85 
IPv. 

5 Physik der Sternatmos phiren, chap. ti, 1938. 


°Ap. J., 43) 103, 1916. 
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c) Ceraski found it to be 
22.93—0.66 + 4.86 mag. brighter than Procyon, 
22.93——0.66 + 6.30 mag. brighter than Polaris, and 
22.93 —0.66 + 3.22 mag. brighter than Sirius. 
With weights 4, 4, and 1, as used by Russell, Ceraski’s value is 
— 26.64 IPv. 

d) W.H. Pickering made to determinations using different stars, 
which give the average — 26.93 IPv. 

The straight mean of these four visual determinations is — 26.77 
IP. 

No direct visual determinations of the magnitude of the sun have 
been published since Russell’s discussion. 

3. Two photographic determinations are available, already men- 
tioned in Russell’s paper. Since the dwarf F and G stars have colors 
almost uniquely determined by their spectral type,’ it appears safe 
to adopt the color index corresponding to the spectral type. This 
type, as determined by Morgan and by the writer, is Gz. The cor- 
responding color index is*® + 0.53. 

E. S. King, from 1o determinations, found m,, = —25.83. Six 
of these depend on comparisons with Polaris, for which King used 
his value m,, = 2.69, as compared with 2.55 IPg. Applying the 
same correction, —0o.14 mag., to the other 4 comparisons (since the 
magnitudes of the comparison stars were determined from direct 
comparisons with Polaris), King’s determination becomes — 25.97 
IPg, or — 26.50 IPv. ; 

The second photographic determination is by Birck. The result, 
as reduced by Russell to King’s photographic system, is — 26.12, to 
which the correction —o.14 mag. should apply as mentioned above. 
Accordingly, we have m,, = — 26.26 IPg or —26.79 IPv. Giving 
double weight to King’s determination, we obtain — 26.60 I Pv from 
the two photographic determinations. 

The preceding values may be supplemented by two other deter- 
minations, one which uses the moon as an intermediate step, and the 
other depending on Pettit and Nicholson’s radiometric work. 

The use of the moon as an intermediate step involves, first, the 
determination of the magnitude difference between the sun and 


7 Cf. Morgan, zbid., 87, 460, 1938. 8 Tbid., 86, 180, 1937. 
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moon, and second, the determination of the stellar magnitude of the 
moon. 

4. A very fine determination of the magnitude difference between 
sun and moon has recently been published by Rougier.? He used for 
these measurements a photoelectric (potassium) cell. He observed 
the sun through an optical device which reduced the brightness by 
14.87 mag., and observed the moon without that device. This made 
the intensities to be measured of the same order of magnitude. The 
reduction was obtained by placing two microscope objectives and 
some diaphragms at proper distances. The reduction constant fol- 
lowed from geometrical optics and from an empirical determination 
of the absorption by the objectives. The uncertainty in the reduc- 
tion constant was estimated to be +o.o1 mag. Both the sun and the 
moon were compared with a standard lamp, which provided a pro- 
visional zero point of the magnitude scale. With this zero point, the 
magnitude of the sun was found to be —17.33+0.03; of the full 
moon, — 3.04; and of the moon at other phases, the values given in 
Rougier’s Table 30." 

Since the relative magnitudes in Rougier’s system seem to be very 
well established, the calibration to stellar magnitudes by means of 
the moon can be made from magnitude determinations at any lunar 
phase covered by Rougier. 

Rougier’s magnitude difference may be combined with 5 magni- 
tude determinations of the moon. 

a) King’s photographic measures," 83 in all, were made by indi- 
vidual comparisons with Polaris. In principle, therefore, one ob- 
servation should suffice to fix Rougier’s zero point, and the whole of 
King’s work should give a good determination if no serious system- 
atic errors are present. 

If all of King’s observations are used, it is found” that the magni- 
tude of the full moon is —11.64 on King’s scale. If only the part 
between the first and last quarter is used, when the observations by 
both authors could be made at greater altitudes (63 observations), 
we find —11.58 + 0.05 (p.e.) on King’s scale. Adopting the latter 

9 Ann. Strasbourg, 2, Part 3, 1933; 3, Part 5, 1937. 

10 [bid., 25 310. 

" Harvard Ann., 59, 63-94, 1912. 12 Rougier, op. vil., p. 323. 
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value, and the IPg magnitude for Polaris (2.55 instead of King’s 
value 2.62), we obtain 


m(full moon) = —11.65 IPg. 


b) Two visual determinations of the stellar magnitude of the 
moon are available, the first of which is by Zollner. (Two still older 
determinations will not be used here.) Z6llner found’ the full moon 
to be 12.31 mag. brighter than Capella. Reducing his observations to 
full moon with the aid of Rougier’s light-curve instead of Zéllner’s 
own curve, we find 12.48 mag. for the difference. Hence, we have 


m(full moon) = —12.29 IPv (Zollner) . 


c) The second visual determination is by W. H. Pickering.’ He 
compared both the moon and Arcturus with a pentane lamp, the 
color of which was made equal to these objects by means of a blue 
solution acting as a filter. The scale was exclusively based on the 
inverse-square law. Reducing the observations® to full moon with 
the aid of Rougier’s light-curve, and using only the observations 
made between first and last quarter,’ we find for the full moon 
—12.63 on Pickering’s scale (with Arcturus 0.24 mag.). For the 
IPv magnitude of Arcturus we find +0.04, so that 


m(full moon),= —12.83 IPv (Pickering) . 


This result is based on 18 observations of the moon but on only 4 of 
Arcturus. The internal probable error of the mean of the observa- 
tions of Arcturus alone is +o0.11 mag. But Pickering’s observations 
are of interest because of their extreme simplicity. 

The agreement between the two visual determinations is rather 
poor. Giving Pickering’s result double the weight of Zéllner’s, the 
combined visual result is 


m(full moon) = —12.65 IPv. 


"3 Photometrische Untersuchungen, pp. 109 and 124, 1865. 
14 Harvard Ann., 62, 62-69, 1908. 5 [bid., Table XI, p. 65. 
© Two observations, Nos. g and 10, were omitted because of Pickering’s comments 


on them. 
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d) Pettit’s radiometric work on the moon as a whole” gives in- 
cidentally a determination of its visual magnitude.” 

Unfortunately, Rougier’s light-curve and Pettit’s curve of the re- 
flected light are not in very good agreement.’? It seems probable that 
at least part of the discrepancy can be attributed to the extrapola- 
tion of the light-curve to full moon, which is necessarily somewhat 
arbitrary, particularly when few observations are made very near 
full moon. Although color determinations of the moon show no 
change during the lunation, only a small variation could appreciably 
affect the difference between the photoelectric and the radiometric 
observations, in view of the yellow color of reflected moonlight. 

On account of this difference, Pettit’s determination will give a 
lower limit to the brightness of the full moon as compared with Rou- 
gier’s system, which we have used in the reductions. (It is clear that 
this fact arises only because we use the magnitude of the full moon 
as derived by Pettit from his own curve. In our previous discussion 
we use Rougier’s definition of the full moon only as a constant con- 
necting the solar magnitude with all the lunar observations between 
first and last quarter. Hence, it does not concern us whether or not 
Rougier’s definition is exactly correct.) 

Pettit’s computed value for the brightness of the full moon is 
— 12.63 mag., which is on the Harvard scale. Since this value is a 
lower limit in the system used here, we shall adopt — 12.75 IPv as 
a probable value. 

e) Calder?® compared the moon photoelectrically with Capella, 
Vega, and Deneb, with a standard lamp as intermediary. In order 
to obtain the photoelectric magnitude of the lamp on the interna- 
tional scale, we apply 0.8” of the color index (IPg—IPv) to the IPv 
magnitudes. With the directions of page 478 we find for the IPv mag- 
nitudes of Capella, Vega, and Deneb, + 0.19, + 0.09, and +1.34, 
respectively; and hence, with the color indices given elsewhere” we 
find the photoelectric magnitudes to be +0.75, + 0.09, and + 1.39. 
These values, in connection with the differences found by Calder, 

17 Ap. J., 81, 17, 1935; Contr. Mt. W. Obs., No. 504. 

8 Ap. J., 81, 33-34. 

19 Ibid., Fig. 4, p. 28. 20 Harvard Ann., 105, 445, 1937. 

21 E.g., Stebbins and Whitford, Ap. J., 87, 252, 1938. 22 Ibid., 86, 180, 1937. 
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lead to — 4.33, —4.35, and —4.33 mag. for the lamp, or — 4.34 mag. 
in the mean. 

Calder measured the moon on two nights, one of which led to final 
results. The phase angle during these observations was +5°50’, so 
that the correction to full moon*s is —o.12 mag. (Calder, with a 
somewhat different definition of full moon, obtained —o.18 mag.). 
Calder’s data then make the full moon 7.75 mag. brighter than the 
standard lamp, or —12.09 IPe. 

We must now reduce the five different values for the magnitude 
of the moon to one system, for which we adopt the IPv system. To 
that end the color index of the moon is needed. The color index of 
the sun has been assumed to be + 0.53 mag.”4 The difference be- 
tween the color indices of sun and moon is found from the following 
measures (the moon is the redder): 


Wilsing, Scheiner, Russell* 0.25 mag. 
Pettitt........ 20 
Danjont. erry 33 
Calder§.... heals pees @289 
Weighted mean ©.22 mag. (6) 


* Ap. J., 43, 126, 1916. 

t Pub. A.S.P., 38, 243, 1926. 

t Ann. Strasbourg, 3, 146, 274, 1936-37. 

§ Harvard Ann., 105, 452, 1937. The color index was found 
to be equal to that of Capella, and hence 0.17 mag. larger than 
for the sun. 


Pettit’s determination has been given double weight and Danjon’s 
visual determination haJf-weight. 

Accordingly, the color index of the moon is + 0.75 mag., corre- 
sponding to dKo, in agreement with Pettit’s result. 

The five magnitude determinations of the full moon are, therefore, 


—11.65 IPg = — 12.40 IPv (King) 
12.65 (Zéllner, W. H. Pickering) 
12.75 (Pettit) 
—12.09 IPe = — 12.69 IPv (Calder) 
Mean = — 12.661Pv (7) 


In forming the mean, double weight was given to the two modern 
determinations. It should be remarked that Pettit’s own reduction, 


23 Rougier, op. cit., Table 30. 24 Cf. p. 432. 
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— 12.63 HS, is very close to the average; but we felt that a small cor- 
rection was needed to allow for the difference between Rougier’s and 
Pettit’s light-curve. 

The value — 12.66 IPv here obtained does not appreciably modify 
Russell’s result — 12.55 HS; the difference is about equal to the some- 
what greater relative brightness of the full moon as defined by Rou- 
gier’s light-curve. 

This result for the brightness of the full moon, together with 
Rougier’s photoelectric difference, 14.29 mag., gives 


m(sun) = —12.66 — (14.29 — 0.8 X 0.22) = —26.77 IPv. 


5. Pettit and Nicholson’s paper on “Stellar Radiation Measure- 
ments’’?> contains a calibration of radiometric magnitudes in terms 
of energy, from which the authors obtained the radiometric mag- 
nitude of the sun, —27.18 mag. From this figure the bolometric 
magnitude is found to be — 27.18 —Am, + 0.62 = —27.01 mag. 
The constant 0.62 mag. is explained on page 451; it reduces the 
quantity (m,—Am,) to the zero point of bolometric magnitudes used 
in this paper.” The average bolometric correction for a G2 dwarf is, 
then, 0.07, as found below (Table 9). Hence, the resulting visual 
magnitude of the sun is — 26.94 IPv. 

The uncertainty of the bolometric magnitude of the sun derived 
from Pettit and Nicholson’s measures is probably less than 0.05 
mag. But, as we shall see later, there may be some scatter in the 
heat indices for main-sequence stars, which also enters in the empiri- 
cally determined bolometric corrections. Hence, the quantity 0.07 
may be in error by a small amount. 

The four groups of determinations of the IPv magnitude of the 
sun thus far discussed are: 


Four visual determinations...................... —26.77 IPv, weight 1 
Two photographic determinations 26.60 ; 
Sun-moon, plus 5 determinations of moon 26.77 2 
Radiometric determination... .... ae — 26.04 3 

Weighted mean............. ee . —26.84+0.04 (p.e) IPv 


2 Ap. J., 68, 279, 1928; Contr. Mt. W. Obs., No. 369. 
> The quantity Am; includes the zenith absorption of the atmosphere as well as the 
loss in the telescope (Pettit and Nicholson). 
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One observation of the magnitude of the sun remains to be dis- 
cussed. Calder*”® compared the sun on three days with his standard 
lamp, which, as we have seen, was in turn compared with three 
stars. The photoelectric magnitude of the sun was found to be 21.46 
mag. brighter than the standard lamp, for which we obtained — 4.34 
IPe. Hence m(sun) = —25.80 IPe = — 26.22 IPv. 

This value deviates 0.6 mag. from the mean value obtained above. 
The photoelectric method used by Calder is, of course, capable of 
yielding very accurate results. The difference cannot, therefore, be 
ascribed to an accidental error in Calder’s work, nor can our mean 
value be in error by an amount much larger than o.1 or 0.2 mag. We 
conclude that Calder’s observation may be affected by a systematic 
error and that it seems best to omit it at the present time. 

In order to obtain the most probable value of the bolometric 
magnitude of the sun, we combine the result obtained from Pettit 
and Nicholson’s radiometric work, — 27.01 mag., with the weighted 
mean, — 26.75 IPv, of the first three groups of observations given 
above. This mean leads to m,,; = — 26.82. On account of the un- 
certainty in the bolometric correction, the relative weights have to 
be different from those assigned in forming the mean IPv value. We 
raise the weight of the radiometric determination from 3 to 5 and 
lower the total weight of the other determinations from 33 to 25. 
The weighted mean then becomes 


Myo1 = — 26.95, Mooi = +4.62, (8) 


whereas 
Mpy = — 26.84, Mpy = +4.73.- (9) 


The third quantity needed for the sun is its spectral class. From 
low-dispersion slit spectra Dr. Morgan and the writer independently 
found G2 on the conventional scale. 


SUNSPOTS 
As mentioned below, only for three stars in addition to the sun 
have radii and luminosities been determined with such accuracy as is 
needed for furnishing reliable determinations of effective tempera- 
ture. For that reason it is fortunate that the sun provides a second 
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determination for a type much later than G2. Since Minnaert and 
Wanders have shown that the umbra of sunspots is in radiative 
equilibrium,” it follows that the relation between surface brightness 
and spectral type will be the same as for ordinary stars in radiative 
equilibrium. Since, further, the value of the surface gravity is very 
nearly the same for G and K dwarfs, the relation found will be valid 
for dwarfs. 

The effective temperature of the umbra of medium to large-sized 


. 4/ . . . 
spots is 5713 V0.42 = 4600° as determined from the total radiation,”® 


4450° from the surface brightness,?? and 4700° from the degree of 
ionization.*° We adopt 4600”. 

The spectral type of these spots is usually given as Ko. We as- 
sume that this means that the type is within the limits G8 and Ka. 
With the limits indicated, the value has been plotted in Figure 4. 
An accurate comparison with dwarf K-type spectra would be of 
great interest. 

TEMPERATURES OF EARLY-TYPE STARS 

1. Partly on the basis of theoretical work by Unséld** and Panne- 
koek*? and from Pannekoek’s recent review of the methods of tem- 
perature determination,*3 we derive the following conclusions. 

Effective temperatures depend in a somewhat irregular way on 
stellar colors; the relation is not yet fully accounted for by the the- 
ory of the continuous absorption coefficient. The relation is further 
complicated for the late-type stars by the presence in the spectrum 
of numerous lines and bands and for the early-type stars by space 
reddening. Nevertheless, colors will be useful for giving interpola- 
tory values between certain points in the spectral sequence for which 
effective temperatures may be derived. 

Temperature determinations based on the Boltzmann function 
show certain anomalies which make this method uncertain at the 
present time. 

2. For the O and B stars the best determination may be made 


27 Physik der Sternatmos pharen, p. 363, 1938. 

28 Thid., p. 361. 31 Zs. f. Ap., 8, 32, 1933; P- 225, 1934- 
29 Tbid., p. 363. 32 M.N., 95, 529, 1935. 

3° Ibid., p. 366. 33 Ap. J., 84, 481, 1936. 
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from a study of the maxima of certain spectral series, rediscussed on 
the basis of the theory of the continuous absorption coefficient, as 
Pannekoek has done. His Table IV*4 contains the necessary informa- 
tion, which seems particularly valuable since each element gives an 
independent temperature determination in which possible variations 
of abundance with spectral type enter as accidental errors. Of course, 
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Fic. 1.—The stellar temperature scale for types O-G 


the method is independent of space absorption, which affects most of 
the early-type stars. 

Pannekoek’s values for log g = 4.4 (main-sequence stars) are 
plotted in our Figure 1. The accordance between the different ele- 
ments is very satisfactory. Even if in the future it should be neces- 
sary to revise the theory of the continuous absorption coefficient, it 
seems unlikely that the Variations with temperature of the correc- 
tions could be so large as tO\{ffect appreciably the positions of the 

34 Ibid., p. 496. Double weight was given to the column headed ‘‘E. G. Williams.” 
In addition Dr. Morgan has informed the writer that the maximum of Si* should be 


placed at Bo rather than at Ao if the classification near Ao is based on the strength of 
the K line and He 4471. Cf. Morgan, Pub. Yerkes Obs., 7, Part 3, 1935. 
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maxima. The influence on the observed maxima of a gradual transi- 
tion with increasing temperature of absorption lines to emission 
lines should be considered, but it is probably small because of the 
great differences between the ionization and excitation potentials of 
the series used between types Bo and O8. If the effect is present, the 
temperatures found are somewhat too high. 

Pannekoek’s temperatures for log g = 2.4 (giant B and O stars) 
are systematically 3000° lower. 

In Figure 1 are also plotted four temperature determinations of A 
stars by Pannekoek, based on the intensity of the K line. According 
to Pannekoek, there is a possibility of a small systematic error com- 
mon to all four points. 








TABLE 1 
Star Spec. 1o0oomr 
ie, SEO... B3, B3 7"4+t+ 0:8 
6 Aur : At, Al 37° Ea 
RS CVn Note F4, G8 48 +12 
44 BooB... dG2p 7900«(UE CS 
a GemC dK6+, dK6+ 7 i, an ae 








3. Five eclipsing binaries are known which show two spectra and 
have a measurable parallax. They are collected in Table 1. 

For 44 Bootis B no spectroscopic orbit has been determined as 
yet, and the light-curve is still somewhat provisional, owing to the 
proximity of A. The parallax of RS CVn is still very uncertain. The 
other three systems are discussed below, Castor C on pages 455- 
458. 

Pannekoek has already discussed the system of 6 Aurigae,*> but 
our rediscussion leads to a somewhat lower value of the temperature 
and a somewhat greater uncertainty of the result obtained. 

Stebbins found 6 Aurigae at maximum 0.35 mag. fainter than 
8 Tauri. According to Zinner’s compilation, the magnitude of 
8 Tauri is 1.98+0.04 (m.e.) P.D. and that of 6 Aurigae 2.26 + 0.06 
P.D. From these data, and from Table 2, page 478, we find for 
8 Aurigae at maximum on the IPv scale, 2.10 + 0.04 and a value 


% Ap. J., 84, 506, 1936. 
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slightly brighter than 2.04+0.06, respectively. We shall adopt 2.07 
IPv, which happens to be identical with the HD magnitude used 
by Pannekoek. The trigonometric parallax being 07%037+07004 
(p.e.), the absolute magnitude of the average component is + 0.66+ 
0.24 (p.e.). 

There are two main sources of uncertainty in the derived tempera- 
ture: (1) the unknown value of the limb darkening; (2) the uncer- 
tainty in the absolute magnitude, owing to the parallax. 

Shapley has given two photometric solutions of Stebbins’ data, 
U (uniform) and D (fully darkened). Considering these extremes, 
as well as an intermediate case, we have for 7 = 07037: 


/ 














| | 

| U 3(U+ D) | D 
ie at a ©.147 0.153 0.159 
cos 7 fs 0.218 0.223 0.228 
log R/RO 0.414 0.431 0.448 
‘iP | IO150 g80o 9500 





The last line follows from equation (2), page 430, in connection with 
the table of bolometric corrections (Table 6). 

The intermediate solution will be adopted on the basis of Panne- 
koek’s theoretical work on limb darkening,** which indicates a value 
of u = 0.5 for 10,000° in the visual region of the spectrum (selenium 
cell). A slightly different value of u would not change the result ap- 
preciably. 

The second source of uncertainty is of more importance. With 
the intermediate solution a probable error of + 960° in T, is found, 
or +0.04 in log T,. (Obviously, the uncertainty in log 7, is larger 
than +0.24/10, since for this temperature the visual brightness in- 
creases more slowly with increasing T, than the bolometric bright- 
ness.) 

According to Baker, the spectral type is more advanced than that 
of Sirius (Ao) and is nearly equal to that of a, Geminorum (A1). 

We have, accordingly, for Ar: 


log T, = 3.99 + 0.04 (p.e.). (10) 


This value is plotted in Figure tr. 


3° M.N., 95, 733, 1935- 
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4. The system of pu, Scorpii has recently been discussed by Rud- 
nick and Elvey.*7 At maximum, y, Scorpii was found to be 0.60 
mag. brighter than w,.. The magnitude of the latter star is 3.92 PD 
(Zinner), or 3.68 IPv, so that at maximum 4, is 3.08 IPv. The clus- 
ter parallax is o0%0074 (Kapteyn). The probable error of this paral- 
lax depends largely upon the probable error of the PCG proper mo- 
tion, which is 11 per cent of the proper motion itself in each co-ordi- 
nate. The probable error of the parallax is, therefore, at least 11 per 
cent of its value, and the average absolute magnitude of the com- 
ponents, —1.82 IPv, has a probable error of at least +0.24 mag. 

Miss Maury found the photographic intensity of B to be about 
0.6 or 0.7 of A. Since J;/J, = 0.613 (Rudnick and Elvey), we may 
assume that the components have equal radii. In that case a/a = 
ay;/a = 0.380 and b;/a = 0.336. From these data we find by the 
method used in the next article** that R/a = 0.346, and the adjusted 
value of ai/a = 0.370. With this value of R and Miss Maury’s value 
of a sin 7 = 9,534,000 km, together with 7 = 62°0, we find R, = Ry 
= £3970. 

Adopting the bolometric corrections of Table 6 below, which are 
uncertain by probably not more than o.1 mag., we find 16,000° as the 
average effective temperature of the components. The minimum 
value of the probable error of this result (owing to the uncertainty 
of the parallax alone) is +1900°. In logarithms: 


log T, = 4.20 + 0.05 (minimum p.e.) . (11) 


This value is shown in Figure 1 at B3, the average spectral type ac- 
cording to Miss Maury. 

5. An eclipsing system of still earlier type (Bo+), which will be- 
come of great interest when the spectroscopic orbit and the mass 
ratio have been determined, is SZ Camelopardalis, for which Wes- 
selink published an exceptionally accurate light-curve. The sys- 
tem being a member of the cluster NGC 1502, a fairly accurate 
determination of the absolute magnitude will be possible. With a 
distance modulus of 12.0,*° the absolute IPv magnitude of the pri- 
mary is —4.8, if Am = 1.0 is assumed. 

37 Ap. J., 87, :553, 1938. 

38 Cf. the treatment of V Puppis, p. 500. 39 Dissertation, Leiden, 1938. 


4° Trumpler, Lick Obs. Bull., 14, 156, 1930. See also Ap. J., 86, 187, 1937. 
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Since, in general, the parallax is the main source of error in an ab- 
solute 7, determination, it is of interest to know the effect of a cer- 
tain error in the parallax. In Table 2, for six different values of 7, 
the probable error of T, and of log 7, has been computed under the 
assumption that the probable error of the parallax is 10 per cent of 
its value. This corresponds to +0.22 mag. in the distance modulus, 
which for distant objects is really the quantity desired (space ab- 
sorption). It may be hoped that for galactic clusters the probable 
error may be reduced below this value. The probable error of log 7, 
for very high temperatures follows simply from the Rayleigh-Jeans 
law: 


My =C—2.slogT. (7.—- ©). 


It follows from Table 2 that for the range 10,000° —40,000° (Ao—O8) 
the percentage error in 7, is about the same as that in the parallax. 


TABLE 2 


PROBABLE ERROR, IF P.E. OF 7 IS to PER CENT 











| 

a In Te In log Te 
10,000 : | +. 030° +0.040 
15,000 1650 047 
20,000 | 2150 046 
25,000 | 2000 040 
30,000 . | 3400 049 
40,000 +4500 048 
oe) ; | +0.088 





For much higher temperatures the error in 7, would be double that 
of the parallax, whereas for temperatures of G and K stars the per- 
centage error will be less than half that of the parallax. 

6. Stars earlier than O8 are not covered by the maxima discussed 
by Pannekoek. From the publication by E. G. Williams, which 
was extensively used by Pannekoek, we may tentatively derive one 
more maximum, that for He*. A plot of Williams’ measured inten- 
sities shows a steady increase in strength of He* toward the earlier 
types up to 07, while for the only O06 star the strength is consider- 
ably less again. The maximum is probably near O7. According to 


41 Ap. J., 83, 296, Table VI, 1936. 
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Pannekoek, the temperature corresponding to this maximum is 
44,000° for log g = 4.4. 

Another source of information is Beals’s study of the Wolf-Rayet 
stars.” Two stars are mentioned which have not only a Wolf-Rayet 
spectrum but also an absorption O7 spectrum. These stars, HD 
193793 and HD 211856, have mimimum temperatures (determined 
by Beals using the Zanstra method) of about 59,000°*3 and 65,000°.44 

The Zanstra method as applied to stars having relatively small 
shells (Wolf-Rayet stars, Be stars, etc., contrary to planetary nebu- 
lae), was discussed theoretically by Gordeladse,** who pointed out 
that the observed emissions in such objects are only partly due to 
absorptions from the ground state and that absorptions from higher 
levels contribute appreciably. As Dr. Wurm mentioned to the writ- 
er, this fact will be of importance when the observed emissions are 
strong, which will cause many atoms to fall into intermediate levels. 
As Gordeladse states, the temperatures derived by the strict applica- 
tion of the Zanstra method will in this case lead to too-high tem- 
peratures. The average of the two minimum temperatures found by 
the Zanstra method, 62,000°, may therefore be somewhat too high. 

Another estimate may easily be made, as Dr. Wurm suggested, 
from the circumstance that certain emissions are present (He*, C**, 
and Ctt+ in HD 193793, W VI; and Het, N++, and N**t in HD 
211856, W6 or W7). This would lead to temperatures of the order of 
50,000°-55,000°, found by comparison with Pannekoek’s tempera- 
tures for the absorption lines. More definite results could only be 
obtained by a detailed treatment of the behavior of different atoms 
in diluted radiation. 

Gordeladse’s remarks also apply to the emission B stars. Tem- 
peratures derived for Be stars by the unmodified Zanstra method, 
as given by O. Mohler,*® may be in error. In addition, a selection 
effect enters; only stars having emissions above a certain strength 
are discussed, whatever the cause of these emissions may be. This 
probably accounts for the small scatter in the temperatures found 

? Pub. Dom. Ap. Obs. Victoria, 6, 95, 1934. 

43 Op. cit., Table 8, p. 139. 

4 Op. cit., Table 8 and p. 143. 

4 Zs. f. Ap., 13, 48, 1936. 4 Pub. U. Mich. Obs., 5, 43, 1933- 
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by Mohler (which average 18,000°) and the absence of a notable 
dependence on spectral type (Bo-B8). 

7. From Figure 1, which shows all the temperature data collected 
for the O, B, and A stars, mean values may be read, which are found 
in Table 3. The values for 06 and Os, given in parentheses, are mere 
extrapolations. The bolometric corrections for these temperatures, 
derived in the next section, are added. 








TABLE 3 

Spec | log T; B.( | Spec. | log T, | EC. 
Os. | (4.9) 856 WRB Bice ac | 4.27 —2.01 
O06 | (4.8) 4.9 | Bae yay ohne | 4.23 1.79 
O7 } 4.7 43 || Bs | 4.19 | 1.58 
O8 | 4.60 3.8 | Bo 4.16 L.A2 
Og 4.50 3.24 | B8 | 4.09 1.04 
Bo |} 4.40 2.70 | Ao 4.03 0.72 
Br. | w36 2.48 || A2 3.986 0.52 
Ba. 4.33 —2.22 | AS | 3.93 —o.31 








BOLOMETRIC CORRECTIONS 

1. Before attempting to fix the temperature scale of the later 
spectral classes it is advantageous to consider the bolometric cor- 
rections, B.C. = m,.; — mpy. Three computations of this quantity 
as a function of JT, have been made, all assuming black-body radia- 
tion. Hertzsprung*? made the computation for the lower tempera- 
tures, using the sensitivity curve of the eye; Pike‘* did the same for 
the higher temperatures. Finally Hertzsprung has computed the 
bolometric correction for photovisual instead of visual magnitudes.” 
From his data we derive the bolometric corrections given in Table 4; 
we denote with Am, bolometric corrections computed on the assump- 
tion of black-body radiation. It is seen that these corrections to 
photovisual magnitudes are very similar to those for visual magni- 
tudes, although photovisual magnitudes are more nearly monochro- 
matic. 


47 Zs. f. Wiss. Phot., 4, 49, 1906; Eddington, The Internal Constitution of the Stars, 
Pp. 139. 
48 M.N., 89, 539, 1929. 
49 Zs. f. Wiss. Phot., 30, 173, 1931. 
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It is well known that the stars do not radiate as black bodies. The 
values Am can therefore not be used; they need a correction 
which we denote by 6m. This correction can be evaluated only if 
the whole spectral energy-curve is known. Two approximations are 
available to this ideal case: (1) theoretical spectral energy-curves, 
computed on the basis of a certain assumed composition of stellar 
atmospheres, and (2) radiometric observations of low-temperature 
stars (for high temperatures too much of the energy is radiated in 
the inaccessible ultraviolet). 

2. Spectral energy-curves have been computed by Unsdéld*° and 
Pannekoek.*' Unséld adopted for the stellar atmospheres the same 


TABLE 4 
i 3 : ——— 
T Am | E | Am 

2,000 —4.69 || 6,000°... —0.02 
2,500. 2.84 || 8,000. 0.06 
3,000 nr. 74 || 10,000. 0.28 
3,500 1.07 || 12,000 0.55 
4,000 0.64 |} 15,000 0.97 
4,500 0.36 =|} 20,000 —1.61 
5,000 —o.18 | 

|| 











abundance of hydrogen, one-third by weight, as had been found in 
steliar interiors (B. Strémgren, Eddington). This ratio corresponds 
to 14 hydrogen atoms for each metal atom. Burkhardt, in discussing 
the intensity jump at the head of the Balmer series, later** considered 
50:1 for this ratio, whereas Pannekoek adopted 1000:1. We have 
used Pannekoek’s computations in preference to Unsdéld’s because 
for T = 10,080° Unsdld predicts the intensity jump at 3647 A to be 
0.62 mag., against Pannekoek’s 1.20 mag. The observations® give 
about 1.2 mag. (Although this agreement with Pannekoek’s predic- 
tions does not necessarily confirm the high abundance of hydrogen 
used by Pannekoek, in view of numerical uncertainties™ in the theory 

50 Zs. f. Ap., 8, 241 ff., 1934; Physik der Sternatmos pharen, pp. 131-138, 1938. 

AD. J., 84, 484, 1936. 52 Zs. f. Ap., 13, 56, 1936. 


53 Handb. d. Ap., 7, 460. See also Kienle, Erg. d. exakt. Naturwiss., 16, 462, 19373 
Unsild, Physik der Sternatmos phiren, p. 134, Fig. 530. 


54 Unsdld, of. cit., pp. 135-136. 
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of the continuous absorption coefficient the agreement is important 
for the derivation of the bolometric corrections). 

In spite of this agreement on the strength of the Balmer con- 
tinuum, the corrections 6m obtained from Pannekoek’s data cannot 
be considered as final. It is true that for temperatures higher than 
6500°-7000° the theoretical relation between color temperature (de- 
fined by the slope of the spectral energy-curve) and effective temper- 
ature is in fair agreement with observation for such widely different 
compositions as have been adopted by Uns6éld and Pannekoek; but 
for the solar temperature the difference between the two tempera- 
tures is predicted with the wrong sign by both authors, a difficulty 
which has not been cleared up.** *° 

For three effective temperatures—10,080°, 5040°, and 3150° 
Pannekoek has computed spectral intensities for most of the spec- 
trum; for other temperatures, only at three wave lengths, 4000, 
5000, and 6000 A. A complication arises from the fact that the re- 
lation between surface and effective temperatures, 7! = 2 7}, no 
longer exactly holds. Pannekoek gives the necessary corrections, 
from integrations over the computed spectra, for the three tempera- 
tures stated. In the absence of detailed computations we adopt in- 
terpolatory corrections for the other temperatures. The result of 
these corrections is to increase the total bolometric correction by 
values varying between 0.08 mag. (7, = 5040°) ando.25 mag. (7, = 
25,200, the highest temperature considered by Pannekoek). 

Applying these small corrections to the computed spectral inten- 
sities of Table 1,°”’ we find, for \ = 5500 A, the 6m values given in 
Table 5. The true bolometric correction to photovisual magnitudes 
is, then, B.C. = Am + 6m. 

It so happens that 6m reaches the maximum value of +-o.10 mag. 
for 7, = 6600° when Am is zero. In order to avoid bolometric cor- 
rections with different signs, we subtract ++-o.10 mag. from all the 
6m values before adding them to Am. The bolometric corrections 
thus obtained are given in Table 6 for temperatures up to 50,000”. 


55 Cf. also Minnaert, Zs. f. Ap., 12, 260, 1936. 
6 During a symposium at the Yerkes Observatory (June, 1938) Dr. Wildt an- 
nounced a mechanism of continuous absorption that may remove the difficulty. 


57 Op. cit.. p. 486 
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The correction 6m has been extrapolated for temperatures higher 
than 25,000°, assuming that it slowly decreases. 

For low temperatures, 6m is not well determined by Table s. 
Furthermore, the values given are based on a smoothing process for 


TABLE 5 
| eee = a 
i? | bm 1] Ta im 

- 
3150 —1.26 || 8400 =. 16 
5040 —=6..00 10050 21 
5000. +0.01 1} 12600 29 
6200 +0.09 || 16800 | 46 
7200 +0.00 |} 25200 —0O.44 


TABLE 6 
THEORETICAL VALUES OF BOLOMETRIC CORRECTIONS 


(Temperatures in 1000°) 


r B.( r, | B.< 
2.0 — 32M; 13 — 78 
3.5 re 14 1.35 
4:0. ae te | 2s r.51 
4.5 0.63 |} 16 1.66 
5.0 0.34 || 17 | 1.50 
5-5 O:r7 |) 18 1.94 
0.0 0.00 | 19 2.06 
6. s. 0.00 |} 20 2.18 
| | 
(fae, 0.01 |} 22 2.40 
7% oO. 12 | 25 2.69 
8.0 0.22 |i 30 3.12 
Q.0o. 0.40 i} 35 2.5 
10. 0.57. || 40 | 3.8 
, | 
II 0.78 1} 45 | 4.1 
12 —o.98 | 50 —4.3 


| 
| 
| 
| 


absorption bands and lines that can be only roughly correct; also, 
as we have seen, for the solar temperature an incorrect value of the 
color index is predicted. Fortunately, it is for the low temperatures 
that the second source of information regarding bolometric correc- 
tions is important. 

3. The radiometric observations by Pettit and Nicholson** pro- 
vide extremely valuable information about both the bolometric cor- 


58 4p. J., 68, 279, 1928; Contr. Mt. W. Obs., No. 360. 
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rections and the temperature scale. The theoretical work discussed 
above indicates that the observed energy-curves of low-temperature 
stars are displaced toward the red with respect to the Planck curves, 
that correspond to their effective temperatures. The displacement 
may be about to or 15 per cent. If this result is correct, the 
temperatures derived from water-cell transmissions and heat in- 
dices, being essentially color temperatures, should, of course, be too 
low; but the bolometric magnitudes obtained would still be correct, 
since the reductions to ‘‘no atmosphere,” Am,, were computed with 
a Planck curve corresponding to the observed energy distribution. 
Furthermore, the corrections Am, (which include instrumental 
losses) vary between 0.4 and 0.6 mag. for spectral types between F'5 
and M6, indicating, first, that nearly all the energy reached the in- 
strument (so that deviations from a Planck curve cannot have en- 
tered seriously) and, second, that a small error in the assumed tem- 
perature hardly affects the derived bolometric magnitude. 

From Pettit and Nicholson’s data, bolometric corrections were de- 
rived for all stars of classes Ao-M6. For Fo—-M the results are given 
in Table 7. The radiometric magnitudes, m,, were corrected by 
—Am,; then magnitudes on the IPv scale were obtained®’ from Zin- 
ner’s compilation or Potsdam Publications, Volume 17; for Groom- 
bridge 34 my own measures were used, and for B 4342 the value was 
taken from Mrs. Gaposchkin’s compilation. Variable stars with 
ranges not less than 0.2 mag. are indicated by parentheses. 

The differences, (m, — Am,) — my,, give the bolometric correc- 
tions apart from the zero point. The latter was made equal to that 
of Table 6 by a comparison of the corrections derived by the two 
methods for stars of types dFo-Go. This is the only interval where 
both the theoretical computation and the radiometric determination 
should be approximately correct. The comparison is shown in Table 
8; the temperature scale of Figure 1 was adopted. The difference 
F2 — Go is nearly the same as that found by both methods. There 
is a small discrepancy for the A stars, probably because the devia- 
tions from the Planck curve were accounted for in the first computa- 
tion but not in the second; the difference is in the right direction. 


59 The corrections are given in Table 2, p. 478 
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We shall adopt Table 6 for stars earlier than Fo and Table 7 for 
stars later than Fo. 

With the zero point just defined (which makes the theoretically 
determined bolometric correction zero at 6600°) the constant correc- 


TABLE 7 


Star IPv Sp. | M B:¢. Star IP\y Sp. Mpv B.C. 
6B Cas 2.37 dF2+1.6—0.02 6 And 3.42 gK3 + 0.5 —0.87 
a Tri 3.47 dF4+2.0— .07° 53 Eri AB 3.94 gK1 + 1.3 0.50 
a CMi 0.45 dF3+2.8— .07 BCnc 3.607 gK4 — 0.4 1.21 
a Per 1.90 cF4—4.4— .12. aHya 2.18 gK4 — It 1.12 
§ Peg 4.13 dF6+2.6— .02 ¢ Cep 2.56 cKi — 4+ 6.34 
7 Boo 2.82 Fo+3.1— .16 a Tau t.o2 gKk5 — 0.7 1.46 
e Leo 3.09 Gri-2. — .23. y Dra 2.41 gK5 Gree 1.31 
B CVn 4.26 dGo+4.4— .18 & Cyg 3.82 cK4 —4 1.50 
d Ser 4.40 dGi+4.3— .o9 61 Cyg. 5.43 dK5 +78 0.94 
y Cyg 2.29 cF7—5. — .13 61 Cyg B.| 6.09 dKo + 8.5 1.06 
a Aqu 3.18 cFo—3 - .47 B And 2.22 gKs5s+ + 0.3 1.63 
t Per 4.08 dGo+3.7— .o1 Cir218 6.52 dK6 + 8.2 1.04 
a Aur Sa 0.94 gGi1—o.1— .35 B 4342 7.57+ \*dK4+ [+ 7.3 0.944 
B Com 4.27. dFg+4.9+ .03 a Sco 1.30 cMo — 3.9 2.54 
n Peg 3.09 gG2—1+— .48 Gr. 34 8.09 *dMi+ +10.4 2.08 
n Psc 3.78 gG6 o. |— .36 a Cet 2.81 gMo — I 1.96 
8 Her 2.80 gG6—0.7— .51 a Ori (0.96) cMi+ — 44+ (2.55 
B Cet 2.18 gG6+o0.8— .48 7 Leo 4.78 gM1 ° 1.66 
e Vir.. 2.90 gG6+0.7— .35 B 2935 7.50 *dMe2+ +10.4 1.98 
r Cet 3.53 dG6+5.9— .10 mVir eee gM2 ° 1.98 
o Tau 3.75 gG6—0.4— .36 rAqu 3.83 *gM2 0.5 1.89 
y Tau 3.80 gG8+o.g— .38 8B Peg (2.62) gM3 — 1+) (2.25 
a UMaAB1.96 gGg—0.5— .64 55 Peg 4.60 *oM1 + o+ 1.48 
e Boo A 2.66 Gogo-—2. — .61 15 Tri 5.58 *9M2+ ° 2.06 
8B Aql 3.72 Got3.2-— .18 yw Gem 3.10 gM3 I 2.36 
6 Eri 3.57 dGo+3.8— .32 6 Vir 3.62 gM3 ° 2.31 
8B Gem 1.23 gGo+1.2— .48 30Psc 4.62 *oM 3 + 1: 2.53 
o. Eri A 4.38 dK2+5.9— .15 p Per (3.65) |\*2M4 ° (2.92) 
y UMa...3.12, gK1+0.8— .65 51 Gem.../(5.18) gM4+ — 0+ (2.79) 
a Boo 0.04 gK2—0.3— .84 _ B 3100 5.39 gM4+ ° 2.95 
t Psc 1.55 gK2+1+— .37 13 Sge (5.52) |*2M4 ° (2.68) 
a Ari .2.03 gK2+0.4— .59 62 Sgr 4.57 *oM4 ° 2.83 
e Eri.....3.69 dK2+6.1— .15 56 Leo (5.88) *gM5 — I (2.99) 
70 Oph A..'4.20 dKi+5.7— .05 aHerA. (3.52) Ms -—2 (4.27 
a B..5.90 dK4+7.4— .53. RLyr....(4.3+) *gM4t+ — 1+ (3.7+) 
y And A. .'2.33 gK3—1.4— .95) 45 Ari (6.02) *gM5 o (—4.28) 
@ SORE sw: 2.76, gK2'-+1.0—0.62 





tion to the Mount Wilson values of m, — Am, is +0.62 mag. This 
same zero point has been adopted in Table 7. 

The spectral types in Table 7 are mostly due to Dr. Morgan, who 
has recently obtained a very homogeneous and accurate set of spec- 
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tral types for the brighter stars.® A small systematic deviation for 
the G stars with respect to the Mount Wilson system has been re- 
moved, but for the K and M dwarfs the types are left unchanged. 
The system for the giants is practically identical with that used at 
Mount Wilson. A few determinations, indicated by an asterisk, are 
my own. For a Orionis, Morgan obtained cM3; but since it is possi- 
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Fic. 2.—Empirical bolometric corrections for types F—M. The three curves repre 
sent the dwarfs, giants (M,=0.0), and supergiants (4/7, = — 4.0), respectively. 


ble that the spectral type varies slightly, as well as the magnitude, 
the average, M1-++, of all determinations available has been adopted. 

The absolute magnitudes are based partly on up-to-date trigo- 
nometric determinations and partly on spectroscopic data, some of 
which were kindly supplied by Dr. Morgan. 

Figure 2 shows the relation between bolometric correction and 
spectral type, taking into account the absolute magnitude. Three 
lines are drawn in the diagram, indicating the mean positions for 


6° Mostly unpublished. 
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main-sequence stars, stars with M,, = 0.0 (giants), and stars with 
M,v = —4.0 (supergiants), respectively. The curve for M = 0.0 is 
well determined; the normal points, indicated by crosses and inde- 
pendent from each other, are seen to fit the curve very well indeed. 

The absolute-magnitude effect is very clearly shown, more so than 
in Pettit and Nicholson’s original diagram. This result justified the 
rediscussion of the accurate radiometric observations with other im- 


TABLE 8 











Spec. | No. Theoretical Radiometric 
al | 
NE sa xia entnaicee | 11 | —0.53 —0.42 
dF2. | 4 oI .O4 
dGo 5 | —o.10 —o.06 





TABLE 9 


EMPIRICAL BOLOMETRIC CORRECTIONS FOR DIFFERENT Mpy 

















| | = 
Sp — — M =0.0 | M=— $.0]| Sp | Manin 5s M =o0.0 M =—4.0 
quence | } quence 

| | | - 
F2 —0.04 | —0.04 | —0.04 | K4 | —0.55 | —1.11 | —1.56 
F5 04 08 0.12 || K5 | 0.85 b.45 1.86: 
F8 05 17 0. 28 || K6 | r.¥4 |. a Perse are 
Go | 06 25 0.42 || Mo ee E55 2 2 
G2 07 31 0.52 || M1 1.70 E:72 2.6 
G5 | 10 30 0.65:|| M2 | eg 1.95 3.0: 
G8 | 10 47 0.80 | M3 | (2238) 2.26 | —3.6: 
Ko | 11 54 ©.93:|| M4 | (2.7) i a Paneer 
K2 BG 72 1.20 || Ms 1— (3.1:) | —3.4: 
K3 | O..35. |) —Ors0th —b3541 | 

| | 

! | 








proved data. The absolute-magnitude effect is in agreement with 
that shown by colors; it disappears in the early F stars. 

Table 9 gives the bolometric corrections for the three absolute 
magnitude groups, read from the curves of Figure 2. This table may 
be used for stars for which no radiometric data are available. For 
the M stars the increase in bolometric correction is so rapid, how- 
ever, that a direct determination of the radiometric magnitude 
should be made for all stars where the bolometric magnitude is need- 
ed. This, of course, is also true for the K supergiants, particularly 
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when the absolute magnitude is uncertain. But for earlier types 
Table 9 should be useful if a good spectral type, as well as an approxi- 
mate value of the absolute magnitude, is known. The uncertainty 
in the bolometric correction should then not exceed about o.1 mag., 
which is good enough for most purposes. 

No radiometric observations are available beyond dM2. The 
color index increases continuously for the later dwarfs to 2.0 mag. for 











I 2 3 4 5 6 7 8 9 10 11 Mpy 

BC: T T T T T T T T T 

— 

2.0 }— le — 

a ye aa 

I.5- f od 

LL . A 

r.O;— 3 

e 

0.5 . 4 

h ’ i =! 
ota? wre ! 





0.0 


Fic. 3.—Empirical bolometric corrections for main-sequence stars 


Wolf 359 (M8), and there can be no doubt that the bolometric cor- 
rection will also increase. We shall extrapolate the bolometric cor- 
rection values for types later than dM2; we shall see on page 487 
that the values obtained are not inconsistent with observations of 
colors and spectra. 

Because of the very small scatter on the lower part of the main 
sequence,” the bolometric correction may be plotted directly against 
M,,. Figure 3 shows the result, which is also tabulated in Table to. 
Unfortunately, the relation is not yet very strongly determined. 


% Morgan, Ap. J., 87, 589, 1938; Kuiper, ibid., p. 592. 
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Table 10 will be useful when the spectral type is uncertain. Unre- 
solved binaries will, of course, affect Table 10 but not Table o. 


TEMPERATURES OF K AND M STARS 

1. For low-temperature stars, the spectral energy-curve is far 
better determined by radiometric observations (heat indices and 
water-cell absorptions) than by color indices. The latter define only 
a small part of the short-wave slope of the spectral energy-curve, 
which, moreover, is crowded with absorption lines and bands. 

Yet there are indications that temperatures derived from the ob- 
served spectral energy-curve (uncorrected for absorption lines) may 


TABLE 10 


“MPIRICAL BOLOMETRIC CORRECTIONS FOR MAIN SEQUENCE 














Mpv BC. Mpv B.:€. 
3.0. —0.05 70 —O.5! 
4.0. 05 8.0 0.95 
§.@.% (ols) g.o I.40 
O:67- =6.32 10.0... —EOg 





be systematically too low, by perhaps to or 15 per cent (Pannekoek). 
We could determine this percentage empirically if for some late- 
type stars absolute temperature determinations were available. 

This, fortunately, is the case. For the dwarfs the eclipsing binary 
Castor C may be used; for the giants, the measurements of angular 
diameter. 

2. Castor C.—This important star was discovered by Adams and 
Joy” to be a spectroscopic binary, showing two spectra. Van Gent® 
found that the system is an eclipsing variable with two equal mini- 
ma. The orbital elements were derived by Joy and Sanford,°4 and 
an improved light-curve was subsequently published by van Gent.®s 
Van Gent also rediscussed the dimensions of the system. Later a 
few visual observations made in 1923 were published by Haas.” 

Van Gent in his second paper used 1166 exposures made on twen- 

6 Pub. A.S.P., 32, 158, 1920. 


63 B.A.N., 3, 121, 1926. 6 B.A.N., 6, 97, 1931. 
64 Ap. J., 64, 250, 1926. 0 ALN., 248s £10; 1935- 
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ty-one nights. The probable error of one exposure is +0.046 mag., 
so that the light-curve has considerable weight. The two minima 
have equal depth, which agrees with the equality of the two spectral 
types as observed at Mount Wilson. In his photometric solution 
ran Gent assumed two identical spherical stars. 

The assumption of sphericity is very nearly correct. We may use 
Roche’s model in computing the ellipticity of the components, for 
two reasons: (1) for stars of not too large a mass the concentration 
of matter toward the center will not be very different from Edding- 
ton’s standard model (p./p = 54); (2) the value of r/ais small (0.16). 
Taking the periods of rotation to be equal to the period of revolution, 
we find,” with r/a = 0.16, that the ratio between the two equatorial 
axes is 1.0063. Hence the coefficient of the ellipticity correction is 
0.0068 mag. Only 0.23 of this value, or 0.0016 mag., enters in the 
part of the light-curve near the minimum that has actually been 
used. 

The assumption of equality in size is not quite correct, since the 
Mount Wilson observers found the ratio of the intensities of the 
spectra to be about 5:4. The surface brightnesses being the same, 
the radii are about in the ratio V5 : V4. 

Van Gent obtained r/a and sin 7 from a least-squares solution in 
which the weights of the equations of condition were carefully de- 
termined. We shall examine whether his solution may be retained 
for the average componeht even if the two radii are slightly unequal. 
The results will have also a more general interest.” 

For two circular disks of equal and uniform surface brightness, 
‘ach with radius 7, the fraction of light lost is 


I icemeaiomes 
a (arc cos D — D V; — D*) : (12) 
- ; 
where D is the projected distance between the centers, expressed in 
terms of 2r. D varies, therefore, between D = o (eclipse central) and 
D = 1 (first or last contact). 

67 Cf. p. 496 of the next paper. 

68 Since along the main sequence the radius varies very slowly with absolute magni 
tude, the case of nearly equal radii is of special importance. The conclusions reached 
here explicitly are in agreement with Russell’s statements, Ap. J., 35, 330, 1912. 


rT 
a 
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Suppose, now, that the radii are slightly different but with their 
sum unchanged: 7, = r(1+e) and r, = r(1—e). Then the fraction 
of light lost, neglecting higher than second powers of e, is 








[x — D*}3/2 
) ; (13) 


: (are cos D — D V1 — D2 — ¢ D 


This expression is valid for all values of D except very small ones 
(D?>>é€). Hence, (13) and (12) should be compared only for par- 
tial eclipses. 

The first power of €is lacking in (13), as might have been expected 
from the symmetry of the problem. That means that the solutions 
(12) and (13) differ only in &, ay 2 ~ é. It is not at once 
evident from the equations that this also holds for D values close to 
unity (beginning and end of eclipse). Series development gives 


(13) — (22) 
” * eialeaalaas $e[1 t+a(1—D)+....], (14) 
ifs 
. ‘ II 7V 2 
in which a = —- — _— = 0.06. 
128 384 


For Castor C, r7/r? is about 5/4; hence € = 0.056 and & = 0.003. 
The small correction (14) is approximately equivalent to a small error 
in the grating constant used, amounting to 0.005 mag. per magni- 
tude. This grating constant has a slight uncertainty. Van Gent 
adopts 0.923 mag. and Oosterhoff, more recently,? 0.95 mag. I have 
therefore examined what change is caused in the solution if the larger 
value of the grating constant is used. It appears that r/a is not seri- 
ously changed but that 7 increases by 0° 2 or 0°3. 

It is permissible, therefore, to use van Gent’s solution for the 
mean component, except that the inclination should be increased by 
about o°2. This solution, made for uniform disks, is r/a = 0.16005; 
i = 86°5. 

The amount of darkening, and its dependence on wave length, 
has not yet been studied empirically; this would be particularly 


°°: BAN, F-E8S;-1033- 
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promising in view of the fact that Pannekoek’s theoretical work™ 7 
predicts a strong dependence on wave length. For the average com- 
ponent of Castor C, log g = 4.61 (c.g.s. units); from Pannekoek’s 


Table 26a7 we tind 





ree ‘ | 7000 | 6000 | 5000 4000 | 3500 | 3000 A 
fo: 2.965 £07 Leos 0.33 0.13 | 0.035 
ut 0.73 0.66 0.51 O25 O.115 | 0.034 





Verification of this theoretical result will be most interesting. In the 
meantime we may adopt the theoretical value for 4200 A, u = 0.3. 
A fully darkened solution, kindly made by Dr. A. J. Wesselink, gave 
r/a = 0.169 and i = 85°2. For u = 0.3 we may therefore adopt 
r/a = 0.163 and 7 = 86°%o. 

The apparent magnitude of Castor C at maximum was measured 
by the writer on two nights and found to be 8.93 IPv. The probable 
error is estimated to be +0.05 mag. The parallax being 0%073+3, 
the absolute magnitude of the average component is M/,, = 9.00. 

The spectral type, as determined by Dr. Morgan as well as by the 
writer, is K6+; the star does not show the 770 bands in the visual 
region (< Mo) but is later than 61 Cygni B (K6). According to Ta- 
ble 9, the bolometric correction is —1.28 mag. Table to gives for 
M,, = 9.00, B.C. = —1,40 mag. We adopt —1.32; hence Mj. = 
7.68. The probable error of this quantity is estimated at +0.14 
mag. (of this, +0.05 is due to the apparent magnitude, +0.09 to 
the parallax, and +0.10 to the bolometric correction). 

The luminosity of the average component is, therefore, 3.06+0.14 


mag. fainter than the sun, or log L = —1.22+0.056. 
The value of a sin z found by Joy and Sanford"! is 2,695,000 km. 
With 7 = 86°0 and r/a = 0.163, we have log R/© = —0.198 (R = 


0.633©). From equation (2), page 430, we now obtain 
log T./O = —0.306 + 0.014 + 3(0.198) = —0.207 + 0.014. 


Or 


log T. = 3.550 + 0.014, and 7, = 3550 + 110. (15) 


70 M.N., 95» 733, 1935. 
7 Pub. Amsterdam, No. 4, Tables 26 and 26a (Addendum), 1935. 
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3. Using Michelson’s interferometer method, Pease has obtained 
angular diameters for several stars. A full account of the methods 


Measured diameters 


"2 
72 


and results was published by Pease in 1931. 
were given for 7 stars, obtained with the 20-foot interferometer on 
the roo-inch telescope. The result for the faintest star, a Herculis 
(3.5 mag.) was preliminary?’ and was omitted in a subsequent un- 
published paper by Dr. Pease, which was kindly communicated to 
the writer. The star will also be omitted here. For a Orionis the 
mean of all g determinations listed by Sanford,74 070411+0%0016 
(p.e.), will be adopted (when mirror separations were given only, \ = 
5650 A was used to find the diameters). The data for the 6 stars are 


TABLE 11 

Star d’’(Obs.) IPv Mo] log Te Sp. 
a Boo 07020 0.04 —o.80 3.61 gK2 
a Tau 020 1.02 —oO.44 3.58 gKs5 
@Ori:.... O41 (0.96) —1.59 3.54 cMi+ 
B Peg. O21 (2.62) +0 .37 3.40 gM3 
o Cet (M) 047 var. —o.14+ 3.38 gMo+ 
a Sco 0.040 1.30 —1.24 x. /Sy cMo 




















found in Table 11; for a Orionis the diameter should be better known 
than for any of the other stars. 

A particularly important feature of the 20-foot interferometer is 
the auxiliary apparatus for measuring the influence of seeing on the 
visibility of the fringes; it is probable, therefore, that the diameters 
obtained with it have no serious systematic errors. 

The computation of effective temperatures from diameters, near- 
ly free from assumptions about the spectral energy-curve, is made 
by using the radiometric observations of Pettit and Nicholson. 
These authors have already indicated the method ;7 for our purpose 


we write the formula in the form 
log T. = 2.703 — 3 log d” — {ymbo1 . (16) 


2 Erg. d. exakt. Naturwiss., 10, 84-96, 1931. 

73 Ann. Rept., Mt. W. Obs., 1925. 

74 Contr. Mt. W. Obs., 20, 270; Ap. J., 77, 117, 1933- 

3 Contr. Mt. W. Obs., 16, 480; Ap. J., 68, 300, formula (4), 1928. 
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The constant 2.703 follows at once by inserting the values for the 
sun. 

The quantities d” (obs.) of Table 11 have been derived from the 
interferometer measures by assuming uniform disks for the stars. 
In the original paper by Michelson and Pease” the influence of dark- 
ening toward the limb was examined and the results were tabulated. 
But a recent paper by Diinnweber” considers darkening in the cus- 
tomary analytical form; the effect is tabulated against the ordinary 
darkening coefficient. We shall, therefore, use Diinnweber’s re- 
sults. 

The order of magnitude of the correction factor to d’’ is found as 
follows: for the stars in Table 11, log g = 2 approximately (c.g.s. 
units). For this value of g, T. = 3350° (5040/7. = 1.5), and \ = 
5650 A, Pannekoek’s Table 26a gives C = 1.1 (u = 0.53). With 
this value of C (or 8 in Diinnweber’s notation) the factor with which 
the diameters should be multiplied is 1.07. After preliminary tem- 
peratures were found by using this factor, a second approximation 
was obtained for each star separately; the factor is 1.09 for a Boo; 
1.08 for a Tau; 1.07 for a Ori, a Sco, and 6 Peg., and 1.025 for o Ceti 
(Max.). 

The theoretical question whether the darkening toward the limb 
could be quite different from that computed by Pannekoek, on ac- 
count of the possible presence of very large atmospheres (VV Cephei, 
¢ Aurigae) has not been examined. If the effect would correspond to 
a diminished effective surface gravity, Pannekoek’s tables show that 
then only a small decrease in darkening would result. 

The data on magnitudes and spectra are taken from Table 7. The 
effective temperatures computed by (16), after correcting d” as in- 
dicated above, are given in the fifth column. The probable error of 
the result for a Orionis, which corresponds to the internal probable 
error of d’’ (obs.) found above, is +0.008 in log T, (T. = 3470°+ 
70 ). 

The six temperature values of Table 11 are shown graphically in 
Figure 4. A discussion of the results is given in the next section. 

% Contr. Mt. W. Obs., 10, 86 ff.; Ap. J., 53, 254, 1921. 

7 Zs. J. Ap., 13, 104, 1936. 
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COLOR TEMPERATURES; THE STELLAR TEMPERATURE SCALE 
No attempt will be made here to discuss all the many papers 
in which original measures of color or gradient are given or in which 
such measures are rediscussed. Our present interest is limited to 
finding a set of color temperatures that may be calibrated to effective 


temperatures by means of the results of the preceding sections. The 

















Fic. 4.—The stellar temperature scale for types F-M. Black dots: dwarfs; open 
dots: giants. Circles are absolute determinations; squares are based on reduced color 


data. 


necessity of such a calibration has perhaps not always been apparent, 
in view of the fact that, for instance, the temperature scale given in 
Russell, Dugan, and Stewart’s textbook Astronomy, which is based 
on colors, shows values consistent with other information. But it 
should not be overlooked that the zero point of that scale was ad- 
justed to the effective temperature, not to the color temperature of 
the sun. If the color temperature (6700°+) is used, and if the color 
index adopted for the sun (0.57 mag.) is left unchanged, the color 
temperature of Ao stars is raised from 11,200° to 14,500° (the gradi- 
ent @ = 1.14). This value is in good agreement with the four re- 
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cent determinations for Ao stars made in the same general region of 
the spectrum; three of these are spectrophotometric: 


Greenwich 1932*......¢@) = 1.07+0.05 (p.e) 
Greenwich 1933f........ ©.Q1+0.05 
ere 1.15+0.05 
IVEY: <faetets s's'6 o = 1.04, or T,, = 17,000 
* M.N., 94, 488, 1934. 
t Ibid., p. 505. tErg. d. Exakt. Naturwiss. 16, 446, 1937. 


The fourth is the photoelectric determination by W. Becker.” His 
value is 16,400°. The spectrophotometric values may be expected 
to be somewhat higher than those based on colors, since the absorp- 
tion lines, which crowd toward the violet, are avoided in the spectro- 
photometry. 

That the high color temperatures for Ao stars are, in fact, to be 
expected theoretically was shown by Unséld”? and, more explicitly, 
by Pannekoek.* 

For the classes earlier than M1 we shall base our further discussion 
on two sources of color temperatures: (1) W. Becker’s recent deriva- 
tion,’* which is of considerable weight for two reasons: (a) the 
photoelectric colors are accurate and homogeneous, and (b) a large 
number of stars (738) is included; (2) Brill’s temperature scale,” 
which is based on most of the important series of determinations of 
color temperature published previous to 1932. 

The two scales of color temperature, expressed as c,/7T, are given 
in the second and third columns of Table 13. The spectral types in 
the first column are on the same system as those in Table 7: the 
giants are on the Mount Wilson scale, and also the dwarfs of type 
arlier than K2. For the later dwarfs the scale adopted by Dr. 
Morgan and the writer” is used because of its applicability to very 
faint stars.*? In this system the M stars (giants and dwarfs) are 


8 Veriff. Berlin-Babelsberg, 10, Heft 6, p. 8 (Table 5), 1935. 
99 Zs. f. Ap., 8, 246, Fig. 4, 1934. 

80 M.N., 95, 532, 1935: 

8: Handb. d. Ap., 5, 193 and 199, 1932. 


82 An M star of photographic magnitude 15.0 can be obtained in less than one hour 
with the spectrograph used by the writer on the 40-inch telescope. 
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classified by the strength of the 77O bands in the red and orange.* 
Whereas the classes thus obtained for the giants appear to be in ex- 
cellent systematic agreement with the Mount Wilson classes, the 
dwarfs differ systematically by —o.1 class at Mo. At M4 there is 
agreement again. 

The spectral types used by Becker and Brill seem, on the average, 
to be nearly on the Mount Wilson system. We may therefore, ac- 
cording to the preceding paragraph, adopt their temperatures as 


TABLE 12 


BECKER’S COLOR INDICES FOR G AND K DWARFS 














Star Spec. Che | Star | Spec. Cpe 
n Cas Go 0.105 61 UMa | G4 0.247 
pw Cas G6 183 || 8 CVn | Go 115 
HR 483 Gi 143 || 43 Com | Fo Ey 
107 Psc Ko 304 | 59 Vir | Go+ (.158) 
HR 511 Ko 300 «|| & Boo | G8+ . 226 
6 Tri G1 124 | n CrB | Go .120 
HR 753 K3 381 || d Ser | Gr .156 
t Per Go 130 || 26 Dra | Gr .105 
x Cet G5 197 || 70 Oph | Kir+ 325 
39 Tau G2 124 || HR 7162 | Go o81 
Aur Go 157 || o Dra G8+ 254 
x' Ori... . Go 115 ! 61 Cyg K5(+) 541 
vs Aur Go 074 || 51 Peg G4 (.179) 
11 LMi | G8+ 276 || HR 8832 K3 402 
15 LMi Go+ 112 || 85 Peg G2 0.160 
47 UMa | G2 0.162 | 
| 





they stand, except for dwarfs later than Ko. Since Becker gives the 
color indices for the individual stars, we shall make a new reduction 
for the G and K dwarfs, using spectral types in the system adopted 
here. The data are collected in Table 12; as we shall see later, they 
are of great importance in a comparison to be made with the sun. 
For binaries observed as one star (70 Oph., 61 Cyg. etc.) the weight- 
ed mean spectral type of the components is given. The writer is 
again indebted to Dr. Morgan for the spectral types used here. Par- 
ticular care was taken to omit all subgiants from Table 12, since the 
colors of these objects are decidedly redder than those of ordinary 
dwarf stars. 

83 Theoretically, these bands have the advantage over those in the blue because they 
are more fundamental. 








464 G. P. KUIPER 


A plot was made of the data in Table 12, and normal points were 
read off at Go, G2, etc. The plot shows the small scatter to which 
Morgan has drawn attention.*4 For this reason the averages should 
be of value in spite of the limited material on which they are based. 

From the mean color indices, temperatures in Becker’s system 
were found by means of Table 4 in Becker’s paper. The results are 


TABLE 13 
THE STELLAR TEMPERATURE SCALE (LAST COLUMN)* 


oi 























| | 
%. | «2/Te | G/T: | Becker | Brill | 4 0 | wi 1 ah | afte |» 
(Becker)} | (Brill) | Corr. | ‘Corr. | (Prov Adop 
ay | | | | | 

Ao: .. } 0 0.88 | 1.060 | 1.44 32 E43 1.33 | 1.34 [10,700 

AES «i ie isep | Asse | 1.83 | CeS3 T.68 | 1.68 | 1.68 | 8,530 

Fo ) SELAS | 1.68 | 1.99 | 1.84 I.g2 | I.Q2 | r.02 | 7,800 

Fs | 1.67 | 2.05 | 2.23 | 2.21 2.22 | 2.22] 2.22 | 6,470 
dGo |; ameSa | 2024 | 2.40 2.39 2.40 | 2.40 | 2 40 | 6,000 
dG2 | 1.95 | (2 35) 2.51 2.51 2.87 4 2.51 2.81% | 5,710 
dGs.. | 2.10 | 2-53 | 2.06] 2.69 | 2.68 2.68 | 2.68 | 5,360 
dKo....| 2.42 | 2 64 2.98 2.80 2.92 | 2.92 2.92 | 4,910 
dK2.. 2.600 ; 3.10 ‘ 3.20 ce 2.30 3.08 | 4,650 
dK5+ | <e98- 1 3.84 3.84 | 4.28 3.92 3.69 | 3,900 
dM2 | | 5.16 4.87 4.49 | 3,200 
gGo... | 2.33 | 2.48 | 2 87 2.64 2.76 | er 2.76 2.76 5,200 
gG5 ie 2.95 3.08 ie 3.10 3.46 3.14 3.10 | 4,620 
phO-.:,| $2490] (3.28 |) 63.55 | 3-44 3.50 3.83 3.52 3.30" | 4,230 
gKs5....| 3.83 | 4.35 | 4.30) 4-38 | 4-85 | 4.58 | 4.32] 4.00] 3,580 
gMo | 4.04 | 4.43 | 4.60 | 4.59] 4.60] 4.81 4.56 4.22 | 3,400 

zi | 

gM2... Js | | | | 5.16 4.87 4.49 | 3,200 
gM4....| | 5.63 | 5.34] 4.90 | 2,930 
gMo6 | | 6.00 ay 5.22 | 2-750 
gMB8e.. .} | 6.38 | 6.09 5.54 | 2,590 


| 


* Interpolations should preferably be made between the recipocal temperatures (cz = 14,320) of the 
preceding column. The spectral types are on the Mount Wilson system, except for the K and M dwarfs 
(cf. pp. 462-463). 


entered in the second column of Table 13 of this paper (dGo-dK5+) 

The data in the third column are based on Brill’s investigation. 
Brill derived first the temperature scale for the giants and then aver- 
aged the differences (giants-dwarfs), expressed in reciprocal tempera- 
ture, taken from 5 determinations.*' Fearing a systematic error with 
respect to the other determinations, we have omitted the second 


84 Ap. J., 87, 460, 1938. 
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source quoted by Brill and have recomputed the differences between 
giants and dwarfs. 

The transformation of color temperatures to effective tempera- 
tures may, to a first approximation, be made by adding a constant 
to all reciprocal temperatures, as is suggested by our remarks on the 
temperature formula in Russell, Dugan, and Stewart’s Astronomy. 
It remains to be seen whether this simple reduction, which involves 
only one constant for the whole temperature scale, is sufficient or 
whether additional corrections are required. 

We reduce the two temperature scales to the sun by adding 0.56 
and 0.16, respectively, to the reciprocal temperatures; this leads to 
c,/T = 2.51 (or T = 5713) for dGz2.*5 The corrected scales are given 
in the fourth and fifth columns, respectively. Except for the A stars, 
where the interagreement is only fair, the two corrected scales agree 
quite well. (The accuracy is most readily judged by remembering 
that the scale of reciprocal temperatures is‘almost exactly double 
that of ordinary color indices.) The average scale is given in the 
sixth column. 

3efore discussing possible further corrections to the temperatures 
in the sixth column, it is preferable to consider the data for spectral 
class M, which is not properly covered by the color data. (Becker 
found even a decrease in color index for the later M giants, undoubt- 
edly due to the increasing strength of the 77O absorption in the 
visual region.) The TiO absorption also vitiates temperature deter- 
minations from heat indices for the later M stars, as was pointed out 
by Pettit and Nicholson.*° Least of all affected are the water-cell 
absorptions obtained by Pettit and Nicholson, which are particular- 
ly significant since they analyze the spectral energy-curves of the 
late-type stars near the maximum. For temperatures much higher 
than 4000°, however, the water-cell absorptions lose their value, as 
is shown by Table IV of Pettit and Nicholson’s paper. 

Accordingly, the color temperatures derived from water-cell ab- 
sorptions, not exceeding 4300°, are added in Table 13, seventh col- 
umn. The data for the giants are simply taken from Pettit and Nich- 
olson’s Table V, whereas the dwarfs (K—M) are again listed separate- 
ly (Table 14). The weighted mean spectral class for 70 Ophiuchi is 


85 Cf. pp. 431 and 438. % Ap. J., 68, 295, 1928. 
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slightly later for radiometric than for photoelectric observations, of 
course. The number of observations made by Pettit and Nicholson 
is given in the last column. 

Three groups are indicated, with mean spectral types dK2, dK5+, 
and dM2. In the first group equal weight will be given to all three 
stars; in the second group the components of 61 Cygni should obvi- 
ously have more weight than the other stars, not only because they 
have been more frequently observed but also because they are bright- 
er. We shall give them three times the weight of the other stars of 


TABLE 14 

Star Spec. W.C. n 

OpPTUAS.6 055s are dK2 oMso 6 
euleri:..:. K2 0.58 = 

70 Oph AB K2 0.63 5 
6r CygA...... ee .5 0.70 7 
61 Cyg B. : K6 0.82 7 
Ci 1278. K6 1.06 2 
ier) ee K4+ 1.03 2 
B a035..... Es M2 1.04 2 














that group. The last star is as faint as 5.4 mag. radiometrically, and 
the water-cell absorption is probably only approximate. We shall 
follow Pettit and Nichelson, who adopt W.C. = 1.14 for dM2,”7 
probably on the basis that the heat indices for dM2 and gM2 are 
equal and that W.C. = 1.14 for gM2. 

With the three average water-cell absorptions, color temperatures 
in Pettit and Nicholson’s system are found from their Table IV (or 
V); the results have been entered in Table 13. 

The six lines in common between columns 6 and 7 of Table 13 in- 
dicate a systematic correction of —o.29 to column 7 if it is to be re- 
duced to column 6. With this correction applied, column 8 is ob- 
tained; equal weights are given to columns 6 and 7. 

The temperature scale of the eighth column has, therefore, the 
following properties: (1) it embodies most of the essential data on 
color temperatures existing; (2) it is homogeneous with respect to 


87 Op. cit., Table 5. 
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the zero point of the color equivalent (this has been obtained by add- 
ing a constant to the three basic systems); and (3) the constants 
make the temperature for dG2 equal to the effective temperature 
of the sun. 

We shall now examine whether this provisional system needs fur- 
ther corrections. The test is obviously twofold: (1) the representa- 
tion of the A stars; (2) the representation of the K and M stars for 
which effective temperatures were derived. 























TABLE 15 

Star Spec. Mpy ¢2/T (Prov.) | ¢2/Te (Obs.) | Obs.—Prov. 
a Gem C dKk6+ +9.0 4.19 4.04 —0.15 
a Boo gK2 —0.3 3.84 3.50 —0.34 
a Tau.. | gks 0.7 4.32 3-79 a 
a Sco.. | cMo 4:0 4.50 4.43 — 28g 
a Ori. . | ¢Mri+ 4+ 4.79 4.14 — .65 
B Peg | gM3 1.0 5.10 4.65 — .45 
o Cet (M). | gM6+ —2: Ceo es 6.06 +0.35+ 

| 





It appears that the A stars are very well represented. The pro- 
visional temperatures are plotted in Figure 1 as crosses; no essential 
improvement can be made. 

The comparison with the late-type stars with known diameters is 
given in Table 15. The fourth column contains the “provisional”’ 
reciprocal temperatures from Table 13, and the fifth the reciprocal 
effective temperatures found on page 458 and in Table 11. The prob- 
able error of the observed reciprocal temperature of Castor C is 
$0.13. 

Table 15 shows that for the lower temperatures our provisional 
scale needs correction. The correction given by Castor C is —o.15+ 
0.13 (p.e.) at c,/T = 4.2; the correction given by the other stars is 
—o0.40+0.08 (p.e.) at c,/T = 4.66. In forming this average, double 
weight was given to a Orionis and half-weight to o Ceti (max.), for 
which the data are, of course, somewhat indefinite. The weight for 
a Orionis would have been increased but for the possible variability 
of the quantities involved. The probable error +0.08 follows from 
the internal agreement of the six objects, with the weights assigned 
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as indicated. If o Ceti is omitted from the average, we obtain the 
correction —0.46+0.06 (p.e.) at ¢./T = 4.57. The absolute amount 
of correction derived from the giants may be somewhat too small, 
because two supergiants are present for which the temperatures may 
be expected to be somewhat lower than for ordinary giants of the 
same spectral type. 

We shall adopt the correction —o.30 for c,/T = 4.2, and —o.40 
for c./T = 5.0. These corrections leave a residual of +0.15+0.13 
for Castor C, and —o0.04+0.08 or —o.11+0.06 for the giants, de- 
pending on which average is used. Since a correction of —o0.29 
had already been applied to the system of reciprocal color tempera- 
tures derived from water-cell absorptions, we find that the total 
correction, both for c,/T = 4.2 and 5.0 (T in the provisional scale), 
amounts to 15 per cent of the original temperature. This solves the 
problem given on page 455: the color temperatures derived from 
water-cell absorptions have to be increased by about 15 per cent in 
order to represent effective temperatures. This result confirms 
Pannekoek’s theoretical conclusions for low-temperature stars.*® 

The data available are, of course, insufficient for showing whether 
this percentage, 15, is constant for the late K and M stars. It is 
reasonable to assume that this is the case to a first approximation. 
For the provisional reciprocal temperatures of Table 13 smaller than 
4.2, we shall adopt smaller corrections and make them vanish at 
c./T = 3.0. For the spéctral types later than Mo, for which the 
temperatures depend on water-cell absorptions only, we shall adopt 
the original color temperatures published by Pettit and Nicholson, 
increased by 15 per cent. This discussion completes the formation of 
the last two columns of Table 13. 

The temperature scale of Table 13 is shown graphically in Figure 
4. The lines representing the giants and the dwarfs show a striking 
similarity to the corresponding lines in Figure 2; this is gratifying 
since the observational data used in the two diagrams are, with one 
exception (dM2), independent. Even details such as the bend in the 
dwarf curve at K2 and in the giant curve at Mo are common to the 
two figures. 

One point of importance has not been fully analyzed: the mean 


88 Ap. J., 84, 490 ff., 1936. 
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absolute magnitude of the giants for which Table 13 gives the tem- 
perature scale. Becker states that in the formation of the averages 
for giants the c stars and red supergiants were excluded. The evalua- 
tion of the precise mean absolute magnitude would require an anal- 
ysis similar to that given for the bolometric corrections (Table 7). 
But the advantage would not be great; it seems far better to use the 
original color measures catalogued by Becker for spectral types A-K 
(provided space reddening may be neglected) and the original meas- 
ures of water-cell absorption by Pettit and Nicholson for tempera- 
tures lower than 4000°. From these measures, effective temperatures 
in the system of Table 13 may be derived. This automatically takes 
account of the absolute-magnitude effect (which is large, and diffi- 
cult to evaluate accurately) and avoids uncertainties of spectral 
type and absolute magnitude. 

The process of obtaining effective temperatures as outlined in the 
preceding paragraph follows immediately from the way in which 
Table 13 was constructed. 

1. If the spectral type is between Ao and Mo, and the star is 
contained in Becker’s Catalogue,*® take the value of the column 
headed ‘‘F.I.” With this value enter a table by Becker for con- 
version into color temperature (the value then obtained is equivalent 
to the second column in Table 13). Add 0.56 to c,/T (fourth or 
eighth column of Table 13); if the result is smaller than 3.0, the 
reciprocal effective temperature has already been found. If the re- 
sult is larger than 3.0, subtract a correction that is 0.00 at ¢,/T = 
3.0, 0.30 at 4.2, and 0.40 at c,/T = 5.0. The final value is the recipro- 
cal effective temperature. 

2. If the star has been measured radiometrically by Pettit and 
Nicholson,” and is of type later than gG5 or dK2, take the water- 
cell absorption from the ninth column,” obtain the color temperature 
from Table IV (or V),” and increase this value by 15 per cent if the 
color temperature is less than 3200°; increase it by 480° between 
3200° and 3600°; add 450° at 4000”, and 420° at 4400°. 

89 Veréff. Berlin Babelsberg, 10, Heft 3, pp. 25-41, 1933. 

9 Jbid., Heft 6, p. 8, Table 4, 1935. 

* Ap. J., 68, Table III (pp. 286-291), 1928; Contr. Mt. W., No. 360. 

92 Op. cit. 
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We shall give two applications of the procedure just described. 
One relates to the giants of Table 15; the temperatures computed di- 
rectly from the color data are found in Table 16. The differences 
between these temperatures and those computed from the measured 
diameters are more irregular than those shown in Table 15, but the 
values in Table 16 should be preferred, since the absolute-magnitude 
effect is now taken into account and uncertainties in spectral class 
are eliminated. If weights are assigned as before (a Orionis double 
weight, o Ceti half-weight), the average of the last column is —o.16 
+0.13 (p.e.). The absolute amount of this deviation is now approxi- 
mately identical to that for Castor C, +0.15+0.13 (p.e.). 




















TABLE 16 

Star €2/ Te (Becker)| ¢2/ Te (W.C.) Av. ¢2/T, tit Diam.—Av 
a Boo 2.62 3.67 3.64 3.50 —o.14 
a Tau 4.2 4.00 4.15 3.79 —o. 36 
a Sco. eee eee 4.67 4.67 4.43 —0.24 
a Ori... 4.62 4.62 4.14 —o.48 
B Peg.. 4.59 4.59 4.65 +0 .06 
o Cet.. 4.90+ 4.90+ 6.06 +1.16+ 








The second application is the derivation of the effective tempera- 
ture of « Aurigae. From Becker’s color index we find 7, = 6870°. 
This figure is somewhat higher than the value 6300° adopted in the 
discussion of this system,%* which was made before the present paper 
was completed. The temperature difference corresponds to 0.37 
mag. bolometric if the dimensions of the star are left unchanged. 
Hence, if the F2 star lies 0.37 mag. above the average empirical 
mass-luminosity relation (which is not very much less probable than 
that the deviation is 0.0 mag.), our dimensions and masses may be 
retained, but the computed luminosities of both components must 
be increased by 0.37 mag. or by 0.15 in log L. The higher tempera- 
ture now found diminishes by a considerable factor the discrepancy 
found by Strémgren® for the ultraviolet intensity. The factor is 


9% Ap. J., 86, 572 and 612, 1937; also 87, 213, 1938. 
94 [bid., 86, 604, 1937. 
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6.6 at 1000 A and 48 at 500A. If the star is slightly reddened be- 
cause of its distance (about tooo parsecs), the factor will be larger. 

In conclusion, some desiderata may be mentioned in connection 
with the results of this paper: 

1. An accurate photometric study of Castor C and 8 Aurigae in 
two or three very different wave lengths. 

2. The discovery and study (photometric and spectrographic) of 
eclipsing binaries showing two spectra, in galactic clusters. 

3. Further determinations of trigonometric parallax for the stars 
in Table 1. 

4. Radiometric observations of Barnard’s star (dM6) for settling 
the question of the temperature scale and of the bolometric correc- 
tions for the late M dwarfs. 

YERKES OBSERVATORY 
June 1938 








THE EMPIRICAL MASS-LUMINOSITY RELATION 
G. P. KUIPER 


ABSTRACT 

After a brief historical introduction the problem is subdivided into two main 
parts. The first, on the temperature scale and bolometric corrections, is treated in the 
preceding article. The second problem, the derivation of the empirical mass-luminosity 
relation, is treated in three sections in this article: (1) the visual binaries, (2) some 
selected spectroscopic binaries, and (3) Trumpler’s massive stars in clusters. An at- 
tempt has been made to obtain the most accurate observational data for all quantities 
entering into the discussion, including magnitudes and spectral types. 

Tables 1 and 7 give the visual binaries used at present. Table 5 shows the derived 
quantities for the stars of Table 1, and Table 6 shows the quantities of theoretical 
interest for the stars of Table 5 for which the accuracy is sufficiently great. Visual 
binaries in the same class as those of Table 1, but for which the data are still incomplete, 
are collected in Table 8. 

The problem of the spectroscopic binaries is only partly treated in this paper. Only 
some representative objects are discussed; they are found in Table 12. Theoretical 
values of the ellipticity and reflection constants are used in the discussion of the ob- 
servations. Trumpler’s massive stars are discussed in Table 13. 

In all three sections the results of the preceding article have been used. The data 
are shown graphically in Figures 1 and 2. 

INTRODUCTION 

The discovery of the mass-luminosity relation has been of great 
importance to the progress of astronomy. The relation has been used 
in statistical astronomy and in double-star astronomy and has been 
a central problem of theoretical astrophysics. Since for most stars 
no direct determination of the mass can be made, the use of the 
mass-luminosity relation is the only method of estimating the total 
mass of the known stars per volume of space—an important dy- 
namical quantity." In binary-star statistics the observable Am 
and a rough knowledge of the absolute magnitude can now be used 
in obtaining a statistically useful determination of the mass ratio, 
which is of great cosmogenetic interest.? Investigations, such as a 
recent one on e Aurigae? show the need of having the relation well 
established. But particularly the theoretical importance, both in ab- 


1 E.g., Oort, B.A.N., 6, 285, Table 34, 1932. 
2Pub. A.S.P.,.47, 17 fi., 143 ff., 1935. 
3 Ap. J., 86, 574, 1937. 
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stract form and in numerical form for particular stars, is clear from 
Eddington’s work and from more recent developments.* § 

1. Historical.—The idea that the mass and the luminosity are cor- 
related seems to have developed gradually as the knowledge of paral- 
laxes increased. Even today the knowledge of masses is to a large 
extent limited by the knowledge of accurate parallaxes. The fact 
that one of the best-known stars, Sirius B, does not follow the general 
mass-luminosity relation must have delayed the discovery of that 
relation. 

Halm was probably the first to state explicitly the existence of a 
statistical relation between intrinsic brightness and mass.° His con- 
clusion was essentially based on the correlation of mass with spectral 
type and of spectral type with luminosity. The first relation was 
established mainly from double-line spectroscopic binaries. The re- 
sult was therefore partly accidental, because double-line spectro- 
scopic binaries are very rare among giants, so that Halm was es- 
sentially dealing with main-sequence stars only. In fact, the relation 
was not recognized in Russell’s paper of 1913,’ although in 1914 Rus- 
sell found evidence® for a definite correlation, which was obtained 
by comparing absolute magnitudes derived from hypothetical paral- 
laxes with absolute magnitudes derived by methods of stellar sta- 
tistics. 

Hertzsprung, in 1918, gave the relation? 


log n= —o.06(M, — 5) ; 


in fair agreement with modern data. He also gave the corresponding 
formula for the mass ratio of a binary derived from Am, and the 
first formula for dynamical (not hypothetical) parallaxes. Van Maa- 
nen shortly afterward emphasized that the mass-luminosity relation 
is independent of spectral type."° 

4B. Strémgren, Erg. d. exact. Naturwiss., 16, 497, 513, 1937. 

5 §. Chandrasekhar, Introduction to the Study of Stellar Structures, ‘Astrophysical 
Monograph,” chaps. vii and viii, 1938. (In press.) 

6 M.N., 71, 638, 1911. 7 Observatory, 36, 327, 1913. 8 Pop. Astr., 23, 340, 1914. 

9 4.N., 208, 96, 1919. Already in 1915 Hertzsprung had given the same relation in a 
less explicit form (A p. J., 42, 115, n. 2). 


10: Pub, ASP. 3312335 1910: 
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Many other investigations followed, most of which have been re- 
viewed by Lundmark." Of all these, Eddington’s results are un- 
doubtedly the best known, particularly because Eddington used the 
relation in connection with his well-known theoretical investigations. 

The most recent study of the subject is that by P. P. Parenago,” 
which was received after most of the present investigation was com- 
pleted. In some ways the two investigations run parallel: for in- 
stance, in the use of bolometric corrections derived from radiomet- 
ric observations and the estimation of weights of the individual mass 
determinations. But the present study uses several improved orbits, 
mass ratios, and Am values and does not use spectroscopic parallaxes. 
Minor discrepancies will be found by comparing Parenago’s Table 2 
with our Tables 1 and 5. The treatment of the spectroscopic binaries 
is entirely different in the two papers: Parenago uses spectroscopic 
absolute magnitudes whereas we use the stellar temperature scale. 
A number of objects are used here which are not considered by Pa- 
renago; among these are Trumpler’s stars. 

2. The present investigation is limited to those individual stars for 
which the three fundamental quantities—the mass m, the radius R, 
and the luminosity Z—may be obtained from the observations with 
fair precision. The derivation of the radius, in addition to the other 
two quantities, does not require any extra information, since radius 
and luminosity are simply related by the effective temperature, 
which is needed in any case before the eclipsing binaries can be used. 
In view of the fact that all three parameters are required in theoreti- 
cal work, we shall give them explicitly as the final result of the in- 
vestigation. 

Since all the emphasis is laid on accurate data for individual stars 
(as required in theoretical work), no attempt is made here to supple- 
ment the results by a statistical treatment of other data of smaller 
individual accuracy. Such a statistical treatment would be useful 
for determining the trend of the mass-luminosity relation for the 
giants and the O stars, where few data of high individual accuracy 
are available. 

A close co-operation between theory and observation will certainly 
1" Handbuch der Astrophysik, 5, Part II, 1933. 
12 A.J. Sov. Union, 14, 33, 1937- 
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be the quickest and least wasteful method of solving the fundamen- 
tal problems of stellar structure and evolution.’ In this connection 
it would seem that the visual binaries of Tables 1, 7, and 8 deserve 
preference in double-star observations and in determinations of 
parallax and mass ratio. At least four to six determinations for each 
object will be needed in order to make the mass determinations suffi- 
ciently accurate. Photometric and spectroscopic measures of double- 
line spectroscopic binaries that are also eclipsing variables are ob- 
viously the second group of observations especially needed. Finally, 
investigations on the stellar temperature scale are required for the 
computation of luminosities from radii and vice versa. 

Many important direct or indirect contributions to the knowledge 
of stellar masses have recently become available: Schlesinger’s Cata- 
logue of Parallaxes; Boss’s General Catalogue of 33342 Stars, contain- 
ing many mass ratios for visual binaries, based on all the available 
meridian positions; many new orbits, of which we want to mention 
especially those derived with the help of photographic measures 
(Strand, van den Bos, Hertzsprung); and several new mass ratios by 
van de Kamp. For the discussion of the eclipsing binaries, of special 
importance is Pannekoek’s rediscussion of the temperatures of O 
and B stars, derived from maxima of spectral series. Data for very 
interesting systems, such as ¢ Aurigae and VV Cephei, have recently 
been published. Last, but not least, Trumpler’s discovery of very 
massive stars in galactic clusters should be mentioned. These many 
advances make it promising to assemble and discuss the data now 
available. 

We have already referred to the stellar temperature scale as an 
integral part of the subject under discussion, being needed in the 
conversion of radii into luminosities and vice versa, according to the 


= (YE) 0 
Lo \Ro)  \T.o)’ ; 


Before this formula can be used, the relation between M,;, and L has 

to be known. This relation is given by the bolometric correction. 
Hence, the subject falls into two parts: (1) the discussion of the 

stellar JT, scale and the derivation of the bolometric corrections; and 


formula 
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(2) the discussion of the binary systems, with suitable observational 
data, and of Trumpler’s stars. The first problem is treated in the 
preceding article; the second naturally falls into three parts: the 
visual binaries, the eclipsing binaries, and Trumpler’s massive stars. 


THE VISUAL BINARIES 
If the orbital elements and the parallax of a visual binary are 
known, Kepler’s third law 


3 log a’ — 2 log P — 3 log x’ = log =m (2) 


will give the total mass, 2m, of the system. If the orbit is well de- 
termined, and hence the error in (3 log a’’ — 2 log P) is small, the 
error in 2m will depend mainly on the error in the parallax; we 
then have 3A(log 7”) = —A(log 2m). If the probable error of the 
parallax equals one-tenth of the parallax itself, the probable error 
in log 2m will be +0.13. This value corresponds to about +1.4 ab- 
solute magnitude in the region of the mass-luminosity relation cov- 
ered by the visual binaries of our list. The computed absolute mag- 
nitude is also affected by an error in the parallax, and a certain com- 
pensation of errors occurs when the data are used in a diagram corre- 
lating mass and luminosity; the value +1.4 mag. is then reduced to 
+1.2 mag. As true deviations from the mean mass-luminosity curve 
larger than 2 mag. are very rare indeed, it is obvious that data with 
probable errors larger than +1.2 mag. are of no individual value. 
Hence, all systems for which the probable error of the parallax ex- 
ceeds 10 per cent are excluded from our discussion. 

The remaining binaries are collected in Table 1. It would be of 
great importance for the knowledge of stellar masses if parallax ob- 
servers would pay special attention to these objects, particularly to 
those for which the probable error of the parallax still exceeds 
+0%004. Of not much less importance are the stars in Table 8, which 
may be used in the near future. 

The first three columns of Table 1 need no explanation. The 
fourth column gives, on the international photovisual system, the 

"3 We have, furthermore, excluded stars with parallaxes smaller than 0’045. No at- 
tempt is made here to use averages for binaries of different spectral types—a problem 
considered recently by Professor Russell. 
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apparent magnitude of the combined light of the binary. Most of 
these data are based on Zinner’s catalogue (which is on the Potsdam 
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TABLE 1 


VISUAL BINARIES FOR THE MAsSs-LUMINOSITY RELATION 











Name | a 1goo | 6 1900 
159 | oh32m2|—25° 10! 
Cas | © 43.0/+57 17 
Eri 1 36.0/\—s6 42 
Eri B¢ 4 10 7]— 7 49 
Aur } 5 og 3] +45 54 
CMa 6 40 7) —16 35 
Gem 7 28.21|+32 06 
CMi 7 34.11+ § 29 
Pup 7 47 1|—13 38 
UMa II 12.9|+32 06 
Vir 12 36.6/— 0 54 
UMaA 13 19 g|-+ss 27 
Cen I4 32 8) - 60 25 
Boo 14 46.3/-+19 31 
Boo .. 15 00.5|+48 03 
10075 16 24.5|+18 37 
Her 16 37 5 | +31 47 
352 16 50.1|— 8 og 
6416 |17 11.5|—46 32 
6426 I7 32.2134: 53 
6510 17 25.2/— 9° 59 
Dra 117 34.0|+61 57 
Her BC 117 42.5|+27 47 
Oph j18 00.4/-+ 2 31 
Her 118 03.2/+30 33 
7162 18 53.3|/+32 40 
Equ 21 09.6/+ Q 306 
Cyg |21 10 8) +37 37 
60 122 24.4;/1T57 12 
Peg |23 50 9| +26 33 














| 
| sp | Spec | Am | F | a” IRef ,0| x” (Trig.) 

5.61 G7 ©.3 | 25.00 | 0%670 | 2,A loto72+ 6 
3.571 Go 5 3-74 |§26.0 |12.534 | 3,B | .182+5 
5.21 Gs 5.06 |251 8.31 4,C 161+7 
9.37| Bo, Mse 1.48 |247.9 6.894 | 1,AB}] .202+3 
0.19} G4, F4 0.15 |104.02d} 0.0536] 1, A 0632 
1.52] A5 10.06 | 49.94 | 7.62 1,A 376+ 2.6 
1.59] Ao, Ar | 0.86 |340 5.84 1,B 07343 
0.45] F3, 10.3 40.2 4.26 1,A 2901+4 
5.27 G | 0.7 23.34 | 0.69 1,A o61+4 
3.77| Fo, Fo | 0.43 | 59.86 | 2.5355] 1,4 138+6 
2.98} Fo, Fo | 0.03 [171.37 | 3.746 | 3,4 o89+7 
2.43} A2, A2 lo 20.54d| 0.0115] 1,A 045 

—o.20} Gy, Ki | 1.37 | 80.09 |17.593 | 5,4 -750+7 
4.51] G8, Ks _ | 1.96 |149.95 | 4.884 | 3, A 147+ 6 
4.79] G1, G5 | 82) |219.5 3.609 | 3,B .079+5 
6.82 | 0.08 [287.1 2.198 | 1,AB} .o5s9+6 
2.93] G1, | 2.80 | 34.42 | 1.349 1,A -II0ot+5 
9.20) M4 | O.r | 1.72 | 0.19 6, AB! -148+4 
5.54| Ko, | 2.6 |242 4.04 7, BC| .132+6 
5.82] Ks, | 1.2 42.2 1.83 1,A .147+6 
5.22] G6 | 0.14 | 46.0 1.06 1,A 050+ 4 
5.22|) Gr, | 2.61 | 76.06 | 1.55 8,B 066+ 5 
9 80} M4 | 0.1 43.02 | 1.287 | 1,A -109+6 
3.97) Ki, Ky | 1.68 | 87.85 | 4.558 | 3,4 -196+4 
5.13] Fs | 3-49 | 54.7 1.13 | 9,4 049+ 5 
5.17] Gc |} 1.78 | 61.8 1.24 |10,A 058+ 5 
1.50 F6 10.0 | 5.70 | 0.26 4, A .060+ 4 
3 73! “oO, | 2.82 19.8 0.94 |11,A 04644 
9.64] M4+, M6] 1.5 44.52 | 2.362 | 1,A 256+4 
5 Sol Ga, ... | 3.5: 20.46 | 0.82 or ogo + 6 














EXPLANATION OF NINTH COLUMN (REF.) 


. W. S. Finsen, ‘‘Catalogue of Visual Binary Star Orbits,’’ Circ. Union Obs., 91, 24, 1934. In most 


cases the same orbit is given in R. G. Aitken’s recent book, The Binary Stars, p. 284, Table 1. This last 
table is complete to September, 1933. 
this section was written. 


> 


3 
4. 
5. E. Hertzsprung, B.A.N., 7, 330, 1936. 
6. 


7. 
is a parabola: m = 0.0707; g = 


W. H. van den Bos, Circ. 
. K. A. Strand, Annals Leiden Obs., 18, Part II, 1937. 


Union 


Baize’s catalogue (Bull. Astr., 10, 273, 1937) was received after 


Obs., No. 98, 1938. 


W. J. Luyten, “Investigations of Binary Stars,’’ Pub. Minnesota, 2, Part I, 1934. 


G. P. Kuiper, unpublished. The orbit is nearly circular in appearance; i = 0°, e = 0.04. 
H. N. Russell, A.J.,45, 95, 1936; also W. S. Finsen, Circ. Union Obs., No. 95, 230, 1936. Finsen’s orbit 


27175. 


Russell’s orbit represents the measures better than does Finsen’s. 


In computing the mass, Russell’s orbit has been given double weight. 
8. P. Rudnick, A.J., 43, 164, 1934. 
9. W. J Luyten and E. Ebbighausen, ibid., 45, 54, 1036. 
1o. R. C. Huffer, Ap. J., 80, 270, 1934. 
11. W. J. Luyten and E. Ebbighausen, A.J., 44, 110, 1035. 


system),' corrected by the amount —o™22 + 0.11C.'5 Other mag- 
nitude data were first reduced to the PD system with the aid of Zin- 


4 Veriff. Bamberg, 2, 1926. 


SF. W. Seares, Trans. 1.A.U., 4, 142, 1932; C. Payne Gaposchkin, Harvard Ann., 


89, 105-107, 1935. 
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ner’s tables, and then to the IPv system.” The magnitudes of 0, Eri- 
dani BC, BD —8°4352, and Kriiger 60 are my own determinations. 
For uw Herculis BC, the reduced Harvard magnitude is 9.79 mag., 
and my own measures give 9.82 mag. 

For convenience, the correction IPv — PD = —0.22 + 0.11C, 
which seems to be well established, is tabulated against spectral type 
in Table 2. The color index C is defined as IPg — IPv, which 


makes C = 1.08 for gKo if C = o for Ao. 




















TABLE 2 
— i ] a ee RENEE T RS | — ae 
Spec. | IPv—PD ] Spec IPv —PD | Spec | IPv —PD 
| 1] 1] | 

0... | —o.24 || dGo —0.17 | gGo =0:..15 
Bo | 24 | Gs 16 Gs 13 
B8.. 23 | Ko | .14 Ko 10 
Ao.. | 22 K2 12 K2 09 
Ax: 21 5. 09 .5 05 
AS | 20 dMo —0.05 gMo —0.03 
Fo.. .20 
F5 | 0.18 | 








The fifth column of Table 1 gives the spectral type either of the 
two components separately or of the combined light. For A stars 
the HD classes were used, and for most of the other stars spectral 
types were obtained by Dr. Morgan or the writer."’ In particular, 
the writer wishes to acknowledge Dr. Morgan’s classification of the 
spectra of the components of Capella. Three-prism spectrograms 
were used which were compared with similar spectrograms of other 
F and G giants. Dr. Morgan estimates Am(pg) = 0.0, which is in 
good agreement with the small visual difference found interfero- 
metrically. 

For a Centauri B the published spectral type is K5; but since 
this value led to an abnormally large radius for this star, the writer 
has asked Dr. Shapley’s permission to compare the Harvard spectra 
of this star with those of 61 Cygni (dK5), € and 4o Eridani (dK2), 

6 The valuable compilation by C. Payne Gaposchkin (Harvard Obs. Mimeograms, 
3, 1938), was received after this paper was practically completed. It was used for 
Table 8 (except y CrA, which seems to be in error). 


17 Cf. pp. 451 and 463 of the preceding article. 
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and 7 Ceti and & Bootis (dG8). The result, K1, is quite well deter- 
mined and makes the star a normal dwarf."* The writer wishes to 
thank Dr. Shapley and Miss Hoffleit for their kind co-operation. 

The sixth column of Table 1 contains the magnitude differences 
between the components. The majority of the determinations are 
my own, but ten Harvard values were available and these were in- 
cluded. 

The seventh, eighth, and ninth columns give the relevant orbital 
data, including the references explained at the bottom of the table, 
and the quality Q of the orbit. Q = A indicates an orbit as good or 
nearly as good as may be obtained for the pair in question; B means 
a more uncertain orbit, usually due to the fact that the measures 
cover the orbit only partly; and C means that the orbit is quite un- 
certain. 

The list of references illustrates the activity in orbital determina- 
tions since the appearance of the ‘“‘Catalogue of Orbits,”’ by Finsen 
in 1934. Of particular interest is Strand’s publication on six classi- 
cal binaries (five of which are used here), in which many photo- 
graphic measures are incorporated. 

The last column gives the trigonometric parallax, in all but a few 
cases the same as in Schlesinger’s Catalogue of 1935. The few more 
recent determinations available were included. 

The data on mass ratios for the binaries of Table 1 are collected 
in Table 3. A most important source of information is the new Boss 
General Catalogue, which contains the mass ratios that can be derived 
from meridian observations. In addition, van de Kamp has recently 
determined several mass ratios photographically, each based on at 
least fifty plates. 

One word should be said about mass ratios determined for binaries 
that are not separated on the plates used, or, if observed visually, 
are not resolved in the meridian circle. These binaries are ¢ Herculis, 
99 Herculis, HR 7162, 7 Cygni, and 85 Pegasi. The result of the 
finite size of the image (caused by imperfect seeing and optical ef- 
fects) will be that in such cases not the brighter star is observed but 
the center of light. The position of this center may be computed if 

18 The result is in agreement with that just published by Miss Hoffleit (Harvard 


Ann., 105, 57, 1937). 
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TABLE 3 

















| | | || | 

Star —— | Authority | Ref. Star =~ | Authority Ref 

| miTm: i] mirm | 
n Cas.. lo 405 | Boss GC I | © Her. «| 0:372 | Boss GC Ee 
| .385 | Strand 2 || 382 | Vande Kamp| 7 
o2 Eri BC} .31 | Van den Bos 3 || 26 Dra. 38 | Boss GC I 
a Aur 443+] Reese, Sanford] 4 || 70 Oph..| .478 | Boss GC } 4 
a CMa..| .282 | Boss GC I || .416 | Comstock : 
326 | Vande Kamp| 5 || .427 | Strand | 2 
a Gem. || 64 | Boss GC 1 |} oo Her 265:| Boss GC 4 
a CMi. 235 | Boss GC | 1 || HR 7162 45 | Vande Kamp} to 
y Vir | 508 | Boss GC | 1 || PAVE. . .328 | Van de Kamp | x4 
¢U Ma Ri} 807 | Hadley | 6 || Kr. 60...) .359 | Van de Kamp} 12 
a Ce | .446 | Boss GC | 1 | 85 Peg.. | 52 Boss GC I 
£ Boo | 461 | Boss GC | 1 || | o.st Van de Kamp |} 13 

44 Boo | 791 | Boss GC Lee | | 


| 





EXPLANATION OF FOURTH COLUMN (REF.) 


1. B. Boss, General Catalogue of 33342 Stars, 1, Appen. II, 1937. ¢ Her: B (observed) 
= 0.313; Am = 3.00; B (corrected) = 0.372. 26 Dra: since the separation reaches 
1"7, the meridian positions must refer to the primary, and B needs no correction. 99 
Her: B (obs.) = 0.223:; Am = 3.40; B (corr.) = 0.265:. According to a personal com- 
munication by Dr. Jenkins, the weight of this determination is small. 85 Peg: B (obs.) 
= 0.48; Am = 3.5; B (corr.) = 0.52. 

Ann. Leiden Obs., 18, Part II, 75, 1937; the probable error is +o0.o019. 

3D AGN qia5.132, F020. 

4. Lick Obs. Bull., 1, 35, t901; Pub. A.S.P., 34, 179, 1922. Value is approximate. 
A.J., 45, 124, 1936; 59 plates were used, dating from 1920 to 1933. 

6. Pub. Michigan Obs., 2, 101, 1915. 

A.J., 44, 83, 1935. From 50 plates between 1919 and 1934, B = 0.323+ 0.015 
was found, which, corrected for Am = 3.00, gives 0.382. 

8. Ibid., 32, 157, 1920. Based on relative micrometer positions of three field stars, 
separated by the present writer from Comstock’s values obtained from meridian 
positions. The probable error is +0.025. 

Ann. Leiden Obs., 18, Part II, 135, 1937. Value derived from radial-velocity 
observations published by Berman. The probable error is + 0.028. 
A.J., 46, 36, 1937. From 54 plates between 1919 and 1936, B (obs.) = 0.288 
0.045 was found, which, corrected for Am = 1.78, gives B = 0.45. 

11. Ibid., 45, 121, 1936. From 53 plates between 1916 and 1935, B (obs.) = 0.257 
0.024. Corrected for Am = 2.82, B = 0.328. 

12. Ibid., 47, 1, 1938. From 57 plates between 1916 and 1936, B = 0.359 + 0.007. 
(The difference between the x- and y-solutions suggests a somewhat larger uncertainty.) 
13. Personal communication. From measures in the x-direction B (obs.) = 0.476 
0.012; from the y-direction B (obs.) = 0.419 + 0.034, weighted mean B (obs.) = 

0.470 + 0.012; Am = 3.5; B (corr.) = 0.51. 
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the magnitude difference, Am, is known. The center will be on the 
line connecting the components, at a distance [L./(L, + L.)]s from 
the brighter star, where L, and L, are the luminosities of the com- 
ponents, measured in the wave length considered, and s is their sepa- 
ration. It follows that the observed (or apparent) values of B = 
m.,/(m, + m,) will have to be increased by L./(L, + L.), in order 
to give the true values of B.” This correction has been applied in 
Table 3, with the details given in the footnotes. In Table 4 the cor- 
rection (B — B,,,) is tabulated against Am. It appears that Am, 
except if larger than 3 mag., must be known quite accurately in 
order that significant mass ratios may be obtained. If Am is well 





TABLE 4 
| | | 
Am B—Bobs Am B—Bobs 
| 
0.0 ©. 500 | vie ©.OgI 
O.5 387 i} 3.0 059 
1.0 235 it 2.8 038 
nS 201 4.0 025 
z.© 0.337 5.o 0.O1O0 


determined, there is no difficulty in obtaining mass ratios even for 
close pairs with small Am. 

We have not attempted to include in Table 3 all mass ratios that 
have been published, but only those that seemed most reliable. Sev- 
eral extensive investigations based on meridian positions have now 
been replaced by the GC mass ratios. Particularly when the whole 
orbit has not yet been observed, as in the case of » Cassiopeiae, 
a Geminorum, y Virginis, and & Bootis, the meridian observations 
should give a more reliable determination of the mass ratio than 
even the best photographic measures extended over a period of less 
than twenty years. This is obvious from the fact that in such cases 
the accuracy involves the square of the time interval used. 

The binaries for which the mass ratios are known yield two abso- 
lute mass determinations each, free from assumptions. The deter- 
mination of mass ratios is therefore most important. 


19 Cf. E. Silbernagel, A.N., 233, 168, 1928; P. van de Kamp, A.J., 46, 37, 1937. 
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If the mass ratio is not known, and the components have nearly 
equal magnitudes, we shall assign half the total mass to the average 
of the two magnitudes. This has been done for nine pairs in Table 1. 
In these cases each binary yields only one absolute mass determina- 
tion. 

Only two pairs, HR 6416 and 6426, are left for which Am is not 
small, and no mass ratio is available from observations. Here we 
shall use the mass ratio corresponding to the observed Am, adopting 
the slope of the mass-luminosity relation found in this paper. We 
shall use the primaries of these two systems only. 

For the systems in Table 3 with more than one determination of 
the mass ratio, relative weights must be assigned to these determina- 
tions. The weights of the two values of » Cassiopeiae are probably 
in ratio about 1:2; hence, B = 0.392 has been adopted. For Sirius, 
van de Kamp’s photographic determination extends over thirteen 
years near apastron, where the curvature of the relative orbit is 
slightly over 1’. The size of the whole orbit, 14”’, is covered by the 
meridian observations (for simplicity, the comparison is made in the 
relative orbit). The value B = 0.295 has been adopted, which gives 
weight 0.7 to the GC determination. 

For ¢ Herculis the agreement is better than could be expected. In 
time the photographic determination will be very accurate. 

Three determinations are available for 70 Ophiuchi, with approxi- 
mate weights 2:1:1, making the average 0.450; this value is very 
well established. 

Table 3 contains two anomalous values: for 44 Bootis and 85 Pe- 
gasi. Taking the duplicity of the companion into account, the slope 
of the mass-luminosity relation would give B = 0.61 for 44 Bootis. 
In view of the relatively small part of the orbit covered by the ob- 
servations, the observed value of B can probably not be considered 
contradictory to the value B = o.6r. 

The case of 85 Pegasi is perhaps more puzzling. The GC value of 
the mass ratio is based on 0.400 from the x-direction (weight 0.71) 
and 0.669 from the y-direction (weight 0.29). Van de Kamp’s de- 
termination was described in note 13 of Table 3, and is of high 
weight. The observations as they stand would indicate that the pri- 
mary is of abnormally low mass and the secondary of abnormally 
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high mass; but if the parallax of the system should be too large by 
twice its probable error, the primary would be nearly normal and 
the large mass of the secondary might be explained by assuming it 
to be double. 

The data of interest in connection with the mass-luminosity rela- 
tion are found in Table 5. Explanation of the columns is given be- 
neath the table. 

We have not attempted to give numerical estimates of the un- 
certainties of the mass ratios used. Qualitative data are given of 
pages 480 and 482. The weights in column 8 will be upper limits for 
the masses in the tenth column; the secondaries will be of lower 
relative accuracy than the primaries unless the components are near- 
ly equal. For the nine objects where the average component is used 
in column 1o, the weight remains identically the same as given in 
column 8, of course. 

Table 6 gives the quantities of theoretical interest for the stars 
of Table 5 having weight 3 or more; a Geminorum and & Ursae Ma- 
joris are omitted because the components are spectroscopic binaries, 
and HR 6426 is omitted because the spectral type refers to a blend of 
the two components. Furthermore, faint and close companions, such 
as that of Procyon, cannot be-used because no spectral type is 
known, whereas for Kriiger 60 B no reliable estimate of 7, can be 
made. 70 Ophiuchi has occasionally been suspected to be a triple 
system, but the evidence has never been entirely convincing, and 
Strand’s discussion of accurate photographic measures gives no in- 
dication of the presence of a third body. 

The masses of Table 6 are taken directly from Table 5; the lumi- 
nosities follow from the bolometric magnitudes of Table 5, in con- 
nection with the bolometric magnitude of the sun, 4.62.2? The value 
A log T, = log T./© follows from the temperature scale of the pre- 
ceding article," except for dwarfs later than M2, for which no radio- 
metric measures are available. For these stars we shall adopt the 
temperature scale for the giants, on the basis of the following con- 
siderations. 

As we have seen, the heat indices (or bolometric corrections) are 
nearly the same for gM2 and dMz2; the color indices are also nearly 


20 Cf. p. 438 of the preceding paper. 2" Table 13, p. 464. 
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the same.” Hence, the effective temperatures are probably nearly 
the same, as we have assumed in the preceding paper. 

When a spectrum (extending from 4000 A to 6600 A) of BD 
— 8°4352(dM4) was compared with spectra of M4 giants, it was found 
that although the red 770 bands correspond (as they should, since 
the classification is based on them), the 77O bands in the yellow, 
green, and blue are considerably weaker. Dr. Morgan has found the 


TABLE 6 




















| | t 
| | - | E 
Star logm | log L | A log 7 log R | logp | loge 
} km/sec 
E ae : | : | | } 
nCas A —o.14 | —0.09 | +0.021| —0.09 +0.27 | 4-47 0.56 
B — .33 | —1.16 1 — .166| —0.25 | $0.56 4.60 0.53 
o2 Eri B == 488i 722-20 | + £303] 2.74 |) 475302 | 7 86 Esse 
, | — -70 | —1.96:| — .305] —0.37 | +0.56] 4.48 | 0.30 
a Aur A 1 + .62 | +2.08 | — .082] +1.20| —2.84 | 2.65 | 0.17 
B + .52] +1.90 | + .066} +0.82 hI) 3.58 6.22 
a CMa A + .37 | Ftp | + .273] +0.25 | —0.23 | 4.31 0.84 
B > ET = 8. | 5 ee | +5.08 | 7.72 7G ae 
a CMi A + .17 | =FO-70 | = 4077) “e122 —O.g6 | 4.06 0.56 
eae | ail | : » | ; 
¢ UMaA + 41 | +1.48 | + .229} +o.28 | —o.28 | 4.29 | 0.86 
aCen A | + .04 | +0.10 | — .o18} -+-0.09 | —0.07 | 4.30 | 0.57 
B — .06 | —0.43 — .078] —o.006 | +0.27 | 4.50 | 0.64 
— Boo A > 2) 8.52) — .off 068) +0.25/| 4.49 | 0.63 
B......] — .12 | —0.83 | — .157| —0.10 | -+0.34 | 4.52 0.61 
¢ Her A ..| — .02 | +0.59 | + .o10] +0.29 | —0.74 | 3.84 0.31 
F ee | | | 
—8°4352 AB. = 445) = 1ago! |) => -290) 0.72) 0700" 4524 0.30 
uw Her BC. 35 | —1.37 | — .290] —o.10 | +0.12] 4.30 0.36 
70 Oph A — .05 | —0.38 | — .078) —0.03 +0.20} 4.46 0.61 
B — .13 | —o.86 | — .136] —o.16 | +0.49 | 4.62 | 0.68 
Kr.60 <A =9,60 | —13.77 | ~0.298| —0.29 | $0.47 | 4.42 0.31 
BUN o.ciec os veewaclel 0.00 | 0.00 | 0.000 0.00 | +0.15 | 4.44 0.64 
| | | 





same phenomenon even more strikingly in Barnard’s star (dM6), 
as compared with Ms5—M6 giants. Another difference between the 
spectra of giants and dwarfs of type M is, of course, the broad ab- 
sorption feature near Ca 4227, present in dwarfs but absent in 
giants. If the effective temperatures of M giants and M dwarfs 
were about the same, the spectral differences mentioned would tend 
to give the dwarfs the larger color index (for two reasons) but the 
smaller bolometric correction. The bolometric corrections are not 
known for the later M dwarfs, but the color indices are, approxi- 


22 Seares, Pub. Nat. Acad. Sci., 5, 232, 1919. 
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mately, from Seares’s determination for Barnard’s star?3 (dM6, 
( = 1.76) and from the approximate value for Wolf 359 (dM8, 
C = 2.0).24 The M giants are known to become slightly bluer (as 
measured by the ordinary color index) with advancing type.” The 
color difference is, therefore, in the direction to be expected from 
the appearance of the spectra. Whether or not the full difference 
would be accounted for (which would mean that the effective tem- 
peratures of giants and dwarfs are equal) would be difficult to de- 
termine. Another approach to the temperatures of the M dwarfs 
would be to estimate the distortion of the color index of Barnard’s 
star by the 770 absorption. Dr. Morgan estimates that the cor- 
rected value may be about 2.0 mag., with a considerable uncertainty. 

With this (uncertain) value, the effective temperature in our sys- 
tem may be computed if color indices for two other spectral types 
are adopted. We take C = 0.00 for Ao, and 1.22 for dK5+ (being 
the mean of Seares’s average for the components of 61 Cygni, 1.18”, 
and W. Becker’s photoelectric color index reduced to the interna- 
tional system, 1.25). With the aid of Table 137" we then find ¢,/T = 
5.19, which happens to be very near to the value 5.22 for gM6 given 
by Table 13. Much significance cannot be attached to this result; 
only radiometric observations of Barnard’s star will solve the prob- 
lem of the effective temperatures of the later M dwarfs. 

Since the temperatures for giants and dwarfs seem to be nearly 
the same for M6 as well as for M2, the temperature scale of the giants 
will be used also for the dwarfs. 

The mean densities, p, and the values of the surface gravity, g, 
given in Table 6 are in c.g.s. units. The variation of g along the 
main sequence is seen to be quite small. The mean densities vary 
from o.5 to 0.6 for Mizar A and Sirius A to about 3 for K and M 
dwarfs. 

The last column gives the Einstein gravitational red shift. It ap- 

23 Seares, Pub. A.S.P., 28, 281, 1916. 

24 Contr. Mt. W. Obs., No. 356, 3, 1928. 

5 For the reasons given, the treatment of the M dwarfs in papers by Opik and 
Gabovits, who assume equal 770 absorption in the visual region for giants and dwarfs, 
and who adopt a constant difference between the color indices (0.3 mag., the dwarfs 


being bluer), is not correct, although the errors in the two assumptions tend to com- 
pensate. Cf. Pub. Tartu, 28, Nos. 3 and 5, 1935-1936; 30, No. 1, 1938. 














488 G. P. KUIPER 


pears to diminish slightly along the main sequence. Only for the 
white dwarfs is the red shift appreciable, although in statistical 
studies dealing with the F, G, and K stars it should not be over- 
looked that the standard radial velocity system is corrected by an 
amount of 0.64 km/sec, the red shift of the sun, in order that the 
sun would show zero radial velocity. The giants of these spectral 
types should, therefore, show a small negative K term. 

The two white dwarfs of Table 6 will be discussed, together with 
similar objects, in a forthcoming paper. 

A theoretical discussion of the masses, radii, and luminosities de- 
rived here will not be given, since Drs. Chandrasekhar and Strém- 





























TABLE 7 

ADS Mpv Sp Mpbol log L | A log T. log R log m Wt 
3264. 3.02 Fo 2.90 +o0.69 | +0.118) +o.11 | +0.07 13 
3483 3.66 Fs 202 + .40] + .054] + .09 | + .04 I 
3135 4.48 F6 4.44 + .07 | + .046] — .06 | — .195] 3 
3169 4.52 F8 4.47 + .06 | + .032] — .03 | — .21 ; 
3475. 4.50 F7 4.54 + .03 | + .039] — .06] — .19 Z 
3210. | 5.65 G5 es —0.37 | —0.028] —0.13 | —0.44 2 








gren intend to undertake such a discussion with the methods they 
have developed. 

The values themselves are shown graphically in Figures 1 and 2. 
Figure 1 contains also objects of lower weight, taken from Table 5, 
and the binaries in the Hyades, given in Table 7. The latter table 
is a revision of a table published previously ;* for ADS 3135 a new or- 
bit was published by van den Bos,??7 who used measures up to date. 
As mentioned in the previous paper, the mean of the six objects in 
the Hyades should give a reliable determination of the hydrogen 
content, but deviations from the mean are not considered signifi- 
cant. In due time the pair 68 Tauri (M,, = +1.4, Mio = +0.9); 
discovered recently by the writer,”* will be valuable in settling the 
question of possible variations of composition with absolute magni- 
tude. 

The companion to Procyon was entered in Figure 1 as a probable 

26 Ap. J., 86, 176, Tables 4 and 5, 1937. 

27 Circ. Union Obs., No. g8, 1937. 2 Pub. A.S.P., 49, 341, 1937. 
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white dwarf for the following reason. Observations show that the 
correlation between absolute visual magnitude and spectral type is 
close for the lower part of the main sequence. The deviation of 
Procyon B from the empirical mass-luminosity relation may there- 
fore be computed under the assumption that the star is a normal 
dwarf. The deviation is in that case 2 mag., or over 0.4 in log m. 
Until other large deviations are found in this region of the mass- 
luminosity diagram, we may conclude that Procyon B is probably 
not a normal red dwarf. 

The dots in Figure 1 define an average empirical mass-luminosity 
relation which is useful in the discussion of such problems as that 
presented by e Aurigae. Since, however, different selection factors 
operate for the different groups of stars used, it is doubtful whether 
this mean empirical relation has any definite physical significance; 
no interpolation-curve was therefore drawn. Some stars, such as the 
Hyades and ¢ Herculis, are seen to be brighter than the average star 
of the same mass; these differences are considered real. Other dif- 
ferences, particularly among the stars of low accuracy, are probably 
spurious. The conclusion, therefore, is that we are not dealing with 
only one mass-luminosity relation but with many such relations, in- 
volving still a parameter different from the mass. Strémgren has 
identified this parameter with the hydrogen (or hydrogen and he- 
lium) content. This conclusion is based on another result, by Stroém- 
gren and Chandrasekhar, that the stars of mass smaller than five to 
ten times the sun are built essentially on the same model, so that 
changes in model cannot be held responsible for the differences in 
luminosity for stars of the same mass. For the more massive stars, 
however, Chandrasekhar has found very striking differences in mod- 
el, which must be largely responsible for the scatter in the mass- 
luminosity relation in that region. 

For the stars with masses smaller than 0.6© the slope of the mass- 
luminosity relation changes;?? the exact slope will depend, of course, 
on the bolometric corrections used, but there is no doubt that the 
change is real. The explanation may be simply that the stars Kriiger 
60 and o, Eridani C have a small hydrogen content, as would follow 


29 The change found now is better determined than in Fig. 1 of Ap. J., 86, 194, 19373 
cf. also Parenago, loc. cit., Gabovits, Pub. Tartu, 30, No. 1, 1938 (see n. 25). 
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from an extrapolation of Strémgren’s tables. A more detailed theo- 
retical investigation of these low-temperature stars will be needed 
before this conclusion can be accepted. If it holds, it may be expect- 
ed that stars will be found which deviate strongly from the mean re- 
lation in Figure 1, toward the fainter side. It is of great interest that 
the white dwarf 0, Eridani B (also) has a very small hydrogen con- 
tent, which follows from the data in Table 6 and Chandrasekhar’s 
theoretical relations. Castor C, connected with the A-type star Cas- 
tor AB, has probably a fairly high hydrogen content. Furthermore, 
the brighter component of u Herculis (A) lies well above the main 
sequence and may therefore be expected to be relatively poor in 
hydrogen. It would, therefore, not be surprising to find BC some- 
what bright for its mass. It is also noteworthy that both components 
of m Cassiopeiae deviate in the same sense from the mean relation of 
Figure 1. 

These remarks are all of a provisional nature, but they indicate 
the direction in which future research may be continued. The accu- 
racy of the observations (parallaxes, temperatures) will set the limits 
to which discussion is profitable. 

In conclusion, in Table 8 the visual binaries are collected which 
promise to be of future interest in the study of the empirical mass- 
luminosity relation. They are of special importance in parallax and 
double-star observations. Spectral types marked with an asterisk 
are my own determinations. Of particular importance is Wolf 424, 
which will provide a significant extension of the mass-luminosity re- 
lation beyond Kriiger 60 B. The value of the spectral parallax of this 
star has been unusually uncertain because of the great steepness of 
the main sequence in this region. A change of M8 to M7 corresponds 
to a change of 1.7 in M, or a factor 2.2 in the spectral parallax. That 
means that even estimates of spectral type with a probable error of 
only one-half of one-tenth will still give quite uncertain results for 
types later than dM6.%° 

The limits set in the compilation of Table 8 are necessarily arbi- 
trary. If the evidence available indicated that the parallax is (or 
may be) in excess of 0050, the pair was included, provided the pro- 
jected separation between the components did not exceed about 


3° Ap. J., 87, 593, Fig. 1, 1938. 
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sixty times the parallax (this limit will already include many pairs 
of periods over five hundred years). For stars of very large parallax, 
the limit was taken somewhat wider, since a much greater relative 
accuracy may be obtained for such pairs (on account of the larger 
angular scale of the orbit and the applicability of the photographic 
method of observation; further, mass ratios may be determined for 
such pairs from relatively short arcs). Boundary-line cases are en- 
tered in parentheses. 

The magnitudes, map, are reduced to the IPv system wherever 
possible; values depending on estimates are indicated by a colon. 
Of the magnitude differences, Am, however, many values depend on 
estimates. Two systems, ADS 1322 and 9476, were excluded as prob- 
ably optical, the spectral parallaxes derived by the author being very 
much smaller than the computed dynamical parallaxes. 

Not mentioned in the table is Ross 614, found by Reuyl to have a 
variable motion; if the pair can be resolved at maximum separation, 
it would be a most interesting addition to Table 8. Other stars of 
particular interest are ADS 9655 and 16644. All these objects de- 
serve a prominent place in parallax and double-star observations. 

Other objects not mentioned are those for which interferometer 
observations could be combined with the spectrographic observa- 
tions. They will be few in number but of great importance individu- 
ally (luminous stars). 

Finally, Ross 492" and a pair announced by Luyten,*” may appear 
to have parallaxes in excess of o%o050. 


SOME SELECTED SPECTROSCOPIC BINARIES 

t. After a discussion of the visual binaries it is natural to collect 
all the data on spectroscopic binaries which have a bearing on the 
problem of the mass-luminosity relation. These data fall into two 
groups: (1) double-line spectroscopic binaries that are also eclipsing 
variables (for these objects the masses and radii may be found and 
the luminosities computed by means of the stellar temperature 
scale); (2) double-line spectroscopic binaries for which the mass ra- 
tio may be measured and the magnitude difference found from 
photometry in the spectrum. 


31 Pub. A.S.P., 47, 282, 1935. 32 Harvard Bull., No. 906, 1937. 
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TABLE 8* 
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ADS 48 
(ADS 61) 
PASP 47,282 
ADS 433 
ADS 490 

— 69°44 


—30°529 AB 
AB, C 

ADS 1733 

ADS 1865 


a For 
468°278 BC 
(ADS 2995) 

68 Tau. 
ADS 3248 


ADS 3321 CD 
HR 1504 
104 Tau 
a Aur Hh 
(ADS 3900) 
—51°1540 
BDS 3112 AB 
AB, C 
HR 2162 
—46 3046 
ADS 6487 
HR 3430 
ADS 7067 
« UMa BC 
A, BC 


10 UMa 
(a. UMa) 
(ADS 
ADS 7284 
vy Arg 


11 LMi 
—60°2911 
Wolf 424 

78 UMa 

42 Com 





(ADS 8841) 
ADS 8887 
ADS 9031 
ADS gogo 
(ADS 9446) 


ADS 9544 
n CrB 

ADS 9716 
wy Ser 








R.A. DECL. m 1pm SPEC. é P Tn Pe SOURCI 
Ab m(d) 
ir dy. sp. 
ohoomy +45°16’ 8.1 0.1 *K6, K6 os 89+5 120 79 930 
© o1.0+57 53 6.0 0.8 G3 106 14+8 53 39 =. 30, QI 
© 07.7+68 47 12.2: 0.3 *M5 7 56+7 
© 26.0+66 42 10.2 1.9*M3 3 103+7 690 
© 30.11— 409/ 5.2/0.6 F7 7 52+8/ (51) 46 930 
I 11.6'\—69 21 | 7.4:| 0.6 Gs 128 52 98 
I 30.4 —30 25 2/0.2 Gs5 1.6 36+9 56 g! 
. 4.¢ I.5 59 03 
2 11.1 —18 42 | 7.8 | 1.0 *K3 2 74 (42) 930 
2 22.5'+ 3 590 | 8.6:| 0.2) *K6 25 61 98 
3 07.5/—20 23 | 4.2 | 3.5) Fs I 72+5 55 | 87 93 
3 38.2 +68 10.8 0.0 *M2+ 0.4 5440 60 
4 00.9 +37 49 7.1 | 1.3 *K2 I 36+6 58 44 930 
419.7, +17 42 | 4.4 | 4.0 A2 1.5 Hyad 25 
4 23.3 +15 56 6.6 1.0 F8 40 Hyad 28 98 
4 30.2 +16 19 11.2:; 2.4* K6 1.7 Hyad 
4 38.6—59 08 | 6.4 | 0.2 Go I 50 03 
5 O1.5 +18 31 | 5.1 | 0.0! Gr o.1| 2+?) 58+8 16 
5 10.0+45 44 10.1 | 3.5*Mr+ 2 73+3 © Intf.=63 
5 14.1j— 3 11 | 7.9:| 2.8% Ky, 4 66+5/} 110}; 58 930 
2;—-51 08 | 9.3:| 0.2)| Ko I 110 93 
6|\—31 03 | 7.6 | 0.4/* Ks < 53+11 79 | (78) 93 
o.1 2 51 93 
6 02.2,—48 27 | 6.5 | 0.4 Gs 2 52 93 
7 14.6;\—46 49 | 7.6:| 1.1| Ko I 73£9| 14 03 
7 §2.9\— © 33 | 8.3:| 5.0% Ka 3 81 (69) 485 
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8 46.0+71 11 | 8.6 | 0.2*K6, K6 3 934 go 6¢ 930 
8 52.4+48 26 10.2 0.0 Mr 35+ 66+6 30 
3:2 aoe AS 6 gI 63 930 
8 54.2\+42 II | 4.0 | 2.0) F2 0.6 72+6 66 
9 01.6|+67 32 | 4.9 | 3.3) F4 1.5 53+6! 40) 50 g! 
9 07.7\+53 07 | 7.2 | 0.0\*K6, K6 19 162+3 |(240) (152) 930 
9 12.0+29 00 7.2 | 0.0*K3 34 60+ 8 57 48 930, 91 
9 26.8|\—40 02 | 3.7 | 1.2) Fo 35 7ot9 61 55 930, OI 
. 
9 29.7;+36 16 | 5.4] 7.5| Ko 2 o5t5 gl 
I 20.3 —61 06 | 7.4:| 1.2 Ks5 234 67+13 104 oI 
2 28.4+ 9 34 11.9 | 0.0*M7(+) I (220) (320) t 
2 56.4+56 54; 5.0 | 2.4 Fo go UMa(?) 47 24 BAN 305 
3 05.1 +18 04 | 4.7, 0.0 F4 26 S77 56 43 | 930, OI 
3 11.9 +17 6.4 | 2.7) Es 5 68+7 94 79 930 
3 18.9+29 45 8.9 | 0.4 *K6 0.6 67+11 19 46 930f 
3 44.5+27 29 7.0 0.4 K6, K6 156 60+9 79 66 gl 
3 58.5 +46 490 | 9.3 | 0.1 *M1, Mr 3 79+6 05 
4 51.6—20 58 | 5.4 2.2) Ks, M2 19 172+4 (200) (174) 03 
5 08.8— 058 6.5 0.0 G8 O.1|' 2t°? Gots , 42 
5 19.1+30 39 | 5.1 |} 0.3} Fo 42 68+8 (57) 44 930 
5 32.5 +40 08 | 6.7 3 *K3 57 46+0 47 54 § 
§ 390-0+ 2 50 5.9 6.3 G5 3 46+10 8g 52 485 





* Spectral types 


erences for the dynz 


marked by an asterisk (*) are my own determinations. The last column gives the ref- 
amical parallaxes; 930=A.J., No. 930; 91, 93, 95, and 98 are the numbers of the 


Union Observatory Circulars; 451 and 485 are the numbers of the Lick Observatory Bulletins. 

+ Two good spectra available for this star give M7 + and M7, respectively. The star is, therefore, defi 
nitely somewhat earlier than Wolf 359. 

t The spectral parallax is about o%06. It seems unlikely that the parallax should be much less than 
o"o05 or 004, as the dynamical parallax would indicate. The value of the latter was confirmed by the in- 
clusion of recent measures. Hence, the orbital inclination must be high (in which case the binary should 
soon close in rapidly), or the mass of the binary is unusually small. Additional trigonometric determinations 


will be valuable. 


§ Derived from Pitman’s orbit, Pub. A.S.P., 46, 196, 1934. 
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TABLE 8—Continued 











1000 PARALLAX 





STAR R.A. DECL. m 1B 4m SPEC d”’ P SOURCE 
Ab r(d) 
lr dy. sp 
+-45°2505 i7hogm2 +45°50’ 9.4 | 0.4*M3+4 <I 143 +6 115 
36 Oph 17 09.2 —26 27 | 4.40.0 Kr, K3 1 .| 179+ 5 52 | 174 93 
ADS 11632 18 41.7+59 29 8.5 | 0.8*M4, M4 16 282+4 (240) (282) 930 
>» CrA 18 59.7 —37 12 | 4.20.0 F7, F7 119 62 51 gI 
17 Lyr C 19 03.6432 21 11.2 0.4*M5 0.2 123+6 a6 
Ross 165 19 41.7 +26 55 12.2  0.9*Ms5+ 0.7 116+11 : 
ADS 12889 IQ 41.8+33 22 | 7.7 | 0.2 *K3 I 48+ 5 57 44 9309 
HR 7703 20 04.6 —36 21 | 5.2 | 5.4 Ky 8 177+9 IQ! ; 
+4°4510 20 34-.6'+ 4 37 | 7.9 | 1.0 *K5 0.7 5744 87 
Furuh 54 20 56.2 +39 41 10.0 1.9*M2+ I 97+8 87 
(61 Cyg) 21 02.4+38 15 | 4.9 | 0.7 *K5, K6 26 2909+ 3 (360) (276) 930 
— 58°7803 21 39.4 —58 08 | 8.6: 0.11 Ma 6.8 49+ 5 54 05 
(uw Cyg) 21 390.7\+28 1714.5 |1.4 F6, F3 <1 46+8 66 44 930** 
HR 8501 22 11.7;\—-54 06 | §.5 | §.1| Go 3 80+7 ae 





Another case (n. {, above) where the dynamic parallax is much too small. 

{ Strand (B.A.N., 8, 206, 1937) has computed a provisional orbit for this star and found the dynamical 
parallax o"117. It is almost certain that this value is too large, probably as a result of the fact that the orbi- 
tal inclination found by Strand is too close to go°. The trigonometric parallax, the Mount Wilson spectro- 
scopic parallax, and my spectral parallax (0"050) agree very well. Also Strand’s doubt whether this binary 
is physically connected with 17 Cygni seems unjustified. Proper motion, radial velocity, and parallax sug- 
gest that the connection is real. 

rhe nature of this case is apparently the opposite of cases t and ||, above. 

** The secondary was classified at Mount Wilson earlier than the primary. 


The first group gives fundamental data for the problem of stellar 
structure, and the second group is useful in a study of the scatter 
of the mass-luminosity relation’ and in giving data that may be 
used statistically in a determination of the mean position of the 
mass-luminosity relation (the mean absolute magnitude may be de- 
termined from the spectral types; the mean mass, by substitution of 
the mean value of sin‘ z into the mean value of m : sin 7). 

A review of the available data showed that for many close pairs 
contradictions seem to be present;*4 it was clear that for many ob- 
jects new observations would be required in order to find the inter- 
pretation of such discrepancies. Particularly close pairs with very 
unequal components are likely to present apparent anomalies; the 
systems with nearly identical components, however, are more easily 
interpreted. In the present paper we shall discuss only a few objects 
of the latter class; the discussion of additional spectroscopic binaries 
will be made later by Dr. Morgan and the present writer. 

33 Cf. Pub. A.S.P., 47, 144, 1935; the scatter found in this manner will be a minimum 
value, since the components of binaries seem, at least statistically, to have common 
properties (hydrogen content). Cf. Harvard Bull., No. 903, 1936. 

34 E.g., Wyse, Bull. Lick Obs., No. 464, 1934. 
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2. It is well known that in many classes of light-curves the ratio 
of the radii cannot be found with any great precision. It is therefore 
necessary to determine the ratio of the total brightnesses from the 
relative strength of the two superimposed spectra; in combination 
with the ratio of the surface brightnesses determined from the light- 
curve the ratio of the radii is found. All our discussions of the radii 
will be made in this manner. Fortunately, the only cases which are 
of interest here are binaries which show two spectra. 

3. Asecond question in connection with the determination of the 
radii arises when the pair is very close, so that the components are 
ellipsoidal in shape. We have a good theoretical evidence that the 
stars possess strong central concentrations; only certain classes of 
very massive stars may be exceptions to this rule.** Hence the 
Roche, or point-mass, model should give a close approximation in 
the determination of the equipotential surfaces; the boundaries of 
the components will, in general, be two of these surfaces. Since for 
the present we shall limit ourselves to pairs of nearly identical stars, 
the model considered will have two equal components with envelopes 
of negligible mass, in which the periods of rotation are equal to the 
orbital period. In this case the equipotential surfaces are given by 
the formula* 

ee ee ae he (3) 

Tr Wo as 

The symbols have the following meaning: is the potential for 
which a constant value is adopted; G is the gravitational constant; 
and M is the mass of each component. The orbit is supposed to be 
circular, with radius a; the «y-plane is assumed to coincide with the 
orbital plane, and « = y = o with the center of gravity; 7, and r, 
are the distances of any point (x, y, z) from the centers of the two 
stars; x, y, z refer to a rotating co-ordinate system. 

With the aid of this formula the writer has computed a family of 
curves, {2 = C, both for the xy-plane and for a plane parallel to the 
yz-plane through the center of one of the components. The results 
were plotted on large-scale diagrams from which the ratio of the 

35 Chandrasekhar, op. cit., chap. viil. 


© E.g., Jeans, Problems of Cosmogony, p. 162, 1919. 
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axes, d,, b,, and c,, were read off for a given a,/a (if a is the distance 
between the components). A certain amount of rounding-off was ap- 
plied to the pointed sides of the stars facing each other, since theo- 
retically the intensity of these points (apart from reflection effects) 
should be low. We are interested in knowing the mean stellar radius, 
R = Wa,b,c,, Whereas the light-curve gives a,, or a, and b,. We 
shall therefore tabulate R/a, against a,/a, so that the former quan- 
tity may be found from the latter, which is observed. The quantity 
b,/a, to be expected from the Roche model is also given, since this 














TABLE 9 

a,/a | R/a: b./a, a,/a R/a; b./a: 

| a | 
0.10 | I.00 1.00 1} 0.35 0.945 0.935 
15 | 0.995 0.995 | 40 O15 go 
20 | 0.990 0.99 | 45 885 86 
25 | 0.98 0.98 | 0.48 0.88 0.855 

0.30 | 0.965 0.90 || 

| 





ratio may be obtained observationally. Table 9 contains the results, 
rounded off to the nearest half of o.or. It should be noted that for 
a,/a@ = 0.43, the two components just touch each other; for larger 
values of a,/a there will be a bridge between them, which, as men- 
tioned above, may be expected to have a low intensity (apart from 
reflection effects), so that in the photometric solution the compo- 
nents will still appear separated. 

Before applying the tabular values to actual stars, we may 
consider the results of Wesselink’s discussion of SZ Camelopar- 
dalis.37 From an unusually accurate light-curve he finds 6,/a, = 
0.9323 +.0020 (m.e.) in the uniform solution, for which a,/a = 0.38, 
a,/a = 0.27, and Am = 1.0. The mean value of a;/a (weighed ac- 
cording to intensity) is 0.35. For the darkened solution he found 
b,/d, = 0.9579 + .o013. The theoretical darkening coefficient for 
this star is about 0.1, and we may adopt the empirical value },/a,; = 
0.935 + .002 (m.e.). This happens to agree precisely with the tabu- 
lar value. (Wesselink had already pointed out this agreement with 
the Roche model.) 


37 Dissertation, Leiden, 1938. 
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4. Finally, the reflection effect should be taken into account in the 
determination of the masses. Since the centers of light on the two 
disks do not project themselves precisely on the centers of mass, but 
are displaced toward the center of gravity, the amplitudes measured 
are found too small even if the components are of similar brightness. 
The reflection effect will increase the bolometric brightness on cer- 
tain areas of the stars; in order to find the photographic increase of 


TABLE 10 


CORRECTION OF MASSES FOR REFLECTION EFFECT 























a:/a Alogm a,/a A log m 
OMIO..6. <: 0.000 6.45... ’ 0.021 
co | eee .OO! AO. <: 032 
20. <.: ‘ 003 45. 046 
Ge... sa .007 (0.50). . (0.066) 
0.30... 0.013 
TABLE 11 


FACTOR TO APPLY TO THE CORRECTION OF TABLE 10 

















Spec | Factor Spec. Factor 
OF: eae | 0.45 BO... a 0.95 
Bo 45 ars... LES 
Bs. ae | =~ dGo I.4 
Ao “65 dKo. rs 
A5 3 | 0.8 dKo.. 255 





intensity, the ratio of photographic increase to bolometric increase 
has to be known, which will depend on the spectral type. The sim- 
plest procedure is to tabulate the correction to log m asa function 
of a,/a for bolometric intensities; this correction will then have to 
be multiplied by a factor giving the rate of increase of photographic 
intensity for a given increase of bolometric intensity. This factor is 
readily found from Table 3,or from Tables 9 and 13, of the preced- 
ing paper, and a table of color indices; the result is found in Table 11. 
The correction to log m for bolometric intensities is given in Table 
10; Eddington’s formula for the reflection effect** was used. 


38 The Internal Constitution of the Stars, p. 212, formula 144.3, 1926. 
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The reflection has, of course, also its effect on the determination 
of a sin 7, and hence on R (involving a) and L (involving a?), as well 
as on m (involving a3). The corrections are of less importance than 
those for the masses, but may be taken from the same tables (10 
and 11), after multiplication by 1/3 for log R, and by 2/3 for log L. 

After these preparations, some close binaries with nearly equal 
components will be discussed. 

5. Castor C.—From the discussion in the preceding article? we 
have log R = —o.198. The ratio of the radii follows from the ratio 
of the intensities of the two spectra, 5:4, and the equality of the 
surface brightnesses. We find log R, = —o.18 and log R, = —o.22. 

The direct determination of the effective temperature gave 
log T. = 3.550 + 0.014;9 from the spectral type K6+ and the 
stellar temperature scale*’ we have log T, = 3.562. We adopt the 


mean, 3.556. Hence," log L, = —1.16, log L, = —1.24; Mpo(1) = 
+7.53, Mioi(2) = +7.73. With 7 = 86°0 the masses are given by” 
log m, = —0.204; log m, = —o0.250. With R/a = 0.16, the correc- 


tion for the reflection effect is about +-0.003, making log m = 
—o.201 and —o0.247, respectively. 

6. B Aurigae.—The photographic magnitudes are nearly equal,* 
and also the surface brightnesses.*4 Hence, the radii are nearly equal: 
log R, = log R, = 0.431.4° The direct determination of 7, gave 
log T. = 3.99 + 0.04; from the spectral type Ar we have log T, = 
4.01; we adopt 4.00. Then log L, = log L, = 1.83; Mio.) = 
+0.04. 

With cos 7 = 0.223 and Baker’s orbit the masses are found: 
log m = 0.376 and 0.370, respectively. The reflection effect found 
from a,/a = 0.15 adds less than 0.001 to log m. 

7. wx Scorpii.—From our former discussion*® we take log R, = 
log R, = 0.73. The direct determination of 7, was log T, = 
4.20 + 0.05. From the spectral type B3 of the average component 


39 P. 458. 

4° Table 13, p. 34, of the preceding paper. 

4* Formula 1, p. 475. 

42 Joy and Sanford, Ap. J., 68, 253, 1926. Computed from K, P, and 7. 

43 Baker, Pub. Allegheny Obs., 1, 173, 1910. 45 P. 442 of the preceding paper. 
44 Stebbins, 4p. J., 34, 126, IgII. 4 P. 443 of the preceding paper. 
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we have log T. = 4.27. We adopt 4.23. Hence log L = 3.35 and 
M1 = — 3-76; this value applies to the average of the components. 
Since the mass ratio for this binary is not known, we shall use the 
average mass of the components in connection with the average 
log L. With 7 = 62°0, the average uncorrected mass is log m = 
1.082; the correction for reflection (a,/a = 0.370) is 0.012, and 
log m = 1.094. 

8. V Puppis.—The mass ratio is not known; therefore we shall 
use the average component. The strength of the companion is about 
0.7 of that of the primary.*’ Roberts’ photometric observations be- 
ing visual, the amount of darkening toward the limb must be negli- 
gible; hence Shapley’s uniform solution® will be used. Then J;/J; = 
1.27 = 1/0.79. This value in connection with the ratio of the total 
brightnesses indicates substantially equal radii, as Shapley has re- 
marked. Shapley finds the mean of a;/a to be 0.438, and b;/a; = 
0.813. The Roche model gives };/a; = 0.87. Since the light-curve 
near the minimum depends much more on 0; than on a;, we shall 
use the value of (3b; + a;)/4 as found by Shapley, and compute R 
by means of the Roche model. The result is R/a = 0.373. With 
Shapley’s value of cos 7 (0.219) and Maury’s value of a sin 7 we find 
log R = 0.821. The mean spectral type being Bz (Maury), log 7, = 
4.31, and the mean log L is 3.85. From P, K, e, and 7, we further 
find log m = 1.249. From the computed value a;/a = 0.409 we find 
A log m = 0.016 and A log L = 0.011 ; the corrected values for the 
mean component are log m = 1.265, log R = 0.83, and log L = 3.86 
(Mio. = — 5-04). 

9. Y Cygni.—As Redman” and Dugan*%’ have pointed out, the 
observations for this binary are not all satisfactorily explained (pos- 
sibility of a third body, etc.), so that some caution should be exer- 
cised in using the results. Yet it is doubtful whether the chief quan- 
tities, R and m, are much affected by these difficulties. 

The spectra are of nearly equal intensity (Plaskett, Redman) and 
the minima are of equal depth (Dugan); hence the radii are nearly 
47 A. C. Maury, Harvard Ann., 84, 173, 1920. 

48 Princeton Pub., No. 3, 1915. 
49 M.N., 90, 754, 1930; Pub. Dom. Ap. Obs. Victoria, 4, 341, 1931. 
5° Pub. Princeton Obs., No. 12, 1931. 
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equal. Dugan found R, = R, = 0.206a = 5.86© (uniform solution, 
which will be used since the darkening must be small). The elliptic- 
ity effect is negligible because the width of the minimum determines 
b;, which is nearly equal to R if R/a < 0.25. 

The writer is under obligation to Dr. Morgan, who obtained an 
accurate spectral class for the star from the plates available at the 
Yerkes Observatory. This was of importance, in view of the dis- 
cordant values published: B2 in the Henry Draper Catalogue, B2 by 
J. S. Plaskett,5" O9.5 by Pearce, and Og in the Victoria catalogue of 
radial velocities of O and B stars. Dr. Morgan finds the type to be 
certainly as early as Og, but probably not much earlier than Og. 
He suggests the value Og as the best obtainable from the data at 
hand. 

From Dugan’s dimensions we now have: log R = 0.768; log L, = 
log L. = 4.51; Myows,2) = —6.6. The masses are 17.1 and 17.3 
(Dugan, Redman), log m = 1.233 and 1.238; the correction for re- 
flection to log m is only + 0.002. 

10. AO Cassiopeiae.—Pearce has made four photometric solu- 
tions, corresponding to Am = 0.25, 0.50, 0.80, and 1.11 mag. He 
estimates the ratio of the intensities from the spectrum to be 2.6:1, 
but this may perhaps be an upper limit and we shall adopt the solu- 
tion for Am = 0.80 mag. The mean value of a;/a is then 0.500; of 
b;/a, 0.478. This is not possible for completely separated stars built 
on the Roche model, for which b;/a < 0.38. If the photometric solu- 
tion is correct, it is probable that the components have a common 
envelope; and since the stars are not identical, it is then further prob- 
able that the system is not in equilibrium. This possibility should 
be borne in mind if comparisons with other stars are made. 

It is in this case uncertain what values for the radii should be 
used; since, according to Table 9, b; and R are never very different, 
we shall use the values of 0;: R,/a = 0.549 and R,/a = 0.406. 

In the case of unequal components, the effect of reflection on the 
constants of the system is slightly more complicated than in the 
case of identical components. With the ratio of the radii given and 
with equal surface temperatures (both spectra are O8.5), Am is 

st Pub. Dom. Ap. Obs. Victoria, 1, 213, 1920. 


s? [bid., 3, 297, 1926. Cf. also: M. Giissow, A.N., 207, 321, 1930. 
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found to be 0.65. In order to make Am(pg) = 0.80, we shall have 
to make the temperatures slightly different: A log T, = 0.03, which 
makes Am(bol) = 1.0. 

The reflection effect on the amplitudes, A, is one-third of that on 
the masses, which is found from Tables 1o and 11. For equal com- 
ponents and a;/a = 0.50, the effect on log K is 3(0.066) + (0.45) = 
0.010. With the intensity ratio 23:1, the effect on log K, will be 
about 0.4 X 0.010 = 0.004, and on log K,, 23 X 0.010 = 0.025. 
Applying these corrections to the observed values, we have AK, = 


220; K, = 248; m,/m, = 0.887; log aO/ = 1.625; log Ri = 1.365; 
log R, = 1.234; log Li = 5.97, log L2 = 5.58; Mbowy = —10.3, 
Myoi2) = —9-3; log m, = 1.634, log m, = 1.582; m, = 43.00; 
m, = 38.20. 


In allowing for the reflection effect we have increased the masses 
by a fair amount and have made them more unequal, which is 
more favorable to the difference in dimensions found for the two 
components (since they are probably in contact) and in better agree- 
ment with the difference in brightness. 

However, perhaps the main uncertainty about this system is the 
value of the inclination, which enters seriously in the determination 
of the masses. This is particularly evident if a comparison is made 
with the ellipsoidal systems, 7° Orionis, b Persei, and ¢ Andromedae, 
discovered by Stebbins. A smaller inclination would obviously in- 
crease the derived ellipticity of the components, which is now found 
to be smaller than is possible for two stars which are just in contact. 
It will remain for a new photometric study of this system to settle 
this question. 

The extremely high luminosity of this system is noteworthy (a 
smaller value of the inclination would even increase the luminosity). 
With the provisional values of the constants, it lies well above the 
mean empirical mass-luminosity relation; even a somewhat de- 
creased inclination would not change that feature. This fact was 
interpreted by Chandrasekhar*s to be a result of strong mass con- 
centration, contrary to the massive stars found by Trumpler and 
even stars like V Puppis, which are less centrally condensed than 
Eddington’s standard model. In that respect AO Cassiopeiae seems 
to resemble the stars of lower mass, whereas Trumpler’s stars are 
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really the stars that deviate from a theoretical mass-luminosity re- 
lation based on the standard model. It will be of interest to dis- 
cover which of the two classes of massive stars is the most frequent 
in nature. 

On account of the uncertainties mentioned in connection with 
AO Cassiopeiae, it is very fortunate that a second system is known 
with similar properties. This system is 29 Canis Majoris. 

11. 29 Canis Majoris.—Spectroscopic orbits have been published 
by Harper, Pearce,54 and Luyten and Ebbighausen.** Pearce has 
measured both components and finds the spectral types both O7; 
the secondary is very faint. Furthermore, A, = 217 + 5 km/sec; 
K, = 288 + 7 km/sec; P = 4.3935d; and e = 0.156 + 0.017. Luy- 
ten and Ebbighausen found a smaller value of e, 0.077 + 0.019. 

The star was found to be an eclipsing variable by S. Gaposchkin,*° 
who derived photometric elements from a light-curve based on pho- 
tographic estimates. Recently Elvey and Rudnick have observed 
the star photoelectrically at the McDonald Observatory. They have 
very kindly made their provisional results available to the writer. 
Two solutions were made, of which one, with k = 0.8, fitted the 
observations best. This solution, in which the orbital eccentricity 


was neglected, is 


asini = 3.01 X 107 a, = 18.10 
m, sin’ 1 = 32.20 b, = 16.60 
m, sin? i = 24.30 b, = 13.20 

k= 0.8 b;/a; = 0.73 
$= 64° m, = 44.30 
a; = 22.70 m, = 33.40 


It is seen that Gaposchkin’s elements are in good general agree- 
ment with those derived by Elvey and Rudnick. The range being 
0.4 mag., the star is photometrically a much better case than AO 
Cassiopeiae. 

Dr. Morgan has very kindly examined the Yerkes plates of this 
star which were measured and discussed by Luyten and Ebbig- 
hausen, and he confirms the presence of the extremely faint com- 

53 Pub. Dom. Obs. Ottawa, 4, 115, 1917. 


54 Pub. Dom. Ap. Obs. Victoria, 6, 50, 1932. 
55 Ap. J., 83, 246, 1935. 5° Bull. Harvard Obs., No. 902, 1936. 
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panion but finds the spectral class of the primary later than S Mono- 
cerotis (=O7 or O7.5) and perhaps a trace earlier than g Camelo- 
pardalis (Og). Morgan classifies it as Og—O8. In the Victoria cata- 
logue of O and B stars, this star is given as Og; Pearce classifies it as 
O7. We shall adopt O8.5. 

From Elvey and Rudnick’s figures we find a,/a = 0.471, a./a = 
0.376. The average of ai/a is 0.424; if we average according to 
brightness, we find 0.441. The theoretical values of b;/a; are, then, 
0.877 and 0.863, respectively, compared with the empirical 0.73. 
We shall, as before, adjust the photometric solution by leaving 
(3b, + a,)/4 the same and by using the theoretical ellipticity con- 
stant. We then find a,/a = 0.400 and a,/a = 0.319. Gaposchkin 
found a,/a = 0.44 and a,/a = 0.35; his ellipticity constant 0.9 is 
equal to the theoretical value. In view of the fact that Elvey and 
Rudnick’s results are still provisional, we shall adopt the average 
values: a,/a = 0.415 and a,/a = 0.33. For the inclination Gaposch- 
kin finds 68°3; Elvey and Rudnick, 64°. We adopt 65°. 

We have not yet used the condition that the total brightnesses of 
the components should correspond to the relative intensities of the 
spectra. It would seem that Am(pg) is about 1.0 mag. Taking the 
small difference in surface brightness into account, we find that the 
ratio a,/a, should be about 2/3. Leaving a, + a, unchanged, so 
that the representation of the light-curve will be little affected, we 
find a,/a = 0.45; a,/a = 0.30; hence R,/a = 0.42 and R,/a = 0.28. 

Following the procedure of applying the correction for reflection 
used for AO Cassiopeiae, we have: correction to log K,, 0.002, to 
log K., 0.010; K, = 218; K, = 294; log a/© = 1.686; log R, = 
1.309, log R, = 1.133; furthermore, allowing for the small differ- 
ence in surface brightness, log L, = 5.84, log L. = 5.39; Miow) = 
—10.0; Myo) = —8.9; log m, = 1.66, log m, = 1.53; m, = 
45:70, Mz = 33.90. 

The striking similarity between 29 Canis Majoris and AO Cassio- 
peiae has been already pointed out by Gaposchkin. But in 29 Canis 
Majoris the components appear to be separated; this fact and the 
deeper minima make this star of greater importance for the mass- 
luminosity relation. 
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In conclusion, the data on spectroscopic binaries thus far dis- 
cussed are collected in Table 12. 

The data are included in the mass-luminosity diagram (Fig. 1) 
and the mass-radius diagram (Fig. 2). The two binaries AO Cassi- 
opeiae and 29 Canis Majoris stand out by their great luminosity; 
Y Cygni appears to be an intermediate case. 


TABLE 12 


SOME SELECTED SPECTROSCOPIC BINARIES 








Star log m log L log R Mpol 
Castor C, —0. 201 —1.16 —o.18 + 7.53 
C.. —0.247 —1.2 —0.22 + 7.73 
6 Aur A +0.378 +1.83 +0.43 + 0.04 
B +0. 370 +1 .83 +0.43 + 0.04 
1 Sco AB +1.094 +3.35 +0.73 — 3.76 
V Pup AB +1.265 +3.86 +0.83 — 5.04 
Y Cyg A +1.240 +4.51 +0.77 — 6.6 
B +1.235 +4.51 +0.77 — 6.6 

AO Cas A +1.634 +5.97 +1.36 —10.3 
B +1.582 +5.58 +1.2: — 9.3 

29 CMa A +1.66 +5.84 +1.31 —10.0 
B +1.53 +5.390 +1.13 — 8.9 

















TRUMPLER’S MASSIVE STARS 

The method, used by Trumpler, of determining the masses of the 
very luminous O and B stars in galactic clusters from the Einstein 
gravitational red shift has, at least for main-sequence stars, extended 
the study of stellar masses to near the upper limit. For the most lu- 
minous of these objects are not appreciably exceeded in brightness 
(and probably not in mass) by any other known stars of classes O 
and B. Trumpler’s method is, of course, capable of being extended 
to visual binaries, although the fact that only one comparjson star 
can be used will be a disadvantage. 

Table 13 gives the observational data for seven stars taken from 
Trumpler’s paper.s? The spectrum of the third star is classed O7 by 
Trumpler, O7 at Victoria, and O7.5 by E. G. Williams. We have 
used the average value. The temperatures and bolometric correc- 


57 Pub, A.S.P., 37, 249, 1935. 
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tions are taken from Table 3 of the preceding paper; the color index 
C was supposed to be —o.2 for all the stars. Log L follows from M,,, 
and M,,.\(sun) = 4.62; log R from log L, log T., and log T.(sun) = 
3-757. It is seen that log m is roughly twice as large as log R; this 
shows that the surface gravity is about the same as that of the sun, 
as was supposed in adopting the temperature scale for log g = 4.4. 
The last column gives the order of reliability of the measured red 
shifts, according to a private communication by Dr. Trumpler. 


TABLE 13 


TRUMPLER’S STARS 






































| | } a 

| ee eee Se eS 
Star I log R | Mpg Sp. log T,| B.C.—C| Myo) | log L | log R| log mi QO 

sea ak ee te ae ae | | ee 
NGC 2244, 15} 7.9] 1 oe a 4{06 4.8:| —4.7 |-—9 | 5 49)0 66: [1 76: | 3 
; : 8} 8.5] 1.13] 4 1/09 4.50 3.0 7-1) 4 69/0.86 |1.99 | 2 
NGC 2264, 60] 14.6] 1.36] 4.7/07(+)] 4.68 4.0 8.71 5 330 82 {2.18 | 1 
NGC 2362, 1] 9.9] 1.19] 6.4/08.5 | 4.55] 3-3] 9-7] 5-73[1-28 |2.47 | 4 
NGC 6871, 2] 9.4] 1-77] 5 7|Oow 4.50 2.0 8 { 5 3et 2S 12.381 6 
Me ee >» S] 25.2] 1238) «§.:41Bo 4.40 25 PG, 1S ort 39 12:60] 5 
NGC 7380, 1] 5 4 0.93/—4 5 hk 4.50] —3.0 |—7 ° 4 " 96 |1.89 | 7 





The data of the last columns are shown graphically in Figures 1 
and 2. The general position of Trumpler’s stars in these diagrams 
is explained by the observational selection for large values of m/R, 
in connection with an intrinsic scatter of the O and B stars in these 
diagrams. 

The spectroscopic binaries Plaskett’s star and HD 698, studied by 
Pearce, for which m, sin’ 7 = 76 © and 1130, respectively, are 
probably objects similar to Trumpler’s cluster stars, whereas AO 
Cas and 29 CMa are very luminous but of only medium large mass. 
The latter stars form a continuation of the less massive stars, both 
in stellar model (Chandrasekhar) and in the mass-luminosity rela- 
tion, which for these stars is approximately 


L = m3} (4 <<m <40). (4) 


Considering the smaller masses only, the relation 


wm 
— 


L=mt (3 <m<2}) 
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fits the data considerably better. This conclusion is consistent with 

the theoretical mass-luminosity relation 
m5) (uB)73 

L=¢.—>— 6 

RI (6) 





in connection with the theoretical result 8 2 1, with the theoretical 
variation of x, and the empirical variation of R, if we assume yu to 
be constant. 

Formula (4) is only roughly correct, and its interpretation is less 
direct since the decrease of 8 with increasing m has also to be taken 
into account. 


YERKES OBSERVATORY 
June 1938 
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ABSTRACT 
“+0 eWdy 


; 1 
A table of the integral — { we 
cai a?-+-(v—y)? 


fABLES FACILITATING THE CALCULATION OF 
LINE ABSORPTION COEFFICIENTS 






is calculated by numerical quadratures 


or values of a up to 0.2, and, with somewhat less accuracy, for values of a up to 0.5. 


The coefficient of line absorption resulting as a consequence of 
natural damping, collision damping, and Doppler broadening is giv- 


en by the well-known expression’ 


j 0°" e“dy 
Kk J-w a+ (v— y)? 


The following notation is used: 


x, Line absorption coefficient at frequency v 
xy Fictitious absorption coefficient at the center of the line, 
when the damping is vanishing [eq. (8)] 


lokT i 





b™ Doppler width 6 = 


7 N Am, 


A Atomic weight in units of mo 
’ Effective natural line width [eq. (9)] 


v Deviation from the center of the line in units of the 


ee ay 


“ar 


Doppler width, i.e., » = 


a Ratio of effective natural line width to Doppler width, 


6’ 


i.e., d= —. We use the same definition of a as do 


be) 


Mitchell and Zemansky. Other authors denote by a a 


quantity equal to twice our a 


— , in frequency units 
: 


(1) 


The integral in (1) can be evaluated with the aid of series expres- 
sions or by numerical quadratures. For a2 0.5 convenient series ex- 


1 E.g., M. Born, Optik, § 93, Berlin, 1932; A. C. G. Mitchell and M. W. Zemansky, 
Resonance Radiation and Excited Atoms, pp. 319 ff., Cambridge, 1934; A. Unsdéld, 


Physik der Sternatmos pharen, § 44, Berlin, 1938. 
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pressions have been given by Zemansky,? who also gives tables of the 
integral for a = 0.5, 1.0, and 1.5. Born? gives series expressions for 
large values of a, and tables for a = 0.5, 1.0, 2.0, and 10.0. Mitchell 
and Zemansky‘ give an expression for small values of a, 0.01 or less, 
and convenient tables for the numerical calculation of the integral 
in this case. Thus, the intermediate region of a-values between 0.01 
and o.5 is not so well covered. At the suggestion of Professor B. 
Strémgren I have, with the aid of numerical quadratures, calculated 
a table of the integral in (1) as a function of v and a, for values of a 
up to a = 0.20. Larger values of a will not normally occur in astro- 
physical applications. In the calculations an accuracy in k,/Kk, of 
one-tenth of 1 per cent, or better, was maintained. 

Two different procedures were used in calculating the integral in 
(1). For moderate values of v (v < 3.5) the following procedure was 
used. The transformation of the integral in (1) to the following form, 
first given by Reiche,’ is well known: 


2, 
K I 
—= — 650 S14 CoS (vx)dx . (2) 
Ko V rJo 


For small values of a it is natural to expand the factor e~* in the 
integrand in powers of a. Retaining only the first two terms in the 
expansion, the expression referred to results 


—1 + a0  (reras| ; (3) 


Mitchell and Zemansky’ give tables with argument v of e~” and of 
the function in brackets. 
The following equation is valid quite generally: 


: , : 
—I + 206 fi ewdx | 
: 0 | 


+ Sf [e~*— 1 + axje~*/4 cos (vx)dx | 
0 ) 


Vir 


Ky om ¢ 92 ae 2a 


Ko Vr 





—— en + = 


Ko Vv WT 





(4) 


2 Phys. Rev., 36, 219, 1930; cf. also Mitchell and Zemansky, oP. cit., p. 329. 

3 Op. cit. Born’s 7 is 1/a, his x is v/a, and ky/Ko is —1/ (a7) times Born’s y. 

4 Loc. cit. Cf. also D. H. Menzel, Pub. Lick Obs., 17, 229, 1931; E. F. M. v. d. Held, 
Zs. f. Phys., 70, 508, 1931. 

5 Verh. d. D. Phys. Ges., 15, 3, 1913. 
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or 


=e + a+ oft) + ¥(a, 2). (5) 

For relatively small values of a this equation is a more convenient 
basis for numerical quadratures than (2), since the product of the 
first and the second factor is then always quite small. With increas- 
ing values of v, however, the calculations become increasingly labori- 
ous, because the largest interval that can be conveniently used in the 
numerical quadratures is inversely proportional to v. Equation (4) 
was used up to v S 3.5, and for a S o.10 up tov S 4.0. 

For moderately large, and large, values of v, equation (1) is quite 
a convenient basis for the numerical quadratures. Generally, for the 
relatively small values of a in question the denominator of the in- 
tegrand in (1) varies in such a way in the neighborhood of y = v as 
to make the equation unsuitable as a basis for numerical quadra- 
tures. When, however, v is large, the numerator is so small for y = v 
as to make the contribution of the immediate neighborhood of y = v 
to the integral practically negligible. 

The two procedures of calculation thus supplement each other. A 
few values were calculated by both procedures so that a check on the 
calculations was obtained. 

For a S$ 0.2 and v 2 5 the integral in (1) is independent of a 
within an accuracy of better than two-tenths of 1 per cent [eqs. (6) 
and (7)]. Thus, a table with one argument, equivalent to the table 
given by Mitchell and Zemansky‘ of the bracketed function in (3), 
suffices for v 2 5. 

The following asymptotic series expression [eq. (4)] was used to 
check the calculations for moderately large values of v (v 2 4): 


co 
Jf [e~* — 1 + axle~*/4 cos (vx) dx 
V r/o 
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The results of the calculations are given in the tables that follow. 
It was not necessary to carry out numerical quadratures for all the 
values given, interpolation being used as far as the resulting inter- 
polated values were deemed sufficiently accurate. 

A table of less accurate values is also given for a up to 0.5. 
The accuracy of the calculation of these values was about 1 per 
cent. The table may be of use in certain problems of theoretical 
physics. 

The arrangement of the tables is as follows. Table 1 contains nu- 
merical values of the third term in (4) and (5) —i.e., ¥(a, v)—togeth- 





























TABLE 1 
| a=0.0o0!I a=0.05 a=o0.10 a=0.20 
| | | 

Ky/k | 1ooy (a, v) Ky/Ko | roo (a, v) Ky/ko | 1oow(a,v)| Kyp/Ko rooy (a, v) 
0.00 lo.99 |+0.01 0.95 Ito 2 ' 90 Ito 93 |o.81 |+3.5 
0.25 | 93 i+ .o1 89 |+ .20 .85 + yr 99 |\+2.9 
0.50 | .77 | 00 75 |+ 10 72 |+ .38 66 |+1.4 
0.75. ey, fole) . 56 |— .O1 55 |— .O4 Sd mT 
DO. 37 .0O a7. |= .8 37 \- 32 3e (Ra 
Los Fei fore) 12 |— og 13 \— 35 1 |—1.3 
2:6. O21 000 030 |— .031 040 |— .13 059 |—0.52 
ee 0032 . 0000 0084] — 0058) O15 |— .02 O27 i— .12 
a O:. OoogI .00000} .0040)/— .0008} .0079/— .0036} .o16 |— .o21 
2.5 00054 .00000| .0027/— .0001| .0053|— .0008] .o11 |— .co60 
LO, .00039 .00000} .0020/— .0001] .0039/— .0003| .0078/— .0028 
Btn. 00030 . 00000 0015| 0000} .0030/— .0002} .0060/— .0017 
SO... 0.00024} Oo ee OO012| 0.0000/0.0024|/—0.0001/0.0048]—0. 0010 














er with two-figure values of x,/ko. The purpose of Table 1 is to 
show the relative importance of the third term for various values of 
aandv. For a = 0.03 the third term is always less than 0.5 per cent, 
for a = o.10 it is less than 3 per cent, and for a = 0.20 it is less than 
g per cent. Thus, our calculations show that when an accuracy of, 
say, 10 per cent is sufficient, the expression e~” + a-¢ (v) can be used 
up to a S 0.2, i.e., considerably farther than was anticipated.‘ Ta- 
ble 2 gives numerical values of k,/x. with arguments @ and v up 


tov = 5. Table 3 gives (1/a) (x,/x.) with argument v for v 2 5 
Table 4 contains two-figure values of k,/k, up to a = 0.5 for v S 5. 


Forv = 5 Table 3 gives two-figure accuracy up to a = 0.5. 
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With the aid of the tables, the calculation of the line absorption 
coefficient within an absorption line originating from the ground 
state runs as follows. Let 7 be the mean lifetime corresponding to the 


TABLE 3 






































Ky/ko ky/Ko Ky ‘Ko 
D B v ns v aw 
MD cig a aid eos 0.0241 10.0 0.00573 15.0 0.00251 
ere... 0197 10.5 00519 15.5 .00235 
6.0.. o164 II.0 00472 16.0 .00221 
6.5 0139 4 .00432 16.5 00208 
F30. OIIg 12.0 .00390 17.0 00190 
os. .O103 2.5 .00305 175 00185 
B20; ; .00903 13.0 .00337 18.0 .OO175 
8.5 .00798 ra .5 .00312 18.5 00166 
9.0 .00710 14.0 .00290 19.0 .00157 
Ors. ; 00636 14.5 .00270 19.5 OO149 
fo) one 0.00573 I5.0 0.00251 20.0 0.00142 
TABLE 4 
Ky lke 
‘ ™ 
a=0.0 a=o.1 a=0.2 a=0.3 a=0.4 a=0.5 
0.00.. I.00 ©.90 0.81 0.74 0.67 0.62 
0.25. 0.94 85 Be .70 .64 .59 
0.50. 0.78 .72 .66 .62 a | 33 
0.75. 0.57 -54 -52 -50 47 -44 
1.00. 0.37 37 37 -37 30 36 
285 0.21 et Py i . 26 (a7 27 
G50: '. 0.10 03 .16 a 4 .19 . 20 
4 ee 0.047 074 .096 .12 Pe é¥4 
2.00... 0.018 .040 .060 .076 .09O0 103 
2256. 0.002 .O15 .027 .038 .049 .058 
3.00 0.0001 .0079 .o16 .023 .030 .037 
3.50. ©.0000 0053 OI .o16 .O21 026 
4.00. ©.0000 0039 .0078 .O12 .o16 O19 
4.50. 0.0000 0030 0060 . 0090 .O12 .O15 
5.00. 0.0000 0.0024 0.0048 ©.0072 0.0096 0.012 























transition between the two stationary states k (upper state) and 1 
(lower state) in question. Then 











Ig. I . Ig »-M Am ’ 
“Bacal ‘ N, = o.N ? 8 
" 7 qx 84/7 5) ; (8) 


7 Tg: 84v/m NRT 
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where g; and q, are the statistical weights of the states k and 1, c the 
velocity of light, b” the Doppler width in frequency units, and N, 
the number of absorbing atoms in state 1 per gram. Further, 


 . I 
= — 2 er + — S, (9) 
ey 21 


where S is the number of broadening collisions® per second, and 


pm = j2kT d | 
NAm, Cc | 
os (10) 


po) = lokT 1| 
NAm, dr} 





b™ and b® being the Doppler width in wave-length and frequency 
units, respectively. Then if \X—\, is the deviation of the wave length 
considered from the center of the line, in wave-length units, 


.* a (11) 
and, further, 
6’ 
a= 75° (12) 


Now, Table 2 gives x,/x, as a function of a and 2, then x, follows 
immediately, since x is given by (8). In making calculations of line 
contours, or equivalent widths, the values of X—X, can usually be 
chosen such that interpolation in Table 2 with respect to v is avoided. 

For a subordinate line the calculation runs quite similarly if it 
is permissible to assume that the atoms in the lower state are dis- 
tributed according to the normal distribution. In that case it is only 
necessary to use the well-known Weisskopf-Wigner equation corre- 
sponding to equation (9). 

OBSERVATORY, COPENHAGEN 


June 1938 


6 E.g., Mitchell and Zemansky, /oc. cit.; Unséld, Joc. cit., and the references there 
given. 








ON THE INTEGRATION OF THE EQUATION OF 
RADIATIVE TRANSFER 


P. SWINGS AND L. DOR 


ABSTRACT 

The equation of the radiative transfer has been integrated by a generalization of 
Spitzer’s method, which makes use of Bessel functions. 

1. In solving the equation of radiative transfer appropriate to the 
problem of the formation of absorption lines, it is necessary to take 
into account the variation of the Planck intensity B, as function of 
the total optical depth within the line. To a sufficient approxi- 
mation, B, may be expected to be a linear function of the form 
a, + b,r, where 7 is the optical depth in the continuous spectrum; 
7 is defined with respect to the Rosseland mean absorption co- 
efficient over the whole spectrum. On the other hand, we can- 
not, in general, expect B, to be a linear function of the optical depth 
i,, The point at issue here is the circumstance to which Strémgren! 
has recently drawn attention: 

Line absorption affects the intensity within the radiated spectrum in two 
ways. First, it cuts down the intensity as would increased continuous absorption, 
the effect depending upon the existence of a temperature gradient in the stellar 
atmosphere. Second, it causes a deviation from Kirchhofi’s law, in the sense 
that emission is reduced, so that there is a further reduction in the radiated 
intensity. 

Spitzer’s recent solution? of the equation of transfer in terms of 
Bessel functions is seen to take into account the second of the fore- 
going two effects ‘‘rigorously,’”’ while the first effect is taken into ac- 
count by assuming for B, a linear function in ¢, (which assumption 
cannot generally be valid). For this reason the perturbation theory, 
developed by Strémgren,' is of greater value inasmuch as Strém- 
gren’s method in fact allows quite a considerable variation in 7,.3 

However, it is of some interest to supplement Spitzer’s analysis 


t Ad. J., 86, 1, 1937. 
2 Ap. J., 87, 1, 1938. 3 Ap. J., 86, 1, 1937, Table IT. 
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by calling attention to the fact that several other possibilities of 
rigorous solutions exist for the equation of radiative transfer. Such 
possibilities are indicated here in a summarized form, no analytical 
discussion or expression of the residual intensity being added, as this 
would proceed along the same line as in Spitzer’s paper. In any spe- 
cial example (i.e., for a definite model and frequency), it will be possi- 
ble to select the analytical case which best approximates the physical 
conaitions. 

2. Using classical notations, the equation of radiative transfer 
may be written 


S(t) _ , a +en) _i1te2 


\ J, (t 





= S(t, B,(t,) }, 
dt? A I + n, v\"y I -}- en, v y) | (1) 
where 
dt, = (1 + »,)dr,. (2) 
We assume, as usual, 
1 + €n,0 
pi es - 
I + en, E y(t,) a, + p,t, ? (3) 
and substitute 
y(t,) ae J, (t,) as (a, -+- p,t,) (4) 
‘ I+ 9, 
Mi) <= S ; 
by) i+. (5) 
Thus, we get instead of (1) 
ay _ 
ae = 3Mb) > yt) - (6) 


Spitzer’s treatment consists essentially in finding adequate sub- 
stitutions which convert equation (6) into a Bessel equation of the 


dy i, y dy = (: ee =) y=0. (7) 


dw w dw 


type 
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Solution I of Spitzer is 


1 = 4 (1 + mt,)~*/*+ , 


with 
m = (s+1)— 


s, A, and D being constants. 
Solution IT is 


j= L+ Me, 


L, M, and u being constants. 
3. Substituting 


z=a+t pt, (8) 
(6) keeps its essential form and becomes 


dy, _ 3A(3) 


dz? 8B 2 





* yx(z) . (9) 


We may try to find a substitution of independent variable w(z) and 
of unknown 
y(z) = u(z) + y:(2) , (10) 


which transforms (7) into (9). 
After some easy algebra, it is found that we need, therefore, 


se ! 4+(2 a) (Se) =° ies 


duw= w dw u dz dw 


ai? Ry 8 nt a (os rdz\_(,4™\{_ 
Bu at) E dz dw’ wdw tS ; = 0, (42) 











du 





Solution I of Spitzer corresponds to 


z=aw; w=2V?; g=1+ mt,, 
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a, 6, and m being constants which may be immediately related to 
Spitzer’s coefficients. 
Solution II corresponds to 


w=ae*; “= const.; s=t. 
The equations (11) and (12) generalize Spitzer’s results; (11) may 
be replaced immediately by a first-order linear equation which 


brings 


dz I 
ana f( 2. du), de ; (13) 
const. + —- . 

“u dzZ w 


We may take, for instance, 





dw 








W 





I 2(a+b s/2 
u=e*; z=at+b,; w= F » @Be?\ i" : 
then \ has the following form 
b?B? nae 
3\ = & + bn? Be¥ + a le ae, (14) 


a, b, A, B, and n, being constants. and the solution is, accordingly, 
y= Cres [cx : T,(w) + G° K,(w)] . (15) 


T,(w) and K,(w) are the usual Bessel functions of imaginary argu- 
ment. 

Obviously, instead of transforming (6) into (7), we could have 
taken the other Bessel equation, or, more generally, any other well- 
known type of differential equation of the second order, namely, that 
providing the orthogonal polynomials. 

4. Equation (6) being linear and homogeneous with regard to y 
and d’y/dt, the substitution 


y = eS udt (16) 
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will transform (6) into a Riccati equation 


du ees 
4 + wv = 3r(t), 


and a solution by quadratures is always obtainable when 


dp 
aX) = = + #’, 
dt 
uw being any function of t. 
The solution is actually 
I 
“u=pe-—-, 


with 
y= eS udt a ferS uatdt ' 


For instance, if 


_fO rt) - LY 
3X(Z) Jo’ 


= 1® or 


f(t) being any function, we have 


fl) +e dt 
y = af(4) + of(t) \ . 
The simple case f(t) = Al” brings 
m(m — 1) 7 
3X(t) a AE and y= c,t™ + cf. 


f? 


5. Let 


y = fl 


(19) 


be a solution of (6); suppose we may write ¢ = F(y). Then (6) be- 


comes 


= = 3\[F(y)] + y = o(y) + y. 


(24) 
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Thus 


dt = [2f{e(y)- y+ dy +a} + dy. (25) 


By quadrature we shall find ¢/ as a function of y. 
In other words, if we take any function ¢(y), we shall get y = f(t) 
by quadrature, and then 


3A(t) = o(y) = gLf(d)] . (26) 
For instance, if 
gly) =—, 


we find 


with 
2A 


en At? + 2Bt+ 2C’ 





A, B, and C being constants. 

6. Summarizing, we may state: (1) Spitzer’s method introducing 
Bessel functions may be generalized by formulae (11) and (12) [Ex- 
ample: formula (15)]; (2) and substitution (16) give a solution by 
quadrature when X(¢) verifies the general relation (17) [Examples: 
formulae (20) and (21)]; (3) another integration by quadrature is 
given by formula (26); this may in some cases give simple solu- 


tions. 


Our thanks are due to S. Chandrasekhar for valuable suggestions. 


UNIVERSITY OF LIEGE, BELGIUM 
DEPARTMENT OF ASTROPHYSICS 
April ro, 1938 











NOTES 


ON COLLAPSED NEUTRON STARS 
ABSTRACT 


Some consequences are discussed of the hypothesis that certain stars and cores of 
stars are composed mainly of neutrons. On the assumption that supernovae represent 
rapid transitions of ordinary stars into neutron stars, the large red shifts observed in 
the spectra of the recent bright supernovae are interpreted as gravitational red shifts. 
The neutron-star hypothesis, in conjunction with the general theory of relativity, leads 
to a theory of critical stellar masses. 

In view of the rapid advances made recently in the observation of 
supernovae, it seems appropriate to give here a brief summary of 
some conclusions of a series of theoretical investigations concerning 
the properties of highly collapsed neutron stars, which it is hoped 
may be ready for publication in the near future. 

The neutron-star hypothesis was first introduced by Baade and 
myself’ in an effort to account for the tremendous liberation of 
energy in supernovae. We suggested the existence of stars and cores 
of stars the average density of which is comparable with the density 
of matter ordinarily encountered only inside of atomic nuclei. With 
the designation ‘‘neutron star’ we do not wish to imply, however, 
that such a star is to be regarded as a giant nucleus composed of 
separate neutrons of precisely the same character as free neutrons. 
We only suggest, that, in contradistinction to ordinary stellar mat- 
ter, in neutron stars already minute regions whose linear dimen- 
sions are larger than 6 = e?/m.c? = 2.8 X 10% cm are electrically 
neutral. 

We here briefly describe some of the properties of neutron stars, as 
well as some new observations of supernovae which tend to support 
the neutron-star hypothesis. 

a) Cold neutron stars, according to present knowledge, represent 
states of lowest energy that matter may assume without being com- 
pletely transformed into radiation. 

1 W. Baade and F. Zwicky, Proc. Nat. Acad. Sci., 20, 259, 1934, and Phys. Rev. 45, 
138, 1934, and 46, 67, 1934; F. Zwicky, Scientific Monthly, 40, 461, 1935. More recently 
the neutron-star hypothesis has also been considered by G. Gamow, Atomic Nuclei, 
p. 234, Oxford, 1937; L. Landau, Nature, 141, 334, 1938; and others. 
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b) According to the general theory of relativity, a limiting mass 
of stars exists for every given average density (Schwarzschild limit).? 
At this limit the energy liberated because of gravitational packing is 


E= (: 7 =} Mc? = 0.58 Mc’, 
3m 


where M is the proper mass of the star. For an average density p = 
10° gm/cm3 the limiting mass is M; = 6.4 X 10*4 gm. The deriva- 
tion of these results which was obtained in discussion with Professor 
R. C. Tolman will be communicated in a joint paper with Professor 
Tolman. 

A star which has reached the Schwarzschild limiting configuration 
must be regarded as an object between which and the rest of the 
universe practically no physical communication is possible. For in- 
stance, the velocity of light on such a star is infinitely small, so that 
it requires light from this star an infinitely long time to reach any 
external point. Also, the gravitational red shift is complete in the 
sense that light originating on the star arrives at any external point 
with the energy zero. It is, therefore, impossible to observe physical 
conditions in stellar bodies which have reached the Schwarzschild 
limit. It should, however, be possible to observe stellar bodies in 
stages intermediate between the ordinary configurations and the col- 
lapsed configurations of limiting mass just described, provided that 
such are accessible. 

c) We may express the limiting mass M, of collapsed neutron 
stars in terms of the universal gravitational constant 


I’ = 6.66 X 10°° gm™ cm; sec’, 


the charge of the electron e and the masses m,., m,, and m, of the 
electron, the proton, and the neutron, respectively. We obtain 


M, = aR3/?m,,, (1) 
where 
e? 
R = —— = 2.3 X 10 (2) 
I'm,m, 


2 See, e.g., R. C. Tolman, Relativity, Thermodynamics and Cosmology, p. 247, Oxford, 
1934. 
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is the ratio of the electrical to the gravitational attraction between an 
electron and a proton and where a is a dimensionless number of 
the order 1. The number of neutrons in the mass M, is N =akR?3”. 
According to an idea first discussed by H. Weyl,’ the ratio R may 
play a role in the determination of large numbers, such as the total 
number of particles in the universe and the ratio between the ‘‘ra- 
dius” of the universe and the radius of the proton. The fact that 
M,, is proportional to R*/? suggests possible observational tests of 
Weyl’s idea, which will be discussed in another place. 

d) If supernovae are transitions from ordinary stars into neutron 
stars, the observation of light-curves and spectra of supernovae 
should furnish us with direct evidence of the neutron-star hypothesis. 
For instance, the surface of the central star of a supernova should be 
exceedingly hot, the acceleration of gravity very high, and light 
coming from this surface should be subject to enormous gravitational 
red shifts. Now, it is a significant fact, first observed by Dr. R. 
Minkowski in the spectrum of the recent bright supernova in IC 
4182, which was at its maximum brightness in August, 1937, that all 
of the permanent features of this spectrum have gradually shifted 
toward the red, until in June, 1938, a stage S was reached when this 
shift amounted to roo angstroms.* On our hypothesis we may tenta- 
tively interpret this red shift as a gravitational red shift. Assuming 
that the central star has a mass M, = 2 X 10% gm, equal to that 
of the sun, the following characteristics for the central star in the 
stage S may be derived: radius r = 74 km; average density 
p =1.2 X 10% gm/cm’. Since the supernova in the stage S was 
about one million times brighter than the sun, it follows that the 
effective surface temperature 7 > 1.8 X 107 degrees absolute. For 
the limiting mass M; = 6.3 X 10*4 gm we obtain r = 3700 km, p = 
1.2 X 10? gm/cm3, and T > 3.1 X 10° degrees. The neutron-star 
hypothesis in conjunction with observations on supernovae may 
therefore lead to some new and far-reaching tests of the general 

3H. Weyl, Raum, Zeit, Materie, 5th ed., p. 277, Springer, Berlin, 1923. See also 
Naturwissenschaften, 22, 145, 1934. 

4 The spectra of the supernovae in NGC 1003 (1937) and NGC 4273 (1936) exhibit 
similar effects. I am indebted to Dr. Minkowski for the permission to make use here 
of these unpublished observational facts, which he will present in detail in a paper 
shortly to be published in this Journal. 
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theory of relativity. Also, the fascinating problem now presents it- 
self of investigating how certain well-known physical processes, such 
as nuclear reactions, will be modified when they take place inside of 
highly collapsed stars in which the very properties of time and space 
are drastically altered. 
F. ZwIcky 
CALIFORNIA INSTITUTE OF TECHNOLOGY 
PASADENA, CALIFORNIA 
August 8, 1938 


SPECTRAL TYPES AND RADIOMETRIC OBSERVA- 
TIONS OF STARS OF LARGE INFRARED INDEX 


In a recent paper Hetzler' announced the discovery of a number 
of red stars with exceptionally large infrared indices, some of the 
values being as large as 10 mag. If one uses these indices as a basis 
for temperature determinations, on the assumption of black-body 
radiation, the large indices will yield temperatures of about 1000”. 
The objects of infrared magnitudes 4 to 6 having this temperature 
would be within easy reach of a vacuum thermocouple, the radio- 
metric magnitudes being of the order of o to —2. For this reason it 
was considered worth while to study the stars radiometrically. At 
the same time a determination of their spectral types seemed de- 
sirable. 

The radiometric measurements were made with a vacuum thermo- 
couple? at the 24-inch reflector with a junction of bismuth against 
95 per cent bismuth plus 5 per cent tin and with a highly sensitive 
galvanometer. Provisional results were obtained with the equip- 
ment not yet developed to its full sensitivity; at the time of the 
measurements a star of radiometric magnitude 3 gave a measurable 
deflection. The measurements consisted of a comparison of three of 
the more promising stars of large infrared index with a Orionis and 
a Tauri. The visual magnitudes of the chosen stars were 11.5, 12.5, 
and 14.0; the first two stars are included in Table 1. No deflection 


tA p. J., 86, 509, 1937. 


2 A more complete description of the radiometric equipment will be given later. 
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of the galvanometer was observable in the case of the selected stars, 
indicating that their radiometric magnitudes must be fainter than 
3.0. 

The spectra were obtained at the 4o-inch refractor with the one- 
prism spectrograph designed by Dr. Kuiper.’ Twenty stars have 
been observed, their infrared indices ranging from 2.0 to 8.0 mag. 


TABLE 1 








1855 
INFRARED 
STAR as ae cas INDEX SpectruMmf 
R.A. —) (HETZLER) 

Se re 19h 26™1 +23° 20/3 8™o0 S3 
V:1 r-27:20 54 23.0 ae Mo 
Wer. I 51.0 57 54.6 6.8 M7 
Se. 22 5:4 Bor Se 6.6 M7 
Vcc... 19 18.5 5 aan ae 6.5 (M8) 
IV:6. 19 17.0 26 52.8 6.2 S3 
V:5§.. I 32.4 55 40.7 6.1 (M7) 
ee 21 47.6 tee fy, 6.0 M7 
V7. EO: £7...2 26 42.6 5.9 M8 
bons... | 2m 40.4 Ce ee 5.6 M7 
FI Lyr. 18 36.4 28 49.0 5.5 Mg 
it:3* 4 54.0 56 25.6 soy M6 
eae ee I 24.2 57 0.9 5.0 R or Nt 
pe C3461 SS 4.5 4.9 M6 
V313.. I 30.2 58 58.5 4.5 Sec 
W206... I 34.6 55 36.8 4.2 Ms5 
RV Cam 4 18.7 ae 4.0 M6 
BD+53°413 oe I 44.3 53 31.6 3.8 M1 
Vers... or 16-0 58 4.4 205 Ms 
BD+22°3792.... 19 41.9 22 24.6 3.2 K4+ 
Waclas..... O 23.2 <a oo a 5% M2 
Per... 2 9.0 Be 29-2 2:0 K1 
EP Lyr.. fo 12.6 BREEN. ONG Gua feats che hevs F8 
ma Lyr,:. IQ 7.2 +26 32.0 M5 

















t The S and M stars are classified very easily. The values in parentheses are uncertain and provisional 


t In the spectral region used, these spectral types look very similar. 


Table 1 gives the spectral types obtained. The stars are designated 
by their field and chart numbers as given by Hetzler, the Roman 
numerals representing the field, the Arabic numerals the chart. 
Columns 2 and 3 give the position for 1855, column 4 the infrared 
index, and the last column the spectral type. The stars marked with 
an asterisk were observed radiometrically. 


3 Ap. J., 87, 592, 1938. 








penance _ —_—— — — ———_ ee _ — - — we 





PLATE XX 





oe *e 
+ - bd o ” e 
> 
° o e ° - e 
. e . ‘e . s . 
o © 
y ¢ 
e es @ « . m Fy € ? e . 
2 . 
o a H 
( * * ab 5 ae ° > ‘ h 
. ° . 
. . e e 3 + 
& nee e e- 
. e e 
> e a e . 8 
o ° . ° +. ° 
. . >. 
. 
ee | 
ef } 
o«| 
ee ad 
s4 Pe °s. 
v 
a ¥ : ae : d 
we 
ee 
ee 
7 « . * ad Ad 











a) A double exposure of the planetary nebula NGC 7027. The right-hand exposure was 
obtained with a focal setting at Ha, while the faint diffuse image on the left was taken with 
the focal setting to the violet of Ha. The exposures were five minutes each 

b) The spiral nebula M 33, photographed on August 22/23, 1938 (on the right), with a six 
hour exposure and with a focal setting appropriate for Ha. 

c) Was obtained with a six-hour exposure on August 21/22, with a focal setting of 
2.5 mm to the violet of Ha. The strongest emission condensation is marked by an arrow 

d) A new Wolf-Rayet star. The exposure on the right is for Ha, and that on the left is 
toward the violet of Ha. The exposure time was twenty minutes each. 


All four exposures were obtained on Agfa Superpan Press film with filter 23A. 
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Table 1 shows that this sample of red stars belongs to classes 
M and S; no new types are found among the stars observed. Accord- 
ingly, the great redness of the stars must be interpreted as being 
caused by heavy band absorption and not by excessively low tem- 
perature. The radiometric observations lead to the same conclusion. 


I am indebted to Professor G. P. Kuiper for the suggestion of the 


problem and for helpful discussions. 
Cari F. Rust 
YERKES OBSERVATORY 
October 15, 1938 


A METHOD FOR THE DETECTION OF SMALL 
OBJECTS HAVING EMISSION SPECTRA 

During the course of an investigation of the colors of knots in 
extragalactic nebulae, a method was devised to detect objects hav- 
ing strong Ha emission in their spectra. Briefly, the method is as 
follows: The suspected object is photographed through a refracting 
telescope on a red-sensitive plate behind a red filter. One exposure 
is made with the telescope focused at Ha, and another exposure is 
made on the same emulsion at a shorter wave length chosen in such 
a way that the images of normal stars on the plate are equal in size 
at the two focal settings. If the object in question has an emission 
spectrum, the image will be sharply in focus at Ha and will be con- 
siderably out of focus at the other wave-length setting, whereas 
the two images of a normal star will be of equal size. 

Plate XX shows some of the results that have been obtained by 
this method with a 6-inch refractor, at the McDonald Observatory. 
All the exposures were made on Agfa Super Pan Press film exposed 
behind a Wratten No. 23A filter. 

Plate XX a is an enlargement of two five-minute exposures on the 
planetary nebula NGC 7027, and shows the effectiveness of the 
method very clearly. For the nebula, the Ha image is quite strong, 
whereas the image taken to the violet of Ha is so weak that it has 
practically disappeared. 
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Plate XX 0 and XXc are enlargements from two six-hour exposures 
on the extragalactic nebula, M33, taken at Ha and to the violet of 
Ha, respectively. The emission character of the three brightest 
knots, first announced by Hubble,‘ is confirmed by the increased 
strength of these knots photographed in Ha light. In addition, two 
fainter knots, hitherto not observed spectroscopically, are found by 
this method to be emission nebulae. 

Plate XXd shows a photograph of a new Wolf-Rayet star dis- 
covered by this method. The position of the star is atg"19™46%7, 
5 + 29°28/4 (1900), about 4’ northwest of the fifth-magnitude star 
2 Cygni. Two objective-prism plates of this object showed Ha, H£, 
and Hy in emission. The fact that the lines Nr and N2 are absent 
from the spectrum indicates that the object is probably not a 
planetary nebula. 

This method of detecting objects having emission spectra has two 
limitations: first, it can be used only on objects of approximately 
stellar size as seen on the photographic plate; and second, it cannot 
be used with a reflecting telescope because of the absence of chro- 
matic aberration. It is expected that this method will prove useful 
for detecting and confirming suspected planetary nebulae and Wolf- 
Rayet stars which are too faint to be observed spectroscopically.’ 

CARL K. SEYFERT 
McDONALD OBSERVATORY 
October 1938 


tAp. J., 63, 236, 1926. 
2A somewhat similar method was used by Barnard (Ap. J., 14, 151, 1901) for 
visual and photographic observations of novae and planetary nebulae. 





ERRATUM 


On Plate XVI in the October issue the pointer indicating the line Ha was 
inadvertently displaced to the left. It should have been placed in line with the 
heavy emission at the extreme right of the reproductions. 











